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PROTOTYPE SUPERCONDUCTING GRAVITY GRAD I OMETER 



PREFACE 


This report describes research on the superconducting gravity gradiometer 
program performed at the University of Maryland from July 1980 to July 1985 
under NASA Contract NAS8-33822. 

This report consists of three main parts: Parts 1 and 2 describing the 

theoretical and experimental work on a prototype superconducting gravity 
gradiometer, respectively, and Part 3 discussing the design of a new advanced 
model of superconducting gradiometer. These three parts represent three 
separate papers being submitted for publication. In addition, two published 
papers, one reporting a null test of the gravitational inverse square law 
performed with the prototype superconducting gradiometer and the other dis- 
cussing space applications of the advanced three-axis instrument, are attached 
as appendices. 

An advanced three-axis superconducting gravity gradiometer, designed 
along the line proposed in Part 3 of this report, is being developed under 
new NASA Contract NAS8- 36165. An associated instrument, a six-axis super- 
conducting accelerometer, is under development with support from AFGL under 
Contract F19628-85-K-0042 . Electronic control of the gradiometer has been 
supported in part by Army Contract DACA-72-84-C-0004. 
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SUMMARY 


A sensitive and stable gravity gradiometer is needed for precision tests 
of fundamental laws of physics and for moving-base gravity survey and inertial 
guidance. Through the five years of research and development at the University 
of Maryland we have demonstrated the feasibility that superconducting tech- 
nology can be utilized not only to lower the intrinsic noise of the instrument 
but also to meet many practical challenges of operating a sensitive gravity 
measuring instrument in a noisy environment. 

A relatively simple prototype single-axis superconducting gravity gradi- 
ometer has been constructed to investigate the basic physics of a superconduct- 
ing gradiometer. At the same time, a detailed analysis of the instrument 
dynamics has been carried out including extensive error modelling. Thorough 
experimental tests of the instrument have shown that the superconducting device 

closely follows the analytical model. The performance level of 0.3 'v 0.7 
-1/2 _ -9 -2 

E Hz (1 E = 10 s ) achieved with this instrument in the laboratory 
without any active control or compensation represents the best reported sensi- 
tivity of any gradiometer to date. The instrument has already been used 
successfully to perform a new null test of the gravitational inverse square 
law. 

Based on the experience obtained with this first instrument and additional 
superconducting technologies developed to improve the performance of the super- 
conducting gradiometer, an advanced design of a three-axis superconducting 
gravity gradiometer has been produced. Incorporated into the new design are 
the concepts of a "superconducting negative spring", "three-dimensional 
residual common mode balance" and a "six-axis superconducting accelerometer." 
Various feedbacks will be applied to control the instrument and the platform. 
This second generation superconducting gravity gradiometer should be able to 
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-4 -1/2 

meet the instrument noise goal, 3 x 10 E Hz , defined for NASA's Gravity 
Gradiometer Mission (GGM) . 

Part 1 of this report discusses the theory of the superconducting gravity 
gradiometer. Although the particular superconducting circuit used for the 
prototype is analyzed here, the method and results are quite general and can 
be readily applied to the new model. Part 2 reports details of the design, 
fabrication, and various tests of the single— axis gradiometer. Calibration 
experiments performed by using gravity gradient signals are described here. 
Part 3 contains the design and theoretical analysis of the new gradiometer. 

It is found that the improved common mode rejection characteristic and the 
ease of g-nulling in the new instrument would enable its operation in a more 
hostile terrestrial moving-base environment. A null test of the inverse 
square law by use of the superconducting gradiometer is reported in Appendix 
I. Finally, the basic concept of GGM is discussed in Appendix II. 

The Maryland superconducting gravity gradiometer project has served as 
an excellent training ground for physicists. A strong Ph.D. thesis was 
produced on the development of the prototype gradiometer and the test of the 
inverse square law (Chan, 1982). Another Ph.D. thesis is nearing completion 
on the development of the six-axis accelerometer. A third student is test- 
ing a transducer concept for the cryogenic gravitational wave detector, 
which came as a spin-off from the gradiometer project. The many challenges 
that one has to overcome to realize the fabulous sensitivity of the space- 
borne gravity gradiometer, orders of magnitude beyond the state of the art, 
and the exotic new technologies being invented to meet them have been inspira- 
tional to the scientific team of the Maryland gravity gradiometer project. 
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PART 1 

THEORY OF A SUPERCONDUCTING GRAVITY GRADIOMETER 


SUPERCONDUCTING GRAVITY GRAD IOMETER FOR SENSITIVE GRAVITY MEASUREMENTS: 

I. THEORY* 


H.A. Chan and H.J. Paik 
Department of Physics and Astronomy, 
Univeristy of Maryland, College Park, MD 20742 


Due to the Equivalence Principle, a global measurement is necessary to 
distinguish gravity from acceleration of the reference frame. A gravity 
gradiometer is therefore an essential instrument needed for precision tests of 
gravity laws and for applications in gravity survey and inertial navigation. 
Superconductivity and SQUID (Superconducting QUantum Interference Device) 
technology can be used to obtain a gravity gradiometer with very high 
sensitivity and stability. A superconducting gravity gradiometer has been 
developed for a null test of the gravitational inverse square law and 
spaceborne geodesy. Here we present a complete theoretical model of this 
instrument. Starting from dynamical equations for the device, we derive 
transfer functions, common mode rejection characteristic, and an error model 
of the superconducting instrument. Since a gradiometer must detect a very 
weak differential gravity signal in the midst of large platform accelerations 
and other environmental disturbances, the scale factor and common mode 
rejection stability of the instrument is extremely important in addition to 
its immunity to temperature and electromagnetic fluctuations. We show how 
flux quantization, the Meissner effect and properties of liquid helium can be 
utilized to meet these challenges. 


*Work supported by NASA under contract No. NAS 8-33822. 
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I. INTRODUCTION 


Highly sensitive gravity sensors are needed to investigate fundamental 


properties of gravitation and to improve accuracies of gravity survey and 


inertial navigation. Torsion balances have been used over two centuries for 
sensitive gravity experiments [1], Spring-mass, pendulum, and free-fall type 
gravimeters [2] have been developed as survey instruments as well as a 
superconducting version of the first type [3], Research to develop detectors 
for gravitational waves of extraterrestrial origin has started over two 
decades ago [4], Cryogenic mass-quadrupole type and laser interferometer type 
detectors are under vigorous development [5], The last two decades have also 
seen dedicated efforts to develop room temperature gravity gradiometers [6-8] 
for moving base survey applications. Superconducting gravity gradiometers 
have emerged more recently as an outgrowth of the superconducting transducer 
work for low temperature gravitational wave detectors [9], 

The extremely weak nature of gravitational interaction poses a challenge 
to the state of the art technology for signal transduction and amplification 
as well as isolation of environmental noise. To compound the problem, the 


gravitational field cannot be distinguished in a local measurement from 


acceleration of the reference frame by the Equivalence Principle. In order to 
separate gravity from frame accelerations, one must resort to a second order 
measurement using the tensor nature of gravitational field gradient or 
gravity gradient". When the platform is undergoing a linear acceleration, a 
differential measurement over a baseline between two proof masses will cancel 
out the acceleration noise, leaving gravity to be detected as the signal. 
Likewise, an angular acceleration can be taken out by combining signals from 
four proof masses as we will see in Section II. Thus, unlike in 
electromagnetism, where a single test charge can be used to determine the 
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field uniquely, a tidal force sensor or a "gradioraeter" is the fundamental 
instrument in gravity which is capable of measuring its field, independent of 
platform motion. True acceleration measurement, In turn, requires removal of 
gravity noise which again calls for the use of a gravity gradioraeter. It is 
therefore not surprising to find that instruments employed in most precision 
gravity experiments, such as torsion balances and Weber-type gravitational 
wave detectors, have actually been special types of gradiometers . 

The acceleration difference along the direction j per unit separation 
along the direction i is defined to be the ij -component of the gravity 
gradient tensor T: 


r 


ij 


<r,t) 


5 2 4)(r,t) 


( 1 ) 


where <Kr,t) is the gravitational potential. A very weak gradient of 1 nm s 
per ra is equal to one Etttvbs (E) unit, defined by 


1 E 2 10 -9 s" 2 . 


( 2 ) 


Many ground-based survey applications call for such high sensitivity. Geodesy 
application in space requires even higher sensitivity at the level of 10 E 
jj z ~ 1/2 jjo] . The extreme weakness of gravitational interaction and the 
practical difficulties associated with balancing out the acceleration noise to 
a sufficient degree have limited the sensitivity of room temperature gradio- 
meters to a level of 1 ~ 10 E Hz“ 1/2 [10]. Major improvements in sensitivity 
and stability are expected of the superconducting devices under development. 

It appears that a superconducting gravity gradiometer of a relatively compact 
design will have a sufficient sensitivity for space applications. 
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Besides possessing low thermal noise and low mechanical drift as direct 


:cs of a cryogenic temperature, the new gtadloraeuer takes advantage 
of many exotic properties of superconductivity. Quantized magnetic flux is 
used as an extremely stable tool to achieve transducer action and common mode 
balance. Operating at liquid helium temperatures, SQUID (Superconducting 


Quantum Interference Device) serves as the most sensitive amplifier of today. 
Superconductivity can be used to make a nearly perfect electromagnetic shield 
and superfluid helium can provide a stable, gradient-free temperature environ- 
ment. Flux quantization can further be used to accomplish stable levitation 


of proof masses against gravity in a terrestrial environment and to enhance 


the gradiometer sensitivity by means of a "superconducting negative spring" 

[ 11 ]. 


Two schemes (current-differencing and displacement-differencing) of 
superconducting gravity gradiometer have been demonstrated by Paik et al [12]. 
Error models were analyzed by Wang [13]. Mapoles [14] has extended the de- 
velopment of a displacement-differencing gravity gradiometer. In this work, 
we have chosen the current-differencing scheme. One advantageous feature of 
the current-differencing gradiometer is the remote coupling of the two dif- 
ferencing acceleration transducers independent of their separation and their 
respective orientations. Therefore, three in-line (or diagonal) component 
gradiometers can be combined together by mounting all three pairs of accelera- 
tion transducers on the six faces of a common cube, with the sensitive axes 
normal to the surfaces of the cube, to form a three-axis in-line component 
gravity gradiometer. Construction of cross (or off-diagonal) component 
gravity gradiometer is feasible by orienting the sensitive axes of the 
acceleration transducers perpendicular to the direction of the baseline. A 
tensor gravity gradiometer to measure all the six components has been 
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proposed [15] as a combination of the in-line and cross component gradio- 
meters. 

While developing a three-axis gradiometer for precision gravity experi- 
ments in Earth's orbit [16], we have completed, with our colleague, a proto- 
type single-axis in-line component gradiometer. This instrument has been used 
to perform a laboratory null test of the gravitational inverse square law 
[17]. The details of this gradiometer development are described in Ref. 18 in 
which the theory of the gradiometer has been given a new formulation with 
generalization and more rigor than the preliminary analysis in Ref. 13. This 
paper (I) is a modified version of this new theoretical analysis with an 
extended error model of the instrument. Paper II presents the construction 
and test results of the gradiometer. Although we confine ourselves to the 
discussion of an in-line component gradiometer with a particular 
superconducting circuit chosen, the methods developed in these papers could 
easily be adapted to cross-component gradiometers and different 
superconducting circuits. 
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II. PRINCIPLE CF GRAVITY GRADIENT DETECTION 


In this section we briefly review the basic principle of separating 
gravity signal from dynamical variables and set a basis for the error model 
developed in Section V. Let the instrument platform be moving with respect to 
an inertial frame with an instantaneous angular velocity i!$(t) and linear 
acceleration a(t). Then, the accelerations of a proof mass observed in the 
two coordinate systems are related by the well-known equation [19] : 


, 2 + 

dt 2Jin 


, 2 ± 


(HU + 8 


dt 


x r) 




(dr\ fd!^ 

(jtV + ( dT'pt * r + 4 


.(3) 


Here the subscripts in and pJl represent the inertial and platform coordinate 
systems in which respective measurements are made. The second and third terms 
on the right-hand side are the centrifugal and the Coriolis accelerations, re- 
spectively. If <|>(r,t) is the gravitational potential in the inertial frame, 
the resulting acceleration of the proof mass with respect to the moving plat- 
form is given by 


g’(r,t) = (4-y) . 

dt px 

= “ $ 4>(r,t) - S$ x x r) - 2$ x (*[£•) x r - a(t) . (4) 

The Coriolis terra produces a force perpendicular to the velocity in the plat- 
form coordinates and therefore drops out when the proof mass is confined to 
move in a single direction. The quantity g'(r,t) is what is measured by an 
accelerometer or a gravimeter undergoing an acceleration. 

It is clear from Eq. (4) that the linear acceleration term, - a(t), can 
be eliminated by a differential measurement over a spatial coordinate Xj : 
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?>g!(r,t) 

9x. 


I »(r,t) . a [„ ( J . 

dx^ 9Xj 5Xj i 


Q 2 x i ] 


dQ. dx . 

r r k\ x 

. \ S ikA dt pA 9x. 

j 


( 5 ) 


Substituting 5x^/Sxj = 6jjj , and introducing the notations of Eq. (1) and 


ag’(r,t) 

r il< r > t > E — sp • 


( 6 ) 


dQ, 

“k^ = ( dt ’ 


( 7 ) 


one finds 


r^(r,t) = r ±j (r,t) - (Q 1 Q j -Q 2 ^) 


"ijk 


\ (t) 


( 8 ) 


Notice that the angular acceleration term is antisymmetric whereas the 
first two terms in Eq. (8) are symmetric tensors. Therefore, one can further 
drop a(t) by symmetrization of F!.: 


.. <? t> = i r*i 

(ij)'‘ r,t; " 2 l 9x. 


9g- 


9x 


= f i j (r,t) - (Q ± Qj - Q 6 tj ) 


( 9 ) 


The centrifugal acceleration term can be taken out in principle by taking 
another spatial derivative: i.e., by means of a third order gravity gradio- 

meter. In practice, one measures with the aid of gyroscopes and removes the 
effect of the centrifugal acceleration by actively stabilizing the platform or 
by compensating the induced error. 

The diagonal component of the tensor, , can be measured by 

detecting the relative acceleration along the in-line direction between two 
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proof masses, separated in the x^ direction. The off-diagonal components 
^(ij) ^ * *) ; hnw^vpr. requires four proof masses because of the symmetri- 
zation. The relative accelerations in the cross directions in two pairs of 
accelerometers, which are separated in the and Xj directions, respectively, 
can be added to yield Notice that one can instead subtract these two 

signals to determine the antisymmetric component : 



( 10 ) 


The angular acceleration of the platform can be obtained by inverting this 
equation : 


a k (t) - e i j k r'^jCr.t) . 


(ID 


A time integration of this vector then gives an alternative means of 
determining the angular velocity • 

A tensor gravity gradiometer with common mode readouts is therefore self- 
sufficient for true gravity detection [15]. On the other hand, the device can 
measure true linear and angular accelerations of a moving platform by removing 
gravity-induced errors. Equation (4) shows that the gravitational field -^4> 
remains as the fundamental error in linear acceleration measurement after 
removing dynamical error terms. The gravity gradiometer comes to the rescue. 
The gradient output V can be integrated over a spatial coordinate to determine 
Therefore, a true accelerometer requires an aid from a gradiometer. 

The symmetrization technique discussed above has been incorporated in the 
rotating gravity gradiometers [6, 7] whereas, in the floated gradiometer which 
has only two proof masses [8], the angular motion of the gradiometer is 
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attenuated by floating the proof masses in a liquid. In the rotating 
gradiometers, the common mode acceleration g' is further rejected by its 
frequency characteristic. The gradient ?, being a tensor of rank 2, is 
modulated at the second harmonic of the rotation frequency whereas the 
acceleration g' is modulated at the fundamental frequency by its vector nature 
[6], A side benefit of this heterodyne detection is the translation of signal 
bandwidth away from the 1/f noise region of the instrument in frequency space. 
The mechanical rotation, however, brings in a penalty: additional, 
dynamically induced errors. From Eq. (8) one can clearly see, for example, 
the devastating effect of the angular velocity error 6Q, which now contributes 
a first order term 0(Q6Q) to the measurement. 

For the prototype superconducting gravity gradiometer, we have chosen a 
non-rotating configuration. The extreme stability of the superconducting 
sensing circuit, combined with the low noise of the SQUID amplifier down to 
low signal frequencies, permits a very high degree of common mode rejection 
without rotation. For orbital applications, however, the superconducting 
gradiometer could be rotated to its advantage by spinning the entire satellite 
quietly. 

The symmetric nature of has been used to construct a gravity 
gradiometer. Further, the trace of this tensor is constrained by the Poisson 
equation: 


£ (r,t) = - V 2 <(>(r,t) = - 4-rtG p(r,t), (12) 

i 

which is a consequence of the inverse square law of the gravitational force. 
This leaves five independent components for the gravity gradient tensor r^j. 
With a three-axis diagonal component gravity gradiometer, the validity of Eq. 
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(12) could be tested by summing the three outputs and comparing the result 
with the local mass density p. This experiment has been proposed as a 
precision null test of the inverse square law [20] , and an early result of 
such an experiment has been reported [17], In an actual experiment, the 
measured quantity is the trace of 

l r’ (ii) (r,t) - - 4itG p(r,t) + 2 Q 2 (t). (13) 

i 

It is therefore important to suppress or separate out the centrifugal 
acceleration terra carefully in such an experiment. 
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III. DYNAMICS OF THE SUPERCONDUCTING GRAVITY GRADICMETER 

The superconducting gravity gradiometer consists of a pair of supercon- 
ducting acceleration transducers and a superconducting inductive load which is 
connected to a SQUID amplifier. The coupling between the transducers and the 
output load is provided by flux quantization. 

The principle of one acceleration transducer element is first discussed 
and its equation of motion is then derived. Each transducer communicates to 
the rest of the superconducting circuit only through one current component 
that flows through the transducer. In a coupled circuit of a pair of trans- 
ducers and a load, flux quantization imposes constraints to the supercon- 
ducting circuit. For detection purposes, the currents at the load , rather 
than the current through the transducers, are the observable quantities at the 
output. The complete dynamical equations are then linearized and expressed in 
terms of these currents, the respective displacements of the proof masses, to- 
gether with the applied common and differential acceleration signals [21], 
Because the gradiometer is a differential accelerometer over a finite base- 
line, the gradiometer must reject common acceleration signals. The conditions 
for a common mode balance is derived. The parameters used to accomplish the 
balance are the persistent, but adjustable, currents stored in various super- 
conducting loops. Further, a wideband common mode balance is shown to be 
possible by iteratively adjusting at least two current components. After 
balancing the common accelerations, only differential acceleration will be 
detected at the load. The transfer function of an applied differential 
acceleration to the corresponding current output at the load is derived. 

Throughout the remainder of this paper, a variable with time or frequency 
dependency will be written explicitly as such functions, whereas the average 
values of these variables will be denoted by the same notation as the 
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function, but with the functional dependency deleted. 

A. Principle of a Superconducting Acceleration Transducer 

The principle of the superconducting acceleration transducer is illus- 
trated schematically in Fig. 1. Analysis of this device as a resonant trans- 
ducer for a gravitational wave antenna and as a sensitive accelerometer has 
been given previously in a different format [9], A superconducting proof 
mass, which is suspended by spring and is confined to one linear degree of 
freedom, responds to an acceleration signal with a displacement relative to a 
sensing coil. The inductance of the coil is then modulated, due to the 
Meissner effect, by the superconducting plane of the proof mass. The coil is 
connected to an output inductor through a superconducting path and a quantized 
magnetic flux is stored in the superconducting loop formed by the sensing and 
output inductors. The current flow through the output inductor is modulated 
as a result of the inductance modulation of the sensing coil. The persistent 
current provides the stability of the transducer scale factor. A low noise 
SQUID amplifier is then used as a dc current-to-voltage power amplifier to 
produce a readout. 

B. Analysis of a Single Acceleration Transducer 

The transfer of mechanical energy to electrical energy is accomplished at 
the sensing coil. In order that the sensing coil converts a displacement to a 
current more linearly within a transducer, a symmetrical pair of "pancake" 
coils are utilized. The coils are located on the opposite faces of the proof 
mass and are connected in parallel (Fig. 2). Each having winding density n^ 
and area A^, the coils are at mean distances of d a and d^ from the respective 
superconducting planes of the proof mass. If the displacement of the proof 


- 15 - 



- 16 - 


Fig. 1. Principle of the superconducting acceleration transducer. 




Fig. 2. Superconducting circuit and its current variables of a single 
acceleration transducer . 
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mass from its average position is x(t), the inductances of the coils are given 
by 


L a (t) = A [d a + x(t)] = L a + A x(t), (14a) 

Ljj(t) = A [d^ - x(t)] = Ljj - A x(t). (14b) 


Here 


A = n^ A l , (15) 

where |i Q is the permeability in free space and 

L a " <L a (t)> » H 5 <V C)> » (16) 

as was noted earlier concerning notation. 

Two current components are needed to characterize the electromagnetic 
state of the two sensing coils. One obvious choice of variables [9, 13] is 
the currents X a (t) and I^(t) which flow through each of the two coils, so that 
the electromagnetic energy in the two coils is 

V EM * k V e > V ') 2 + } V*> < 17 > 

However, in order to make the analysis simpler, a different choice [18] of the 
two current variables can be made with the aim that the expression for the 
electromagnetic energy has as simple a denominator as possible. Such a way of 
choosing variables is a classical analog of the renormalization procedure in 
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quantum field theory. 


Before making the choice, we make two observations. The first 
observation is that the series inductance of the sensing coils is a constant 
19]: 


L s (t > 5 + " L a + h " V 08) 

Consequently, the current I through the superconducting loop of these two 
inductors in series is also a constant because the trapped flux $ ab in this 
loop is quantized. The second observation is that the parallel combination of 
the two sensing coils has a constant denominator in the expression for its 
inductance : 

L a (t) L b (t) 

L P (t) " L a (t) + L^t) “ L a + 1^ lL aM) " (L a " V A x ( c > " A 2 x 2 (t)]. 

(19) 

In fact, even the numerator of this parallel inductance will also be a 
constant up to the first order if the mean spacings d a and d b are matched so 
as to make L a * Ljj. The nonlinearity of the inductance modulation is 
exhibited by the second order term in Eq. (19). 

Expressed in terms of one parameter $ ab and one current variable i(t), 
which flows through the parallel combination of L a (t) and 1^(0, the 
electromagnetic energy has a constant denominator and has thus acquired a 
"renormalized” form: 


$ 


ab 


'EM 


2 L a + *t> 


+ iV° 


i 2 (t), 


( 20 ) 
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It is straightforward to show that the expressions in Eqs. (17) and (20) are 
equivalent to each other [18]. 

The only variable i(t), once the parameter I or is fixed, possesses 
all the dynamical information in the electromagnetic system of the transducer. 
As the inductances L a (t) and Ljj(t) are modulated by the displacement of the 
proof mass, the current i(t) must always split between the two inductors 
according to the inverse ratio of the respective inductances in order to pro- 
duce equal and opposite magnetic flux contributions to the superconducting 
loop formed by the series inductor L a + L^. The net currents I a (t) and I^(t) 
through L a (t) and L^t) are therefore the following linear combinations of I 
and i(t): 


L. (t) 

vo-nr 11 ' 1 - 1 ' <2U) 

a b 
L (t) 

Ib<t> " *'a**~H 1(t> + 1 <21b) 

The force due to magnetic field pressure on the proof mass is given by 

$*(t) $j;(t) 

v - _ L. r a + b l | 

*EM - ax 1-2 L a (t) 2 L^t)- 1 '® a ,\ 

= | [I*(t) - l£(t)] , (22) 


where $ a (t) and $^(t) denote magnetic fluxes in L a (t) and L^Ct), respectively. 
In terms of I and i(t), this force can be rewritten as 


f em 


- A [ I + \ 


L a (t) - L b (t) 
L a + H> 


i(t)l i(t). 


(23) 
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The total force acting on the proof mass consists of this magnetic force 
together th a restoring force of the mechanical spring suspension of the 
proof mass and any externally applied force, f(t), on the proof mass relative 
to the platform of the sensing coils. We will ignore the effects of damping 
term for the time being. When the proof mass m is approximated as a point 
mass located at its center of mass r, the external force becomes 

f(t) = m n • g'(r,t) , (24) 


where n is the unit vector along the direction of the sensitive axis, and 
g'(r,t) is the specific force given by Eq. (4), The equation of motion for 
the proof mass in the platform frame can now be written as 


x'(t) + U> 2 x * ( t ) + ^ [i 


1 L a (t > - L b (t) 

2 L a + L b 


i(t) ] i(t) 


n • g' (r,t) ,(25) 


where u> M is the angular resonance frequency of the mechanical spring. An 
alternative derivation of this equation, in which an electromechanical 
Lagrangian approach is used for the superconducting transducer, is found in 
Ref. 18. 

The dc component of Eq. (25) defines the equilibrium position, x Q , of the 
proof mass. This position is shifted from the relaxed position of the 
mechanical spring by Earth's gravity, g E , and by a dc magnetic force. Substi- 
tuting 


x ' ( t ) = x Q + x(t). 


(26a) 


g'(r,t) = g E (r) + gp(r,t) 


(26b) 
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into Eq. (25), we find 
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*E cos 9 n A ,, 
2 2 '• A 
"“m 


,L • 

+ i_l 
21 
a 



0 i , 


(27) 


In the general case where i = <i(t)> * 0. Here 0 n Is the angle that the 
sensitive axis makes with the upward vertical. The local vertical is defined 
by Earth's gravity vector which is found from Eq. (4) as 


gg(r) *» - $ 4>g(r) ~ x t^ E x (tg + r)] , (28) 

where <j>p(r) , (5^ and ftg are, respectively, Earth's gravitational potential, 
spin angular velocity, and the geocentric position vector of the coordinate 
origin. 

The linearized equation of motion for the ac part, after Fourier trans- 
formation, can be readily shown as 

2 2 

(-w 2 + ) x(w) + i(o>) - g(w) , (29) 

s 

where g(u>) is the Fourier transform of 

g(t) H n • gp(r ,t) (30) 


and 


I' 


I + 


L - 
a 


+ h 


i. 


(31) 


There are two independent superconducting loops in the circuit. The flux 
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quantization condition in the loop L a (t) + 1^(0 has already been used in 
icptcociiLiug j. as a constant* An additional constraint in the dynamical 
equation comes from the flux quantization condition in the loop 1^(0 + L Q 


L|>(t) + “ $ bo‘ 


(32) 


The first order equation becomes, after Fourier transformation, 


Al'x(to) = (L q + Lp) i(w). (33) 

Simultaneous equations (29) and (33) determine the dynamics of a single 
transducer completely. In particular, the acceleration-to-current transfer 
function can be shown to be 


H gi (o>) = 


to 


AI* 


2 L + 

(i) o 


(34) 


Here, o) Q is the resulting angular resonance frequency of the transducer due to 

2 

the addition of spring constant to mw^ by the superconducting circuit: 



0 , 2.2 „ 2 t ,2 

2 A x A I' 

M mL m(L + L ) 
s op 


(35) 


C. Dynamical Equations for the Gradioneter 

The gradiometer, shown in Fig. 3, consists of a pair of the above 
acceleration tranducers and a superconducting inductive load. The coupling 
among these elements is through flux quantization in the superconducting 
circuit which has four independent superconducting loops. 


- 23 - 






The dynamical equations for each transducer is analogous to Eqs. (29) and 
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A Artnl nnVio /tv 1 ! nt* 

fill UUULLAUUUX yUU«J VJ. J- J/ w < 


k * 1 or 2, (36) 

is now used for the respective variables and parameters to distinguish between 
the two transducers. With this subscript, the linearized equations of motion 
for the proof masses in the transducers, being similar to Eq. (29), are 
written as 

, 2 2 y k K 

+ “kM + V w) + = «k (a,) * (3?) 

Here, the geometries of the sensing coils are again assumed to be identical so 
that they can be represented by a single parameter A. defined by Eq. (15). 

Also, the two flux quantization constraints in the superconducting loops + 

have again been used in deriving Eq. (37). Two more constraints, similar 
to Eq. (33), are obtained by using flux quantization conditions in two other 
independent superconducting loops such as 1^(0 + L Q . However, a generaliza- 
tion from Eq. (33) is needed here because the current through L Q is now the 
sum of the current outputs from each transducer. The linearized constraints 
are therefore given by 

A x k (w) = i k (w) + L q [i ^(o>) + i 2 (o>) ] . (38) 

Notice that these two equations are coupled through ij(w) and ^(w). 

The dynamics of the gradiometer, governed by the four coupled equations 
(37) and (38) in the four variables x^(w) and i^(w), has two degrees of 


- 25 - 


freedom. Each of the two transducers couple to the rest of the circuit 
through one single current signal, ij(io) or ^(w). Depending on the signature 
of the applied acceleration components, these current modulations can add or 
subtract at the output inductor L Q . It is therefore convenient to use a new 
set of current variables i d (to) and i c (w) defined by the sum and half of the 
difference of the two transducer currents. In terms of i d («) and i c (w), the 
transducer currents ij(co) and ^(w) can be expressed as 

i k (w) - \ i d (w) + (-l) k i c (u>). (39) 

Likewise, the applied accelerations gj’(t) and g 2 (t) at the proof masses m^ and 
m 2 can be expressed in terms of their differential and common accelerations , 
g d (w) and g c (w), as 


g k (u>) = (-l) k j g d (u>) + g c (w). 


(40) 


Notice that i d (<n) is the actual current flowing through L Q and detected by the 
SQUID. 

With the change of variables in Eqs. (39) and (40), the dynamical 
equations (37) and (38) become 


A 2 ! 2 


• i. j. A. L 1 

(-u, 2 + *4 + snr-> Xk(w) + ~sr 1 1 d (w) + (_1) 


K K 


W, 


”k 


k -ST 

“k c 


(-l) k J g d (u>) + g c (w), 


(41) 


x k (ui) = \v j i d (w) + (-l) k i c (w). 


(42) 
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Upon eliminating the displacement coordinates, we obtain 




+ 2L 


'let 


(-to 2 + a£ d ) — J + (_1)k (_a)2 + V ) IT 1 c (w) 


(~D k J g d («»>) + g c («) = gk^ w ^» 


( 43 ) 


where 

(44a) 
(44a) 

Solving Eq. (43) for the signal current, we find 

i d (w) ** Hg i (o>) g d (co) + Hg i (a>) g c (co), (45) 


» 2,2 a 2 t ,2 
2 = 2 J_ r A St . A *k ^ 

^kM m, k L, _ ^ * 
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. 2,2 , 2 t ,2 


2 _ 2 i r A *k A T k ^ 

^kd ^kM m, ^ L, L, + 2L 
K ks tep o 


Q C 

where H^(w) and Hg^(u) are the transfer functions from g d (u>) and g c (u), 
respectively, to i d (w). These transfer functions of the gradiometer are given 
by 
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The four parameters and o>j cc are, in general, not the two resonances 
of the gradiometer. Nevertheless, comparing the expressions for these para- 
meters with Eqs. (34) and (35) suggests intuitive meaning to these nonobserv- 
ables. Thus, (0^,, can represent the angular resonance frequencies that the two 
transducers would separately have if they were each connected to a short- 
circuited load. A short-circuited load arises in these models because, when 
id(to) = 0 in the gradiometer circuit, the two transducers are driving the 
current i c (w) in a push-pull manner, contributing zero impedance to each 
other. Likewise, can represent the angular resonance frequencies that the 

two transducers would separately have if they were each loaded with an 
inductance of 2L Q but were otherwise decoupled from each other. The 
appearance of 2L Q here is due to the equal contributions of currents y i < j( w ) 
to the signal at L Q . Thus, under the restriction i c (w) = 0, the flux 
modulation produced by each transducer at L Q as seen by the transducer itself 
is doubled by the presence of the other transducer. The transfer functions of 
these separated model accelerometers are given from Eq. (43) as 


<-l>X 1_ 

bg. ( i 2L + L 2 2 ’ 

k C O Kp - U) 

M (-i)Vi^ ! 

Vhi hep «4 C - » 2 ‘ 


(47a) 


(47b) 


D. Common Mode Balance and Gradiometer Transfer Function 

For operation as a gradiometer in the presence of common accelerations, 

c 

the transducers need to be tuned such that H^(to) = 0. No common acceleration 
signal will then appear at the output load L Q , and the coupled acceleration 
transducers are said to operate in a "gradiometer mode". On the other hand, 
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parameters can be chosen such that H^(io) vanishes instead. Then, the output 
will respond only to common accelerations and the system is said to operate in 
an "accelerometer mode". The accelerometer mode operation is a powerful means 
of obtaining a precise calibration of the gradiometer. A detailed discussion 
of the accelerometer mode is found in Section IVB of Paper II. From Eq. 

(46b), the common mode balance is obtained when 


( “i c - “ 2) L i P 


m ' 2 

71 2 TT 

W 2c " W > L 2p 


(48) 


This balance condition can be satisfied at any single frequency by adjusting 
only one persistent current parameter. Balancing over a small frequency range 
near dc appears sufficient for normal low-frequency use of the gradiometer in 
which a) « Wj c , o> 2 c . In the terrestlal environment, however, the 
environmental vibrations occur in a wideband and are very large compared with 
the extremely weak gravitational signals. A wideband balance will help to 
immunize the gradiometer against such environmental vibrations. 

The four current components, Ij, I 2 , ij , and i 2 , do provide more than 
sufficient degrees of freedom to tune for wideband balance. In principle, Eq. 
(48) becomes an identity if 


“lc = “2c 


(49a) 


and 


A If UJ 


(49b) 


are simultaneously satisfied. However, direct matching of o)| tc is not 


- 29 - 


practical because these equivalent resonances are observable only when the 
output load is shorted and when the two transducers are also decoupled 
mechanically [18]. In the experiment, a wideband balance can be achieved 
instead by iterating the balance at two frequencies. Tuning 1^ or I 2 to 
balance out an applied common acceleration at a high frequency (co » Wj c , w 2c^ 
will achieve the condition (49b). This adjustment is followed by tuning of ij^ 
or i 2 for balance at a low frequency (w « Wj c , ( 02 c )» yielding 


1 * X i 


1 * 4 


(49c) 


<0 


lc 


J lp 
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2c 


2p 


The latter operation will, in general, affect the previous balance (49b), and 
iteration between the two balance procedures Is needed. When the conditions 
(49b) and (49c) are both satisfied, the condition (49a) follows. This fre- 
quency-independent balance has been applied in the experiment and is reported 
in Paper II. 


Upon substitution of Eq. (48) into Eq. (46a), H^(to) assumes a simple 


form: 
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r „d , ,,-1 ,2 2 n o 2 lp , t 2 
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.(50) 


Notice that the transfer function of a balanced gradiometer becomes the 
harmonic mean of the transfer functions of the two separated model 
accelerometers, given by Eq. (47a). 

When the common mode balance is not precise, H g ^(co) * 0 constitutes an 
error coefficient. This error, which will be discussed in Section V, is 
obtained from Eqs. (46) as 
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In the foregoing analysis , the transducer has been approximated as a linear 
system. The dynamical equations (25) and (32), however, are nonlinear in 
nature so that the response current i^ta) must contain nonlinear terms in 
general, in addition to the linear terms shown in Eq. (45). These nonlinear 
terras may not drop out even when H ^(m) * 0. The scale factor nonlinearity 
therefore constitutes an important error source, which will be discussed in 
Section V of this paper and in Appendix B of Paper II. 


E. Equivalent Accelerometer Representation 

With two coupled acceleration transducers in the gradiometer, only the 
differential acceleration is the measured signal at L Q . The common accelera- 
tion drops out upon balance. Therefore, a convenient representation of the 
gradiometer is a single (differential) accelerometer which converts a dif- 
ferential acceleration into a signal current. Such a representation involves 
identification of the gradiometer transfer function with the parameters of a 
single accelerometer , as discussed in Section IIIB. These parameters include 
a resonance frequency, a mass and the parameters of a superconducting circuit. 

The two normal mode resonance frequencies of the gradiometer are some 
weighted average of the parameters Wj c and u)2 c f° r common mode and also 
some other weighted average of and for the differential mode. In the 
particular case of wideband balance, ( 0 j c and o» 2 c are equal to each other and 
hence equal to the common mode angular resonance frequency, denoted by w c . 

The differential mode angular resonance frequency, denoted by u) Q , is then 
defined from the singularity of H^(co) in Eq. (50). In practice, the gradio- 
meter need only be balanced at a limited frequency band of the signal. Then, 
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Eq. (50) is only valid in this limited frequency band. Therefore, the exact 
differential mode resonance is determined from the singularity of H^(cj) la 
the more general equation (46a). Yet, as far as signal transduction within 
the limited signal frequency band is concerned, we can still use Eq. (50) to 
define u) Q . 

By defining equivalent circuit parameters to satisfy 
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(52) 


we can convert Eq. (50) into the form of Eq. (34): 
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(53) 


where the effective resonance frequency, co 0 , is related to by 


1 


l o + L p 2 . L o *2 L lp 2 . 
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(54) 


While there is freedom in defining the individual circuit parameters 
which appear on the lefthand side of Eq. (52), one natural choice is to 
identify L Q with the load inductance and Lp with the output inductance of the 
superconducting circuit : 


L 

P 



(55) 


Then, the current parameter I' is defined by Eq. (52). With the aid of Eqs. 
(49b) and (55), this definition of I' in Eq. (52) simplifies to 
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( 56 ) 
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under wideband balance condition. 

The transfer function In Eq. (53), has only characterized the overall 
signal transduction of the gradioraeter. The intrinsic noise of the gradio- 
meter will be given in Section IV in terms of Brownian motion noise and 
amplifier noise. For a single spring-mass system, the Brownian motion depends 
on the mass and the fluctuating force of the spring. Therefore, once the 
effective mass in the equivalent accelerometer representation is determined, 
the Brownian motion of the two coupled acceleration transducers in the dif- 
ferential mode can be derived using this equivalent single spring-mass model. 
This effective mass requires a unique definition in order to give the correct 
Brownian motion noise of the gradlometer and must therefore be derived from 
the dynamics of the gradlometer. 

Dynamically, the gradlometer consists of two coupled masses connected to 
the platform and to each other by three springs. Thus, the homogeneous 
equations of motion of the two proof masses are obtained by eliminating i k ((o) 
from Eqs. (37) and (38): 


uTm k x k (w) = K k x k (a)) + K 12 X 3 _ k (<*>), 


(57) 


where the three spring constants 1^ and Kj 2 are given by 
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Now, we need to find out what displacement variable, x d (to), cause the signal 
current i d (w). Thus, we break up the signal transduction into two 
intermediate steps: g d (w) to x d ((i>) and then x d (c»)) to i d (<o). From Eqs. (41) 

and (42), we obtain 

( —to 2 + i/) x d (oj) = g d (w), (60) 

i d (<d) = L ^ L " x d (a>), (61) 

o p 

where 

L I T ~I f 

x d (w) = -j?- x 2 (w) - (<»>)]• ( 62 ) 

2 p lp 

The effective mass is therefore the inertia for the coupled spring mass system 
of Eqs. (57) - (59) towards the motion x d (w) of Eq. (62). 

We next note that a simple scaling of the two separate displacement 
variables, 

L (-l) k I* 

x£(to) = — jr x k (to), (63) 

\p 

simplifies x d (to) into the form: 


x d (u>) = x^u) - xj(co). 


(64) 


The corresponding scaling required for the masses and spring constants are 
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(65b) 
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because, under these scaling, the dynamical equation (60) is invariant in form 
and both the kinetic and potential energies (within the frequency range of 
common mode balance) remain unchanged. 

Now, the effective mass corresponding to x d (o>) of Eq. (64) is just the 
reduced mass of mj and m£. By using Eqs. (63) and (65a), this effective mass, 
denoted by m, can be shown to satisfy a simple formula: 
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( 66 ) 


In Section IVA, the effective mass defined here will be used in 
conjunction with the equivalent accelerometer representation for finding the 
Brownian motion noise of the transducer. 


- 35 - 



IV. NOISE AND OPTIMIZATION 

In the light of the equivalent accelerometer model given in Section IIIE, 
the fundamental noise terms of the gradiometer will be derived by treating the 
gradiometer as a single accelerometer. In this section, we introduce damping 
and consider its effects. 

A. Transducer Brownian Motion Noise 

The fundamental noise source of a spring-mass system at a temperature T 
is the Brownian motion of the harmonic oscillator at that temperature. The 
force of the spring on the proof mass undergoes random fluctuations. The 
spectral density of such force fluctuations at the resonance frequency u 0 of 
the oscillator is related, by the Fluctuation-Dissipation Theorem, to the 
damping of the oscillator at resonance. 

In applying the Fluctuation-Dissipation Theorem, emphasis is made here 
that measurement of relaxation time x(u) 0 ) or quality factor Q(n> 0 ) 5 u> 0 t(w 0 ) 
gives information on the force fluctuations only at u) Q . The Langevin equation 
is a modification of the equation of motion in Section III by including a 
damping term and an acceleration noise g n (w) term: 

[-“ 2 + t(uT + X(w) = g n (a>) * (67) 

o 

Thus the energy of the oscillator at a temperature T: 

m w 2 <x 2 (t)> = k fi T, (68) 


where 
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<x 2 (t)> = x(oa) x*(u) dw. 


(69) 


is dominated by its spectral density at resonance. Use of Eqs. (67) - (69) 
gives the Nyquist formula: 


) = - k T — r- , 

go it B mt(a) ) 


(70) 


T 

where S (w ) is the spectral density [22] of the acceleration noise g T at uu. 
6 ° °n o 

The force or acceleration fluctuations at a signal frequency w « co Q is 
in general different from that at w Q because the noise could have a compli- 
cated frequency dependence which is governed by the nature of the loss 
mechanisms in the spring. The noise at w due to force fluctuations in the 
spring is obtained from Eq. (70) with t(u) q ) replaced by t(w) : 


S T (o)) = - k T — 7-v- 
g x B mT(o)) 


(71) 


A direct measurement of t(u)) requires shifting the resonance frequency from u) Q 
to a). In principle, one can increase the mass m or use a "negative spring" 
[11] so that the new spring-mass system indeed resonates at to and therefore 
the Q factor of this new system can be measured. From a knowledge of the new 
mass used or of the loss in the negative spring, the dissipation in the 
original spring can be determined. 

The Q factor of the accelerometer depends on both the mechanical and the 
electrical parameters because the electromechanical spring of the accelero- 
meter has contributions from a mechanical spring and two electromagnetic 
springs, in the ratio given by the three terms in Eq. (35), for the total 
spring constant and hence for the stored ac energy at w 0 . The power loss in 
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each of these constituent springs is then proportional to this ratio of stored 
ac energies in these springs divided by a certain relaxation time at w Q . 
Therefore, the following relationship for the relaxation times results from 
the obvious observation that, at <0 o , the total ac power loss is the sum of the 
power losses in the respective springs: 


•c(u>o> 
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(72) 


Here t m (< 0 o ), 'c i (a) Q ) and Tj;.((0 o ) are the relaxation times that each spring 
would have if it is separately resonated at <*> 0 , and y and 0 are the transducer 
energy coupling coefficients via i and I 1 defined by 
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In order to reduce dissipation in all the constituent springs, choice of 
material and geometry, material treatment, surface preparation as well as 
choice of electrical parameters must be optimized. 


B. SQUID Amplifier Noise 

A SQUID can be modelled as an ideal current-to-voltage amplifier having 
an input inductor and two noise generators for its voltage and current 
noise with spectral densities Sy(u)) and S I (w), respectively, at the amplifier 
input (see Fig. 4a) [23]. 

The optimum source impedance, (Sy/Sj) 1 ^ 2 , is much smaller than the para- 
sitic impedance unless L t is tuned out (at one frequency) to noise match 
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Fig. 4(b). An equivalent circuit of the SQUID noise for an untuned SQUID 



the input impedance. Therefore, when the SQUID is used as a wideband ampli- 
fier, the noise is dominated by S x and can be characterized by an "input 
energy sensitivity" : 


Ea(w) = "2 Sj(to). 


(74) 


In a practical gradiometer design, the transducer output inductance Lp is 
usually larger than and a superconducting transformer is desired to bridge 
between 1^ and l^. The equivalent load, representing Lj and the transformer, 
seen by the transducer is an inductor L Q : 
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(75) 


where k t , Lj and L It are the magnetic-field coupling constant, the primary and 
the secondary inductances of the transformer. The ratio of the power being 
detected in Lj^ to the power sent from I^ to L 0 is the forward power transfer 
function H I,II (w) of the transformer and can be shown to be given by a fre- 

p 

quency-independent expression: 
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The equivalent current noise for Sj(to) seen by the transducer is therefore 
given by 
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Now, the circuit in Fig. 4a can be represented by an equivalent circuit with 
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the transducer connected In series with L 0 and a current-to-voltage amplifier 
which haa the equivalent current noise Sj(iu), as shown in Fig. 4b. 

The quantity of interest is the equivalent (differential) acceleration 
noise at the transducer that would correspond to the amplifier noise of Eq. 
(77). The conversion from current noise to differential acceleration noise is 
made by use of the transfer function of the transducer, Eq. (53): 


S 8 (W) 


((/ - to )(L + L ) 0 2 E.(w) 

r O O p -|2 A 


ti i' 


L H 1 ’ 11 
o p 


(78) 


For (o << a) , substitution of Eq. (73b) reduces this amplifier noise term to 


S ^5^- E (w) , 
g m 2 pq A 


(79) 


where 




H 1 ’ 11 

P 


(80) 


is the fraction of electrical energy coupled to the amplifier. The product 0n 
represents the fraction of the total electromechanical energy coupled to the 
SQUID input, or the "SQUID energy coupling coefficient" . 

There are three factors to optimize. The transfer function (75) is 
optimized when the transformer secondary is 


yielding 


( L II>opt 


/(I - kp 


(81) 
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( 82 ) 




opt 


[1 + /(I - k 2 )] 2 


From Eqs. (75) and (81), the transformer primary is chosen according to 


^I^opt 


/(I - kp 


(83) 


The impedance factor L Q /(L 0 + Lp) in Eq. (80) is maximized by choosing L Q » 
Lp. Finally, the frequency factor assumes a minimum: 

2 


(0 


f— ) 

^P'min 


= 0 ) 


(84) 


when I' is chosen such that P = 1/2: namely, 


A 2 !' 2 = J_ 2 

m(L Q + Lp) 2 W o 


(85) 


This last optimization is applicable only when the amplifier noise is 
dominating. 


C. Potential Sensitivity 

The fundamental noise of the gradiometer is the sum of its Brownian 
motion noise and SQUID amplifier noise. In terms of an equivalent gravity 
gradient noise T n , the spectral density of the gradiometer noise is 


Sr(cj) = *4 [sj(«) + S%)] , 


g 


( 86 ) 


where X is the length of the baseline between the two proof masses. 
Substitution from Eqs. (71) and (79) gives the one-sided noise spectral 
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density: 


S r (f) = 


m 


k R T < 

T (f) + 2^ E A (f) ^ 


(87) 


where a factor 2 % has been dropped from Eq. (71) by going from the angular 
frequency domain to the frequency domain. Here, the numeral 4 will become 8 
if m^ is used instead of the reduced mass m. 

In the usual case when > Tj?, the optimum value for I* satisfies an 
inequality : 


A 2 1' 2 

n ✓ opt , 1 2 

0 < m(L + L ) < 2 “o 


( 88 ) 


Toward the lower bound, the power coupled to the SQUID tends to zero. Toward 
the upper bound, the Brownian motion noise term increases because of the rela- 
tively short Tji. If the amplifier noise dominates, the choice is made at the 
upper bound. If this bound is exceeded, the electromechanical spring has be- 
come so stiff that a given acceleration signal applied to the proof mass is 
producing less displacement, thus making it more difficult for the supercon- 
ducting circuit to measure. 

The SQUID input energy sensitivity E^(f) is usually a white noise plus a 
1/f noise at very low frequencies. In addition, a dc drift can be caused by a 
temperature drift of the gradiometer [14]. The drift changes the penetration 
depth and hence the inductances of the superconducting circuit [24]. Tempera- 
ture related drifts can be suppressed by regulating the temperature or be 
balanced out with an improved circuit [25]. This and other errors are the 
subject that will be considered in Section V. The noise given in Eq. (87) 
represents the fundamental noise of the gradiometer which can be reached only 
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when other instrument errors are suppressed sufficiently, 
very low fundamental noise in the gradiometer itself is a 
achieving the desired high sensitivity. 


However, having 
prerequisite for 
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V. ERROR MODEL OF THE INSTRUMENT 


As we have seen in Section II, the gradioraeter measures in general a 
differential gravity signal superposed with specific pseudoforces such as 
centrifugal acceleration and angular acceleration. In addition, the large dc 
gravity bias is modulated by the motion of the platform, resulting in error 
signals which compete with the ac signals under investigation. Thermal and 
electromagnetic fluctuations of the environment can also be coupled to the 
gradioraeter. It is important to have a complete error model of the instrument 
because, for many applications, errors could dominate over the gravity 
gradient signals and therefore must be compensated for to recover the true 
signals. We start this section with a derivation of general expressions for 
the driving accelerations g d (t) and g c (t). We will then elaborate on specific 
error sources. The kinematic and dynamic error mechanisms discussed in this 
section are independent of the particular electrical circuit chosen for the 
superconducting gradiometer. Hence our results have general applicability 
beyond the specific instrument discussed in this paper. 

A. Geometrical Metrology Errors 

Let us consider the case in which the gradiometer is used to measure a 
time-varying in-line component gravity gradient at r in the direction of a 
unit vector n fixed in the laboratory frame: 

r nn (*,t) = ~ (n • $) 2 <fr(r,t) . (89) 

Let r^ and rj represent the actual position vectors for the centers of mass of 

A A 

the two proof masses m^ and m 2 , and n^ and ^ be the unit vectors representing 
the actual sensitive axes of the two component accelerometers, as indicated in 
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Fig. 5. Let p and X represent the average proof mass position vector and the 
baseline vector of the gradiometer, respectively, so that 


(90a) 
(90b) 

Ideally, one wants p * r. In general, one has a gradiometer position 
error 6p : 



* 

P 


-> 

r + 


6p 


(91) 


Likewise, in an ideal gradiometer, n^, ^ and A = X/l are perfectly 

A 

aligned to each other and are oriented along n. In a practical gradiometer, 
however, there are various alignment and orientation errors. The axes 
alignment errors consist of the misalignment between the sensitive axes of the 
two accelerometers : 


6n_ = n 2 “ n i > (92) 

and the misalignment between the average direction of the sensitive axes and 
the direction of the baseline X = X/X: 

6n A = j (n_ + n. ) - A . (93) 

+SL 

These alignment errors will be shown to cause coupling to the gravity gradient 

/V 

output from the common linear acceleration component along the 5n_ direction 
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and from the angular acceleration component along the 6n A x n direction. The 

a+A 

gradiometer orientation error is due to the directions A and (n 2 + n^/2 being 

A A A 

misoriented from n by 6A and 6 n A , respectively: 

+n 


A A 


A = n + 6A , 


(94a) 


j (n 2 + n^ ) = n + 6n A , 
4 +n 


(94b) 


and will be shown to cause coupling from the differential acceleration 


A A 


component along the 6A + 6n * direction. Notice that only one of the two 

A A + n A 

quantities 6JI and fin * in conjunction with the misalignment fin A in Eq. (93) 

+n +A 

is sufficient to define the misorientation of the gradiometer. Introducing 

additional notations in Eqs. (94), however, have the advantage of simplifying 

expressions for misorientation errors by avoiding explicit reference to the 

misalignment parameters, defined within the gradiometer hardware. 

In order to find the specific forces acting on mj and Eqs. (90) are 

substituted into Eq. (4). After Taylor series expansion in 1, one finds 


g'(r lf t) = - $11 - (j • h + jf 4 * ^ )2 - 3l ( f ’ v > 3 + •••] ♦(P.t) 


- &(t) x [&(t) x (p - y)] - a(t) x (p - !) - a(t) , (95a) 


g’(? 2 ,t) = - $ [1 + (| • h + (| • $) 2 + ^ (| • $) 3 + ...] «(.(p,t) 

- £(t) x [£(t) x (p + |)] - £(t) x (p + |) - a(t) , (95b) 


where we dropped the Coriolis terras assuming a rectilinear compliance of the 
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proof mass suspension. In general, p and t are also time-varying due to the 
translational and rotational motions of the platform with respect to an 
inertial frame. 

In Section IIIB, we defined as the acceleration signal the time-varying 
part over the constant background gravity of Earth. Since we are now 
interested in the mechanisms which modulate the Earth ’s bias to produce error 
signals, we consider in this section the response of the gradioraeter to the 
total accelerations: 


g£(t) = \ • g'(r k »t) = g k (t) + n • gg(r) . (96) 

We return at the end to the signal variables g^(t) to connect back to the 
analyses presented in earlier sections. With Earth’s gravity vector g„(r) as 
defined by Eq. (28), the last equality in Eq. (96) can be viewed as a rigorous 
re-definition of g k (t). Here a point-mass (or a spherical-mass) approximation 
of the proof masses has been used to avoid volume integration of g'(r^,t) over 
the finite dimensions of the proof masses. 

The total common and differential accelerations for the gradiometer , 
defined in a similar fashion as g c (t) and gj(t) in Eq. (40), can now be shown 
to be 


g’U) 


(n + 6n *) 


+n 


S^CS.O + \ 6n_ 


!• <P,t) , 

dI 


(97a) 


g'(t) = 6n_ • g' (p,t) + (n + 5n A ) • g' (p,t) , (97b) 

C$ +n dI 

where g' (p,t) and g' (p,t) are the total true common and differential 
Cl Dt 

accelerations acting at p ± 1 / 2 : 
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g^(p,t) = - $[1 + jr (| • ^) 2 + ...] <t>( p) - l(p) , (98a) 


g’ (p,t) = - hd • h + fr 4 

dI 


V) 3 + • • . ] 4> (p) - ^ x (S x i) 


a x $.(98b) 




Here a(p) is the total linear acceleration experienced at 


+ 

p: 


a(p) = !$x(<$xp)+axp + a 


(99) 


and the time dependencies have been omitted on the righthand side of Eqs. (98) 
for notational simplicity. 

Equations (97) through (99) are in the platform frame representation. It 
is straightforward to convert these equations into the laboratory frame 
representation by substituting Eqs. (91), (93) and (94). The relative 

A A 

uncertainties between the two frames, 6p , 5n A and 65, then constitute new 

+n 

error sources. Thus, it can be shown that Eqs. (97) become 


A A 


A A 


g^(t) = (n + 5n A ) • g'(r + 6p,t) + 6n_ • f'(r + Sp,t) • (n + 65) 

+n 


A A 


j- 5n_ • a x (n + 5£) + 6g , 

4 c,t 


(100a) 


g'(t) = 6n_ • g'(r + 5p,t) + 5 (n + 6n A ) • F'(r + 6p,t) • (n + 55.) 
fl +n 


A A 


-5 6n A »axn + 6g . (100b) 

+5 d ,5 

Here g f and are the specific force vector and the specific force gradient 

tensor defined in Eqs, (4) and (8); 6g and 6g are the finite size terms 

c,t d,t 

given in Eqs, (A. l). Equation (100b) shows that the "sensitive axis of the 
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gradioraeter" can be defined by the average of (n^ + and A. 

B. General Description of Errors 

The accelerations in Eqs. (100) contain two types of error sources. The 
first type is intrinsic error sources due to angular motions as given in Eqs. 
(4) and (8), the second type is in geometrical metrology. Further error 
sources arise from non-ideal behaviors of the mechanical and electrical 
springs, scale factor mismatch, incorrect calibration, and residual coupling 
of the gradioraeter to temperature and electromagnetic fluctuations of the 
environment, etc. Equation (45) can therefore be generalized as 


i d (o)) - Hg i (o>) g c (w) + H^ i (o)) g d (w) + Hg?(w) G cc («) + H^(io) G dd (o)) 

+ Hq^(w) G cd ( W ) + H^Go) T c (co) + Hj i (o ) ) T d (o>) + ..., (101) 

where G cc (w), G dd (u>) and G cd (to) are the Fourier transforms of 

G cc (t) = [g c (t)] 2 , (102a) 

G dd (t) = [g d (t)] 2 , (102b) 

G cd (t) = (102c) 

and T c (w) and T d (w) are the Fourier transforms of the common and differential 
temperature fluctuations, T c (t) and T d (t), over the two transducers, defined 
in Eqs. (A26). The coefficients in front of these functions represent the 
transfer functions for the signal variables that they are multiplied with. 
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The centrifugal acceleration and the linear and angular accelerations, which 
are coupled to the gradiometer due to its metrology errors, appear in Eq. 
(101) through errors in g c (w) and g ( j(w) that they cause. 

Notice that Eq. (101) contains the original signal variables, g c (<o) and 
g ( j(u)), which are identical to g^(to) and g^(to) due to the relationships: 

g c (t) = g^(t) - n • g E (r) , (103a) 

g d (t) = g^(t) , (103b) 

which follow from Eq, (96). Equations (103) imply that the errors in the 
unprimed variables are given by the errors in the primed variables : 

6g c (t) = 5g^(t) , (104a) 

6g d (t) - 6 g< j(t) . (104b) 

It is also clear that the Fourier components of the errors are also 
identifical between the two sets of variables* Equations (100) and (103) can 
be combined to obtain the unperturbed acceleration signals : 

g co (<o) = g^ Q (w) = n • [g(r,to) - a(r,u>)] , (105a) 

g do (<o) = g do (u)) = * n * tr(r,d)) — C (to) ] • n , (105b) 

where g(r, to), a(r, to), r(r,to) and C (to) are the Fourier transforms of 
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g(r,t) = - $4>(r,t) - - $4> p (r,t) - ^(r) , 


(106a) 


a(r,t) - a p (r,t) + $ E x x + r)] , (106b) 

f(r,t) = - ^<|>(r,t) = - ^<|> p (r ,t) - ^^(r) , (106c) 

C (t) = £(t) £(t) -Q 2 (t) I . ( 106d) 

Here <J» E (r), Q E and R E are Earth's variables introduced in Eq. (28); <{> p (r,t) 

, *> >> . 

and a p (r,t) are the gravitational potential and the linear acceleration 
measured relative to Earth; and C (t) is the "centrifugal accelerator tensor.” 
The dyadic notation is used in Eqs. (106) to represent rank— 2 tensors. The 
first term in Eq. (105b) is A times the gravity gradient signal r nn (r,oj) to be 
measured. 

It is useful to divide Eq. (101) by A H^w) to obtain an equation of the 

form: 

i i d (w) . 

^nn^’^measured = J “d “ = Ir n n (r true + A J 6g d,i (w) » (107) 

where 6g d ^(<*)) is the equivalent error in the differential acceleration g d (u)) 
due to the i-th type error source. Derivation and discussion of 5g d i for 
various error sources are given in Appendix* Equation (107) can be converted 
into another convenient form: 

6r nn <?.«) - £ r to (J, u ) + £ e^to) * £ t a M a (?,») 

x,m X 

+ | e nA^ Q P,A^ W ^ + £ n^ V w) + ••• » (108) 
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in which dominant error sources have been shown explicity. Here the 
laboratory has been assumed to be stationary with respect to Earth and i5 p (o>) 
is the angular velocity of the platform in the laboratory frame. The first 
two error coefficients can be identified from Eqs. (A6b) and (A19) : 

f ' 'W W 6 " - + 6 *>i * < 109a> 

+n 

e nim (w) = " 8 An 6 mn 6a C (al) " (1 ~ 6 *n> W 6n - + ’ (109b) 

LIA.UI +n 

where 6a r (u) and 6o c (w) are the (dimensionless) calibration error coefficients 
for r(r,w) and C (w) , respectively. It is assumed that the zeroth order part 
of C (w) has been measured independently and removed from the gradiometer out- 
put. Otherwise, 6 o c (w) must be replaced by unity. The common acceleration 
error coefficient is read off Eq. (A21): 

e n* ((0) = " 6 in I " (1 - 6 Jln> J » (109c) 

where 5ag(co) is the (dimensionless) scale factor mismatch between two 
constituent accelerometers. The coefficient for the first-order angular 
velocity induced error is obtained from Eqs. (A10) and (A12): 

&<“> ■ <1 - *»,„> - j5f 0 - 6 kn> E kim «"->k 

+ 3“ l (1 - «kn> E k*n <«" ->k . < l09d > 

k +1 

where 5!$ p = £$ p has been assumed and is the totally antisymmetric tensor 

of rank 3 . 
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Unlike in some conventional gravity gradimeters, matching of proof masses 
ana spring constants are not necessary in the superconducting device because 
the balance is achieved by tuning persistent currents. However, the alignment 
of the sensitive axes of the proof masses ^s_ still important. 

The temperature error coefficient can be determined from the results of 
Section H of Appendix. Depending on the relative magnitude of the two 
competing thermal effects, the error coefficient can be written as 

e n (w) = I h Tc (a>) » ( 109e) 

where h^w) is given by Eq. (A42a), or by 

e n M " I 8ff S<«> S E cos 9 n eTF) S ’ (109e,) 

where E(T) is the Young's modulus of the spring material at temperature T. 

In Paper II, we discuss the observed values or limits of these major 
error coefficients. Once these error coefficients have been determined ex- 
perimentally, the relevant dynamic variables can be measured simultaneously 
along with the gravity gradient by independent instrumentation and the respec- 
tive errors can be compensated for by using the error model derived in this 
section and Appendix. In a three-axis gradiometer, it is possible to suppress 
several error sources by using the geometrical properties of the device [18], 
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VI. SUMMARY AND CONCLUSION 

A gravity gradiometer is a fundamental instrument which can separate 
gravity from frame accelerations. A sensitive tensor gravity gradiometer 
needs to be developed to carry out precision tests on gravity as well as for 
applications in gravity survey and inertial navigation. Since a large common 
mode background has to be removed in gradiometry by differencing signals at 
two or more proof masses, extreme stability is required for the sensitive axis 
orientation and for the scale factors of the component accelerometers. We 
have shown how these challenges can be met at low temperatures by utilizing 
the stability of persistent currents, the enhanced mechanical stability of 
materials and the sensitivity of SQUID amplifiers. An accompanying penalty is 
of course the inconvenience of having to keep the instrument in liquid helium. 
The liquid helium environment, however, can be used further to isolate the 
gradiometer from the fluctuations in the ambient temperature and electro- 
magnetic fields. 

In this paper, we have analyzed a superconducting current-differencing 
gravity gradiometer. Complete dynamical equations have been derived from the 
first principles. Transfer functions, common mode balance conditions and 
procedures as well as a rather extensive error model have been developed from 
these equations. The complicated differential instrument has been reduced to 
an equivalent accelerometer which simplifies the noise analysis and the 
electronic control of the device. In Paper II, we report the performance of 
the superconducting gradiometer in the laboratory, which verifies details of 
the theory developed here. Although our analysis has been confined to a 
specific instrument reported in Paper II, the theoretical methods presented in 
this paper will be a useful guide in analyzing other superconducting inertial 
instruments that may be developed in the future. 
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APPENDIX: DERIVATION (F VARIOUS ERROR COEFFICIENTS 


A. Finite Size Effects 

There are two error sources related to the finite dimensions of a 
practical gradiometer: (1) the finite volume of the proof masses and (2) the 
finite baseline of the instrument. The first effect can be minimized by 
choosing a nearly spherical geometry for the proof masses, whereas the second 
effect is reduced when the baseline length l » |t| is shortened. The finite 
baseline effect has been taken into account in Eqs. (95) by the Taylor expan- 
sion in 1. The finite volume effect has been ignored in the previous analysis 
because of its high order nature. The departure from a spherical geometry 
gives a nonvanishing quadrupole moment to each proof mass or a small octupole 
moment to the gradiometer. Since this will couple at best to the fourth order 
derivatives of 4>(p,t), which is generally small except when the source is 
extremely close to the gradiometer, this error can be ignored in most situa- 
tions. Clearly, the finite volume error is always smaller than the finite 
baseline error. 

The dominant gravity error terras arising from the finite baseline are 
obtained from Eqs. (97) and (98). 

2 

5g (t) - - (A • $) 3 <J>(p , t ) , (Ala) 

c, I 

3 

sg .(t) = - 4 t- u • <e(p,t) . (Aib) 

d, $ 

When <J>(p,t) varies with a characteristic length R, these errors become 
0(Ar/1* ) of the signal. Therefore, the finite baseline error can also be 
ignored when A/R « 1 . 
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B. Misposition of the Gradioneter 


The acceleration errors due to the misposition of the gradfometer are 
obtained from Eqs. (4) and (100): 

6g ^(t) = - n • [^<fr(r) + ^ - Q^I] • 6p - n • a x 6p , (A2a) 

c,p 

6g ^(t) = - i (n • $)[n • ^<t>(r) • 6p) , (A2b) 

d,p 

to the leading order in 5p. These errors are usually negligible because of 
the higher order gradients involved. The displacement 6p(t), produced by the 
linear acceleration a(r,t), modulates the large dc gravity of Earth and 
produces errors that compete with ac signals: 


6g ,<t> - 0 (Mil g ) , 

c,p e 

(A3a) 

6g ,U) - 0 (Mil ir ) , 

d,p E 

(A3b) 


where gg = 9.8 m s - ^, r E = 3.1 x 10^ E and Rg = 6.4 x 10^ m are the vertical 
gravitational acceleration, the vertical gravity gradient and the radius of 
Earth. It is clear that these errors can be ignored in general because 
6p ( t ) / Rg « 1 « 

C. Misalignment and Misorientation of the Gradioneter 

A A 

The acceleration errors caused by the misalignments , 6n_ and 6n A > and 

a A " 

the mlsorientations , 6n A , and 51, are obtained from Eqs. (100): 

+n 

6g a *(t) = 6n a • g'(r,t) + j (6n_ + 51) • f'(r,t) • n 
c,n+A +n 
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(A4a) 


X + 

-r 5n_ • a x n , 


6g a aU) = 6n_ • g'(r,t) + A(6n A + 6A) • f’(r,t) • n 
d , n+A +n 


A y A 

- I 6n , • a x n , (A4b) 

+JI 


to the leading order in 6n_, 6n A » 6n A . It follows from the definitions (92) 

a +A A+n 

- (94) that, in the limits |6n |, |6n A |, |6n A |, |6A| , « 1, 

+A +n 


6n_, 6n a , 6n a , 6A 1 n . (A5) 

+A +n 


Hence these errors cause the cross components of the common acceleration and 
the acceleration gradient to couple to the gradiometer. Written out 
explicitly, Eqs. (A4) become 


5g A A (t) = - 5n a • $<l>(p) + a(p)] - j (6n_ + 6A) • [$$<Kp) + • n 

c , n+A +n 


A A A. 

- 7 " 6n x n 
4 


a » 


(A6a) 


6g a A (t) = - 5n_ • $<J>(p) + a(p)] ~ A (6n A + 6A) • l^$<|>(p) + * n 

d , n+A +n 


A A ^ 

-A6nAxn»a. (A6b) 

+A 

In Eq. (A6b) , the misalignments within the gradiometer are expressed in 

A 

terras of the misalignment 6n_ of the accelerometer axes with respect to each 

A 

other and the misalignment 6n A between the average accelerometer axis 

a a +A A 

(n^ + n^)/2 and the baseline direction A. These errors cause coupling to 
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cross components of linear and angular accelerations which can be large com- 
pared with acceleration gradient signals of interest • The combined error 

A A 

6n * + appearing in the second term of Eq. (A6b) represents the mis- 
+n 

orientation of the gradiometer sensitive axis relative to the theoretical 

A 

direction n fixed to the laboratory frame, which is used for interpretating 
the data. This misorientation produces coupling to cross gradients. 

The coupling to the gradiometer from the dynamic variables of the plat- 
form, a(p,t), &(t) and a(t), are revealed in Eqs. (A6). Since X and 6n^ are 
time-varying in general due to <$(t), the large dc gravity of Earth, <|> E (p), 
will fee modulated to produce errors at the frequency of interest. This 
important error mechanism will be studied in detail in the next section. 

D. Angular Motions of the Platform 

The centrifugal acceleration appears without attenuation in Eq. (105b) in 
direct competition with the gravity gradient signal. This problem is funda- 
mental in any second-order gradiometer, in-line or cross component, as is 
evidenced by Eq. (9), and therefore does not depend on the particular design 
of the instrument. This makes the attitude control or detection as the most 
formidable task in precision gravity gradioraetry • One can easily compute the 
required attitude rate (l5) accuracy from Eq. (105b) for a given sensitivity of 

r. 

When the instrument is in a laboratory rotating at an angular velo- 
city S$ Q (t) with respect to an inertial frame, it is convenient to define 

3(t) = $ (t) + l,(t) . (A7) 

o r 

If uncertainties in $ Q (t) and i!$ p (t) are 6^ Q (t) and 6$ p (t), respectively, then 
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the error in the differential acceleration due to centrifugal acceleration is 
obtained from Eqs. (105b), (106d) and (A7) as 


6g (t)—2A [{n»(!$ + ^ )}{n»(6^ +6i$ )}-(<!i + )»(5<$ 

Or OroPo 


+ 6 & p )]. 

(A8) 


If the gradiometer platform is moving with a velocity v(t) with respect 
to Earth and the gradiometer is in a local geographic orientation, then ^ Q (t) 
is given by 

^(t) = + \ 2 x » (A.9) 

where ^ and are the angular velocity and radius vectors of Earth, intro- 
duced in Eq. (28). In the case when the platform is stationary with respect 
to Earth; i.e. v(t) = 0, Eq. (A8) reduces, after Fourier transformation, to 


6g (w) = -21 [(n • ){n • (to)} - (w)] , (A10) 

d,C E P E P 


where we have assumed the condition Q p « = 7.27 x 10 - ^ rad s - *, which is 
not difficult to satisfy for a stationary platform. 

In a terrestrial laboratory, the modulation of Earth's gravity bias by 
the tilt of the sensitive axes is another important error mechanism, as was 
pointed out in Section C of this Appendix. Since the common mode errors in 


Eq. (A6a) produce second order errors in the differential signal when 
multiplied by the common mode balance error according to Eq. (51), these 
errors will, in general, be dominated by the errors in Eq. (A6b), which is 
raultipled by H^, . 
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When the platform moves at an angular velocity ^p(t) with respect to the 
terrestrial laboratory, the misalignment and misorientation vectors are 
modulated according to 


^ = ^p(t) x 1 , 


(All) 


v A A A A 

where A stands for 6 n_, 6 n A , 6 n A and 6 JL Equation (A 6 b) becomes, upon 

+£ +n 

Fourier transformation, 

A ^p(k)) . v A A ^ (to) . A 

5g a a(w) = - 6 n • — r- — x g (r) - £( 6 n A + SA) • - 77 - — x [r (r) • n] 

d,n+Jl J0) E +n JO) E 


A A 

■Hn , x n • [jw 0 (co) ] , 
+1 


(A12) 


where gg = Tg = -^<J>e and a(u) = jw ^p(u) were substituted. The terra 

A 

arising from the angular modulation of 6 n A produces a second order effect in 

* + * 

Q p (w) and therefore has been dropped. In the case when 


| 6 n_| > | 6 n A + 61 1 x 


IT 


E 


(A13) 


+n 


the second term in Eq. (A12) can be ignored. The third terra becomes 
negligible compared to the first term at sufficiently low signal frequencies: 


(«) 


< (% l/2 



(A14a) 


+1 

Comparing Eq. (A12) with Eq. (A10), one finds that the time-varying tilt noise 
is the dominant source of angular motion induced errors in the frequency range 
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(A14b) 


St? a 

" ‘ is; l 6n -l ■ 

even when 6^ p * $ p , the worst case for the centrifugal acceleration error. 

In Paper II, we will see that the angular motion noise given by Eq. (12) is 
indeed the dominant error mechanism for our prototype superconducting gravity 
gradlometer. 

The foregoing discussion clearly reveals the advantage of operating a 
sensitive gradlometer in space where the gradlometer platform is freely 
falling. If an Earth-pointing reference frame is chosen so that ^ Q (t) coin- 
cides with the orbital angular velocity of the satellite, then Eqs. (A8) and 
(A6b) reduce respectively to 


[6g _(w)] 

d,C 


space = ~ 2 * * V {n * * 625 p (w)] ,(A15) 

^ n (w) 


[6g a a(«)] “ -A (6n A + 6A) • x [f w (r) - ] • n 

d,n+A space +n 




o o 


A A . 

-A 6n A x n • [jio Q (co) ] , 
+A 


(A16) 


where r p (r) is the gravity gradient tensor of (spherical) Earth. Comparison 
of Eqs. (A15) and (A16) leads to a conclusion that the centrifugal accelera- 
tion error could be dominant in most practical cases. 

One can see from Eq. (A8) that 6g (t) has a minimum when the gradio- 


d,C 

meter is in an inertial orientation; i.e., = 0. The worst situation for 

the centrifugal acceleration error is the case when the gradlometer is rotated 
at a frequency high compared to the signal frequency for heterodyne detection 
of the gravity signal. In this case, is the spin angular velocity which is 
large compared to One advantage of the superconducting gravity gradio- 
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meter is that the heterodyne detection is not necessary , unlike in some con- 
ventional devices, because of its low drift and excellent low frequency noise 
characteristics. The superconducting gradiometer therefore permits complete 
freedom in orientation: inertial, earth-pointing, or spin-stabilized [26]. 

Hence, an optimum orientation, which minimizes the overall error budget, is 
available to the superconducting instrument although the angular motion in- 
duced errors are fundamental. 

In Eqs. (A10) and (A15), we considered only the first order errors in 
65p(o)). This is justified when 5<Sp(w) is bandwidth limited. The Fourier 
transform of [6i5p(t)]^ down-converts the centrifugal acceleration noise from ' 
high frequencies to the vicinity of dc, the signal bandwidth. A low-pass 
filter for angular vibrations of the platform is therefore needed to suppress 
this error. 


E. Scale Factor Errors 

There are two types of scale factor errors in a gradiometer: (1) a rela- 

tive error, which comes from the mismatch of the scale factors of the com- 
ponent accelerometers, and (2) an absolute error, which is the error in the 
calibration of the gradiometer transfer function. Since the scale factors are 
determined by persistent currents in the superconducting gradiometer, extreme 
stability is expected in the scale factor match and calibration. The passive, 
iterative, common mode balance procedure described in Section HID allows in 
principle an arbitrarily precise match of scale factors in one direction inde- 
pendent of signal frequency. In practice, however, H^(u)) is not precisely 

zero in Eq. (101), thus producing a common mode rejection error, and H^.(w) 

gi 

also contains a calibration error 6H^. (w). 

gi 

Let us define the (dimensionless) coefficients for scale factor mismatch 
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and gradiometer calibration error by 


6o s (w) = 


H*" . (to) 


6h1. (u>) 


6a r<“> = -/r: — i » 


V w) 


(A17a) 


(Al7b) 


and for the centrifugal acceleration calibration error by Scr c (co), which is 
usually different from Sap(oj) because independent instruments, such as gyro- 
scopes, are used to measure and compensate the effect. The error arising from 
the use of an inaccurate value of R. in the recovery of T 1 nn ( r > t ) from gj Q (t) in 
Eq. (105b) is represented by the second term in Eq, (A17b), With the aid of 
Eqs. (105), the equivalent differential acceleration errors due to scale 
factor mismatch and calibration errors can be written as 


6g d)S < w > 


6g d>c (u>) = * n 


6o s (a3) n • [g(r , oj) - a(r,to)] , 
[Sa^(o)) f(r,w) - 5a c (io) C(w) ] 


n . 


(A18) 


(A19) 


The platform motion terms usually dominate over the ac gravity terms in these 
equations. The second terra in Eq. (A19) constitutes an additional rotation 
induced error which must be added to the list considered in the previous 
section. The linear acceleration terra in Eq. (A18) will be discussed in the 
following section along with other translation induced error sources. 

Now, we briefly discuss 50g(w) which arises from the failure to apply the 
wideband balance procedure. If the balance is attempted at a single frequency 
(i>k, Eq. (48) must be violated in general. Substituting Eqs. (46) into Eq. 
(A17a) and expanding it in a Taylor series, one finds 
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2 K ~ w) ' J b (w 2c 


2 \ 
' W lc ) 


f 2 2 W 2 2 n 

(w. - 0). )(0)„ - 0 ). ) 

1C O ZC D 


6a s (Wb) 


(A20) 


where is a residual balance error at & b . Although not essential, a 

highly symmetric gradiometer is convenient because oj 2c * cj 1c and, therefore, 
the wideband balance can be achieved with only two persistent currents Ij and 
I 2 “ 1^ while keeping ij = i 2 = 0. 


F. Linear Notions of the Platform 

Although several sources of linear motion induced errors have been shown 
in the previous sections, it is instructive to combine these effects here. 

The linear velocity dependent Coriolis force terra has been dropped in Eqs. 

(95) on the assumption that the proof masses are confined to move along the 

A A 

linear directions n^ and n 2 « A non-vanishing residual compliance in the 
directions perpendicular to the sensitive axes would allow the proof masses to 
respond to the transverse components of accelerations, thereby contributing to 
a velocity dependent error. It is therefore important to design the 
suspension spring with high stiffness for all undesired degrees of freedom 
while obtaining a very weak spring in a linear direction. We will see in 
Paper II how this condition is met in the actual design of the superconducting 
gradiometer. 

It has been shown in Section B of this Appendix that the modulation of 
Earth's gravity by a time-varying displacement error 6p(t) can be ignored 
because of the higher order gradients involved. The dominant mechanisms which 
convert the linear vibrations of the platform into gradiometer errors are 

A. 

therefore (1) the sensitive axes misalignment 5n_, which couples the cross 
component accelerations, and (2) the scale factor mismatch 6cf s (<d), which 
couples the in-line component acceleration to the gradiometer output. 
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Combining proper terms in Eqs. (A6b) and (A18), one finds 


6g (w) 
d,a 


[6n_ + 60 g(a)) n] 


a(r,w) . 


(A21) 


The two terms in this equation are related by the condition of the common mode 

balance. The balance is obtained in the real instrument by requiring 

6g (w.) = 0 for a given applied acceleration, say, in the z-direction: 

° > b 

d,a 

a(r,w b ) = a(w b ) z . (A22) 


Substitution of Eq. (A22) into Eq. (A21) leads to 

A A 

6n • z 

6a s (u) b ) = a — + ScJ sO» (A23) 

n • z 


which can be substituted, in turn, into Eq. (A20) to obtain 6o s (w). Notice 
that, even with a wideband balance, 5<Jg(w) ** 6<Jg(a)^) ^ 0 due to the 
misalignment of the sensitive axes and the residual balance error 6 o s q. The 

A 

impossibility of obtaining a balance with an applied acceleration normal to n, 
or (n^ + x\^)/2 to be rigorous, is indicated by the divergence of Eq. (A23), 

A A 

which occurs when n • z = 0. 


Gm Scale Factor Nonlinearity and Dynamic Range 

A nonlinear behavior of accelerometer scale factors arises from departure 
of the acceleration response of the mechanical suspension springs from the 
linear behavior predicted by Hooke’s Law, higher order terras in the modulation 
of the pancake coil inductances, L a and L^, higher order terms in the magnetic 
force F em , and, finally, nonlinearity in the current-to-voltage transfer 


- 68 - 


function of the SQUID electronics. In Eqs. (14), we assumed perfect linearity 


ao i. unctions of the coll spacing. This as sump tic 


4 « *_ 

7 uiup tAWU J-O J.cai J.JT uu u 


rigorous because of the grain structure of the pancake coils and edge effects. 
The higher order terras in F EM , which are apparent in Eqs. (23) and (25), have 
been ignored in our attempt to linearize the circuit equations. Therefore, 
even if the mechanical springs and the SQUID electronics are constructed with 


sufficient linearity, the response of each accelerometer will be somewhat non- 
linear due to the Inevitable higher order response of the electrical circuit. 

In analogy to Eqs. (A17), we define the nonlinearity error coefficients 
h cc ( w )> and hcd^ a> ^ by dividing the nonlinearity transfer functions 

rC gI (w), H^(w) and H^(w), introduced in Eq. (101), by H^oo). Then the equi- 
valent differential acceleration error due to scale factor nonlinearity can be 
written as 


6 Sd,N< w > " h cc (w) G cc (w > + h dd (w) G dd< w > + h cd< w > G cd< u >* < A24 > 


Now, in most practical situations, the gravity terras are small compared to the 
platform motion terms in Eqs. (100) so that 

G cc (t) = [n • a (r,t) ] 2 , (A25a) 

G dd (t) = JL 2 [{n . $(t)} 2 - Q 2 (t)] 2 , (A25b) 

G cd^ b) = A [n • a (r,t) ] f{n • &(t)} 2 - Q 2 (t)]. (A25c) 

It is clear that the Fourier transformation of these nonlinear functions of 
dynamic variables will down-convert the wideband platform noise to the 
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vicinity of dc, where the gravity gradient signal is to be detected. There- 
fore, it is imperative that the vibration and jitter noise of the apparatus be 
low-pass filtered before they reach the gradiometer. Detailed analysis of the 
nooise down-conversion process will be given in Paper II. 

It is difficult to determine the nonlinearity error coefficients analyti- 
cally. They can however be measured from the departures of the transfer 

functions H c . (w) and H d .(w) from the linear behavior. In principle, the non- 
gi g 1 

linearity in the two component accelerometers could also be matched. This 
would not, however, reduce all three error coefficients in Eq. (A24) to zero. 

A more powerful and practical approach to the problem is the linearization of 
the system by means of an electromechanical feedback. A feedback force, which 
is equal and opposite to the detected external force , can be applied to each 
proof mass so that the total force, the "error signal" for the feedback loop, 
is reduced by the inverse of the feedback gain. The reduction of the driving 
accelerations limits the actual displacement of each proof mass to a small 
amplitude, thereby reducing ^g^^Cw). The error coefficients h(w)'s have 
therefore been effectively reduced. The average and differential feedback 
forces are direct measures of g c (t) and g d (t). The gradiometer thus measures 
the specific force instead of the resulting displacement. The negative feed- 
back comes with another important advantage : increase in the dynamic range of 

the device. This so-called "force rebalance" feedback has been successfully 
applied to conventional inertial navigation instruments [27] and gravity 
gradiometers [7], although it is yet to be applied to superconducting gravity 
gradiometers . 

One of the parameters responsible for the high sensitivity of the super- 
conducting gradiometer is the low damping coefficient of the electromechanical 
springs. The low dissipation is necessarily accompanied with high Q 
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resonances of the modes which tend to amplify the acceleration noise at the 
lesuuaiicc nequeiicles. An elegant way to take care of this problem is an 
active damping method [28] by which the Q's are effectively lowered without 
increasing the Brownian motion noise level. The gradiometer response is 
narrowband filtered around the resonances, phase shifted by 90°, and fed back 
to the proof masses to actively drive down the resonant motions. This "cold 
damping" has been successfully demonstrated in superconducting gravity gradio- 
meters (13, 14]. 

H. Nonmechanical Noise of the Environment 

It has been pointed out in the Introduction that superconductivity can be 
used to make a nearly perfect electromagnetic shield. In a practical super- 
conducting shield with impurities, the Meissner effect does not exclude the 
magnetic field completely, but instead "freezes" some trapped magnetic field, 
thus providing a "perfect" shield against time-varying fields. One can 
combine high permeability shields with superconducting shields to attenuate 
both dc and ac electromagnetic fields as well as radiation very effectively. 
Therefore, the electromagnetic susceptibility of the superconducting 
gradiometer can be made negligible. 

The thermal fluctuations of the environment could also be shielded by 
immersing the apparatus in superfluid helium, the Bose condensate phase of He^ 
below 2.17 K. The infinite heat conductivity and the large heat capacity of 
the superfluid provides a stable, gradient free, thermal environment. Even 
normal fluid helium can provide an excellent thermal environment provided its 
vapor pressure is regulated. When the attenuation of the temperature fluctua- 
tions of the environment is not sufficient, one will have error terms that 
couple the temperature noise to the gradiometer output, as shown in Eq. (101). 
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The fluctuation in the gradiometer temperature can cause error signals through 
its interaction with (1) the mechanical part and (2) the superconducting 
circuit. 

In analogy to Eqs. (40), the temperature modulations of the two accelero- 
meters, T^w) and T 2 (w), are expressed in terms of their common and 
differential temperatures , T c (w) and T d (w) : 

TjCu) = T c (w) - j T d (w), (A26a) 

T 2 (u) - T c (w) + j T d (u). (A26b) 

The mechanically coupled temperature-induced error could arise simply from the 
thermal expansion of the baseline: 

5M- «<T) T c (co), (A27) 

where a(T) is the thermal expansion coefficient of the gradiometer body at the 
ambient temperature T. However, a(T) is extremely small for solids at liquid 
helium temperatures so that the dimensional change 5 A(u) is negligible in any 

A 

practical situation [14]. For the same reason, alignment changes 6n ± (w) 
resulting from temperature fluctuations are negligible, demonstrating the 
mechanical stability of the cryogenic gravity gradiometer. 

There is a more subtle coupling mechanism of the temperature noise 
through the mechanical part of the system. It is through the temperature 
dependence of Young's modulus E(T) of the mechanical spring [29]. Since the 
stiffness of the suspension spring changes as a function of temperature, a 
temperature fluctuation will cause a displacement modulation when the spring 
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is under stress either by gravity or by an unbalanced magnetic pressure. This 
can be seen iti Eq. (25) for an accelerometer. mien n - £„ ( r) r 0 so that * 

Ci O 

2 

0, the temperature induced modulation of produces a first order effect in 
displacement* It is convenient to move this term to the righthand side of Eq. 
(25) and define an additional effective driving acceleration: 


6g T (t) = - [wJ(T + 6T) - u£(T)) x q . 


(A28) 


The proportionality of the spring constant im. to E(T) can be used to rewrite 

M 

Eq. (A28) as 


6g T (t) = - w^(T) x Q H ^T(t). (A29) 

Substitution of Eq. (27) into Eq. (A29) and Fourier transformation leads to 

8 T (») ■ [« E COS 9 n + i (I + i jS-pjT ‘Xliliyr T< “>- <A30) 

a b 

It is clear from this equation that this effect can be eliminated by choosing 
i such that the quantity in the square bracket vanishes; i.e., by compensating 
the gravitational force with a magnetic levitation force. 

In the actual operation of the gradiometer reported in Paper II, i = 0 
was chosen so that the magnetic pressure term in Eq. (A30) was negligible. 

The temperature effect on the gradiometer can be found then by adding 

8>> “ Se cos 0 n e(Y) T c (u) » (A3 la) 

*d (w) “ «E cos 9 n f(TT if T d (w) (A31b) 
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to g c (w) and g d (co) in Eq. (41). It is assumed here that the two suspension 
springs of the proof masses have an identical Young's modulus E(T). The 
temperature-induced current output can be written as 


i^(w) = H^Cu) T c (oj) + H^Cu) T d (a)), 


(A32) 


where 


H Ti (w) = H gi (a)) g E cos 9 n eTW fl ’ (A33a) 

H Ti (u) = H gi (a)) g E cos 9 n * (A33b) 

Notice that the common temperature fluctuation T c (w) is balanced out to the 
same degree as the common acceleration . For a properly balanced gradioraeter , 
therefore, this effect is expected to be negligible as long as i = 0. 

The second effect of the temperature noise, which acts on the supercon- 
ducting circuit directly, comes from the temperature dependence of the 
’’penetration depth” X(T) of the magnetic field in the superconductor: 


X(T) = 


X(0) 

[1 - (T/T q ) 4 ] 1/2 


(A34) 


where T 0 is the critical temperature for the superconductor. For niobium, 
X(0) = 5.0 x 10 -8 m and T Q = 9.2 K. The variation of X(T k ) of the proof mass 
and pancake coils due to a temperature fluctuation 6T k (t) of the k-th 
acceleration transducer. 


6X(T k ) = X(T k + 5T k ) - X(T k ) * ^ 6T k , 


(A35) 
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modulates the superconducting inductances of Eq. (14) according to 


W c > ” Mca + + S) 6X(T k ) + A x£(t) , (A36a) 

yt) - Lfcb + Ml + O 6A(T k ) - A x£(t) . (A36b) 

Here, 


0 < C < 1 (A37) 

is the modulation efficiency of the penetration depth of the pancake coils. 

Approximation is made here that the sensing inductances are reasonably well 

matched so that distinction of C for each of these coils is not necessary. In 

Eqs. (A36), x^(t) is the dynamic displacement produced in self-consistent 

response to the current modulations resulting from the temperature-induced 

inductance modulations. Given the penetration depth modulations, the dis- 

T T 

placement response x^(t) and the current response ij(t) are solved from the 
dynamical equations of the gradiometer. The dynamical equations and their 
solutions are modifications of those given in Section III. 

Modifying Eqs. (18) and (19), the series and parallel inductances of the 
sensing coils and 1^ now becomes 

(A38) 
(A39) 

Flux quantization through the loops L ka + gives 


W*> " \b + 2 A(l + 6X < T k>> 

2 2 

Ifcp(t) = Lfcp + A <1 + C) 5\(T k ) - ^ A x£(t). 


- 75 - 


I k (t) - Ik 11 - 


2A(1 + C) 6A(T k ) 


(A40) 






Substituting Eqs. (A38) - (A40) into the flux quantization condition for the 
loops 1^(0 + L Q , which is generalized from Eq. (32), and taking Fourier 
transformation, one finds a modified version of Eq. (42): 


xj(to) + «J(1 + C) 6X(T k (u)) = ~ J id (w) + (_1)k 1 c (a)) » (A41) 


where the coefficients are defined by 


T _ ^ka ^kb ^kp **k 
“k = 




K * 


(A42) 


The homogeneous part of Eq. (25) generalized for k = 1, 2 is then solved using 
Eqs. (A38) - (A40) again, resulting in 


* 2.2 

, 2 . 2 A \ 


T 1 T k^^kT 

■> + ? #■> * (-D -5* i c r C») 


* ^ + v 1 

“khcs 


(1 + c) 6X(T k (w)), 


(A43) 


which are identical to Eqs. (41) with the driving gravity signal replaced by 

the 6X(Tj c ) term on the right hand side. 

The relationships between the currents and 6X’s are obtained by 

eliminating x k from Eqs. (A41) and (A43). The resulting pair of equations are 

the same as Eq. (43) on the lefthand side and have equal coefficients of 
T 

i c (u)) due to the common mode balance condition of Eq. (48). Therefore, the 
T 

i terras are readily eliminated, yielding 
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ij(b>) - ^(wXl + r,)[ <x£ (u 2 2x - (O 2 ) 6\(T 2 (w)) - " w 2 ) SMT^a)))], 

(A44) 


where 


2 _ 2 
^kX ‘ “kM 


2 2 
A ^ 

+ — ~ + 


m, 


vh. 


S 


* v \ 
\ \ 


(A45) 


The temperature induced output current can now be expressed in the form of Eq. 
(A32) with the identification: 


H Ti (w) = H gi (w) h Tc (a,) > 

(A46a) 

H Ti (t0) “ H gi ((o) h Td (w >» 

(A46b) 


where the temperature error coefficients are obtained from Eqs. (A26), (A35) 
and (A44) as 


h Tc (w) = [a 2 (o) 2 ^ - w 2 ) - - to 2 ) ] ( 1 + O » (A47a) 

h Td (u>) ■ y “ w 2 ) + aJXu) 2 ^ ~ « 2 )](1 + . (A47b) 


In terras of these error coefficients, the equivalent differential acceleration 
error due to temperature sensitivity is 


Sgd >T ( w > 


It is apparent from Eqs, 


= h Tc (w) T c (<o) + h Td T d (w). (A48) 

(A42) and (A47a) that the effect of T c (w) is 
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partially balanced by the match between pancake coil inductances, I^ a and L^ b , 
and by the common acceleration balance condition (48). On the other hand, 
T d (a)) is usually negligible because the good thermal conductivity of the 
gradiometer body keeps the entire instrument in thermal equilibrium. 

It has been pointed out [18] that an exact common temperature balance 
could be achieved by adjusting the fourth persistent current Z 2 with respect 
to 1^ , I 2 and ij that are used for the wideband common acceleration balance. 
The adjustments are iterated such that both the common acceleration balance 
condition (48) and the common temperature balance condition, 

a 2^2\ “ 01)2 ^ “ a l^ w l\ ~ w 2 ) » (A49) 

are simultaneous satisfied. In practice, it will be easier to couple a 
separate superconducting loop, which senses only temperature, to the SQUID and 
adjust the persistent current in this loop to obtain the temperature 
balance. The advantage of the latter scheme is that the temperature and 
acceleration balances can be performed independently. 

A slowly varying temperature of the environment, if uncompensated, can be 
an important source of a dc drift of the gradiometer. Such a drift, however, 
does not produce a random walk of the output, but the error is bounded because 
the dc level of the output is locked to the temperature of the gradiometer 
which is self-regulated to a large extent by the liquid helium itself. 
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PART 2 


DEVELOPMENT OF A PROTOTYPE SINGLE-AXIS 
SUPERCONDUCTING GRAVITY GRADIOMETER 


SUPERCONDUCTING GRAVITY GRADIGMETER FOR SENSITIVE GRAVITY MEASUREMENTS: 

II. EXPERIMENT* 


H.A. Chan, M.V. Moody and H.J. Paik 
Department of Physics and Astronomy 
University of Maryland, Collge Park, MD 20742 


A sensitive superconducting gravity gradiometer has been constructed and 
tested. Coupling to gravity signals is obtained by having two superconducting 
proof masses modulate magnetic fields produced by persistent currents. The 
induced electrical currents are differenced by a passive superconducting 
circuit coupled to a SQUID (Superconducting Quantum Interference Device). The 
experimental behavior of this device has been shown to follow the theoretical 
model closely in both signal transfer and noise characteristics. While its 
intrinsic noise level is shown to be 0.07 E Hz -1 ^ 2 (1 E = 10 -9 s -2 ), the 
actual performance of the gravity gradiometer on a passive platform has been 
limited to 0.3 ~ 0,7 E Hz ^ 2 due to its coupling to the environmental noise. 
The detailed structure of this excess noise is understood in terms of an 
analytical error model of the instrument. The calibration of the gradiometer 
has been obtained by two independent methods: by applying a linear accelera- 

tion and a gravity signal in two different operational modes of the instru- 
ment. This device has been successfully operated as a detector in a new null 
experiment for the gravitational inverse square law. In this paper, we report 
the design, fabrication and detailed test results of the superconducting 
gravity gradiometer. We also present additional theoretical analyses which 
predict the specific dynamical behavior of the gradiometer and of the test 
platform, and compare the results with experiments. 

*Work supported by NASA under contract No. NAS 8-33822. 
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I. INTRODUCTION 


The Equivalence Principle prohibits distinction of "gravity" from 
"acceleration" by a local measurement. The same principle, however, implies 
that a sensitive accelerometer , or a gravimeter , can be used as an approximate 
detector of gravity in an environment where the platform accelerations are 
small. Thus, gravimeters, which measure absolute and relative magnitudes of 
the gravitational acceleration vector, usually from a stationary platform, 
have been employed as geophysical survey instruments [1], This method of 
measuring gravity can be extended to a moving platform by independently 
determining and compensating for the platform accelerations with the aid of a 
geodetic navigation device such as the Global Positioning System (GPS) 
satellite network [2], The latter measurement is clearly of a global nature. 

An alternative, more direct, approach to a rigorous determination of the 
gravitational field is by means of a gravity gradiometer , which measures com- 
ponents of the tidal force or the Riemann tensor [3]. Here an almost local 
measurement is made by monitoring relative motions of proof masses separated 
by a short, but finite, baseline. With the torsion balance as the predeces- 
sor , several advanced types of gravity gradiometers have been developed 
recently [4]. 

The superconducting gravity gradiometer, to be discussed in this paper, 
has been developed with the primary motivation being to perform precise tests 
of gravitational theories in both terrestrial and space laboratories. A pre- 
liminary version of a null test of the gravitational inverse square law has 
been carried out with this instrument [5], While the instrument has served as 
a prototype for a three-axis superconducting gravity gradiometer [6] which is 
under construction, the system has undergone a thorough theoretical and ex- 
perimental analysis to test the physics of the instrument. The results could 
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easily be adapted to a more complicated device. 

In Paper I, a detailed theoretical analysis of the current-differencing 
gravity gradiometer was presented. This paper (Paper II) discusses the 
design, construction and test results of the actual instrument. These test 
results demonstrate that the physical device is described by the theoretical 
model in great detail. In addition, we discuss here the principle and 
operation of two elegant experimental techniques, which have been applied to 
the superconducting gravity gradiometer: (1) accelerometer mode calibration 

and (2) vibration isolation by means of pendulum action. Also, the error 
model of the instrument developed in Paper I is extended and applied here to 
the specific instrument configuration chosen. Paper II represents an 
extension of the experimental work reported in Ref. 7. 

The gravity gradient signal sought for is usually a small fluctuation 
over a relatively large background of the Earth’s gravity and the seismic 
activity of the platform. Therefore, precise mechanical alignment of com- 
ponents and extreme stability of scale factors, as well as high sensitivity, 
are required for a practical gravity gradiometer. Advantages of a super- 
conducting instrument in these regards have been pointed out in Paper I. In 
this paper, we show how quantum-mechanical properties of superconductors are 
utilized specifically., along with careful mechanical design, to meet the 
challenge of constructing a highly sensitive gravity gradiometer. 

Notations employed in Paper II follow those of Paper I, unless explicitly 
stated otherwise. 
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II. DESIGN AND FABRICATION OF THE SUPERCONDUCTING GRAVITY GRAD IOMETER. 

A. Gradloaeter Hardware 

The gradiometer consists of two accelerometers, each with a proof mass, a 
suspension structure, and sensing coils. The suspension structure consists of 
mechanical springs which confine the proof mass motion along the desired 
direction. 

A folded cantilever suspension has been employed to produce a suspension 
with a linear spring constant which is weak, in the direction of motion, yet 
relatively rigid for other degrees of freedom. This type of suspension was 
suggested by Prof. Daniel DeBra at Stanford University and has been employed 
in the displacement-differencing gravity gradiometer [8]. The design contains 
eight folded cantilevers in an 8m symmetry (to be defined later) to suspend 
the proof mass at the center (Fig. 1). Each cantilever consists of three 
joined diaphragm strips, of length Jt g and thickness t g , formed at the bottom 
of one radially oriented groove. The "fold” of the cantilever is at the 
radially outward ends of the strips. The radially inward ends of the strip 
are the ends of the spring. The middle strip of width b g holds the center 
moving mass. The other two strips of width y b g link to a bulk mass which is 
a solid continuation of the outer rim of the suspension. This design allows 
the cantilevers to undergo pure bending without stretching in order to achieve 
linearity in the spring constant. There are four folded cantilevers on each 
of two parallel planes with 2it/8 rotation plus inversion symmetry (8m 
symmetry). Thus, motion of the center mass is linearized along the 
cylindrical axis. 

The entire suspension structure is machined out of a single piece of 
niobium (Nb) in order to insure mechanical precision and to obtain a high 


- 86 - 


Symmetry: S 8v 
(5m) 

8: 2j rotation 
8 

m: Vertical 
Mirror Plane 



Proof 

Maes 



Direction 
of Motion 



Double 

Diaphragm 

Zero Length 

Springs 

Suspension 

Scheme 


inches 


Fig. 1. 


A cut-away view of proof mass suspension. 




quality factor of resonance. Eight grooves, four from each side, are first 
milled out, leaving a thick diaphragm at their opposite sides. The metal is 
then stress-relieved at 800°C. The grooves are then electric-discharge- 
machined (EDM) to thicknesses of (0.35 ± 0.01) mm. Next, the suspension 
structure is chemically polished with a mixture of hydrofluoric acid (HF) , 
phosphoric acid (H 3 PO 4 ) and nitric acid (HNO 3 ) in an ultrasonic cleaner. 

After the polishing, the Nb is partially recrystalized by annealing in a 
vacuum oven at 1100°C. The designed pattern of slots is then cut with a wire 
EDM. Finally, the finished proof mass suspension structure is chemically 
polished again. 

The center moving mass of the suspension is loaded with two mushroom- 
shaped Nb masses which screw into opposite ends to form one proof mass. This 
added load serves to increase the mass m^, lower the suspension frequency f^, 
and provide two superconducting planes which modulate the inductances of the 
sensing coils. The mechanical suspension frequency of the proof mass is de- 
rived from the formula for a loaded beam and is given by 



<2 *V 


3 

E b t: 


3 s s 

2 ' 


( 1 ) 


Here E g = 1.03x10^ N m - ^ is the Young's modulus of Nb and n g = 8 is the 
number of springs used. Also, t g = 0.35 mm, i g = 19 mm, b g = 2.5 mm, and m^ = 
0.40 kg, giving f fcM = 20 Hz. 

The sensing inductors are "pancake -shaped" coils wound out of thin (0.076 
mm in diameter) Nb wire in a single layer. The coil form is made of machin- 
able glass -ceramic, Macor (Corning Glass Works, Corning, New York), whose 
thermal expansion coefficient matches closely with that of Nb down to cryo- 
genic temperatures. A continuous length of insulated Nb wire is wound uni- 
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formly on the flat coll form in a spiral shape. Low viscosity, transparent, 
TRA“ CAST BB30Q2 epoxy ( TEA- CON, Inc = , Medford, Massarhnfipt t s) is used to bond 
the Nb wire. The coil winding procedure used here is a modification of that 
described in Ref. 9. 

Two such pancake coils are located next to the two plane surfaces of the 
proof mass in each accelerometer. To a high accuracy, the inductance for each 
pancake coil, 


L(t) = n Q A L + > (2) 

changes linearly according to Its average spacing from the superconducting 
ground plane, and is Insensitive to any rocking or transverse modes of the 
proof mass. The coil forms are mounted rigidly on Nb holders using Nb screws 
and GE 7031 varnish. The two Nb holders are mounted rigidly to the outer rim 
of the proof mass suspension using titanium (Ti) screws to form a shielded 
accelerometer. Ti is closely matched in thermal expansion coefficient with 
Nb. The spacing between the coil and the proof mass is adjusted by adding Nb 
spacers cut out of Nb sheet. Parallelism between the coil and proof mass 
surface is desirable to allow the use of a small spacing between them. 

Two accelerometers are mounted on the opposite faces of a precision Ti 
(Ti V4 AA6 alloy) cube to form a single-axis gradiometer (see Fig. 2). The 
surfaces of the cube have parallelism and orthogonality within 50 ppm. The 
same degree of parallelism is kept for the mating surfaces of all the gradio- 
raeter parts. The Ti material at the center cube is continued down to an 
aluminum-to-titanium joint. An aluminum (Al) base is tight fitted onto the Ti 
at room temperature. Differential contraction tightens the joint further as 
the assembly is cooled down. A similar joint is made at the upper end. The 




baseline of the gradiometer, which is defined by the distance between the 


centers of mass for the two proof masses, is A — 0.16 tti , 

For most of the tests reported in this paper, one triagonal axis of the 
mounting cube is aligned with the vertical, as is shown in Fig. 2, so that the 
sensitive axis of the gradiometer is tilted by an angle 


0 U = tan - * /2 


(3) 


from the vertical. This "umbrella orientation” has the advantage of allowing 
measurements in three orthogonal directions to be made by rotating the gradio- 
raeter around the vertical axis by 120° increments. The fact that the large dc 
gravity bias is matched along the sensitive axes is an important scientific 
reason to prefer this orientation when a careful three-axis measurement is 
called for. The above procedure has therefore been used for the null test of 
the gravitational inverse square law [5], which will be discussed in Section 
VI. 


B. Superconducting Circuitry 

The superconducting circuit for the gradiometer is shown in Fig 3. lq a , 
Ljlj, and I- 2 a , are pancake coil inductors for the two proof masses, raj and 
m 2 , respectively. S^ a , S^, S 2 a , S 2 ^ and Sj are current leads from the 
current supply outside the cryostat. Hj, H^g, H 2 , and are shielded heat- 
switch resistors. When currents are sent through these resistors, appropriate 
sections of the superconducting path are switched into their normal state. 

and R 25 represent extremely low resistance current paths (3xl0 - ^ Q). 
is the input inductor of the SQUID amplifier. Lj and Ljj are, respectively, 
the primary and secondary of a high coupling transformer. 
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The current Ij can be stored in Lj a + Lj^ by turning Hj on and then off 
while the current is sent across Sj a and S lb . The inductance L la + can be 
determined by measuring the flux trapped in it after a known current has been 
stored. The flux is measured by the time integral of the voltage decay across 
the leads S^ a and as Hj is turned on. This method has been described in 
Ref. 9. There is also a heat-switch, Hg, between Ljj and that will be 
turned on in all operations involving current changes. This heat-switch 
serves to protect the Josephson junction of the SQUID ring against any surge 
current induced in the input coil 1^. The currents Ij and I 2 in the two 
sensing loops can be fine tuned by letting the stored flux leak very slowly 
across an ultra low resistance path after a section of the superconducting 
path has been turned normal. If the heat-switch is turned on, the current 
I 2 in will decay across R 2 g with a decay time t = (L 2a + L 2 ^)/R 2 a - 

3 x 10^ s, which is extremely long. This current will rise or decay even more 
slowly if a second current, slightly greater or smaller than I 2 , is being sent 
across the leads S 2a and S 2 ^ while H 2 g is kept on. An additional current ij 
can be added in the loop formed by a parallel combination L^ a / /L^ and Lj- (see 
Fig. 3 of Paper I) by turning both and H 2 g on and then off while a current 
is sent across Sj and either Sj a or Sj^. The decay of 1^ and I 2 during the 
time required to store i^ is negligible because of the long decay time t. The 
directions of these persistent currents are as defined in Fig. 3 of Paper I. 

The superconducting circuit is shielded inside Nb. The outer surfaces of 
the Nb holders for sensing coils are shielded with additional Nb covers to 
form junction boxes. Circuit components such as superconducting joints, a 
superconducting transformer, shielded heat-switches and ultra low-value 
resistors are mounted in these junction boxes. 

The superconducting joints and transformers are similar to those used in 
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References 9 and 10. A pair of Nb lead wires to be joined are twisted 
together, wrapped around a brass screw, sandwiched between two Nb sheets, and 
clamped together under a brass washer. They are clamped hard by tightening 
the brass screw into a piece of insulated brass block. The Nb wires are 
chemically cleaned with concentrated HNO 3 before they are joined. Flatness 
and parallelism of the brass block surface relative to the screw head surface 
is found helpful for making very tight joints without breaking the thin Nb 
wires. A fresh section of the lead wire is used to make the joint. Joints 
made this way consistently have critical currents in excess of 3 A through 
several thermal cycles [11]. The transformer makes use of a Nb sheet to guide 
the magnetic field through the primary and secondary coils. 

Shielded heat-switches are made with a small cylindrical Nb cup which 
shields the superconducting circuit from the noise of the heater resistor. 

The heater element is a small 500 Q Allen-Bradley resistor which is glued with 
GE 7031 varnish inside the Nb cup. Shielding from the resistor is continued 
with superconducting lead-tin (PbSn) tubing which starts from inside the Nb 
cup to enclose the lead wires of the heater. The sections of Nb wire to be 
affected by the heat-switch are wound non-inductively on the outside surface 
of the Nb cup and are glued with the GE varnish. The cup is then glued on a 
piece of fiberglass which is held rigidly inside the junction box. The Nb 
wires from the heat-switches are heat-sunk to the Nb mass of the junction box. 

The shielding of the heat-switches provides significant improvement in 
preventing electrical pickup and rf interference at the input of the SQUID 
amplifier. The shielding has a drawback, though, in possessing a higher heat 
capacity which prolongs the response time of the heat-switch. By making the 
heat-switch compact, we managed to obtain an on/off response time of about 0.2 
s with 4 mW supplied to the heater. 
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The ultra low-resistance heat-switch makes use of the radial conduction 
across a very thin cylindrical wail of copper. we simply solder Lwo lengths 
of copper-clad Nb wires together with PbSn solder. The copper coating gives 
extremely low resistance between the superconducting PbSn on its outside wall 
and the Nb wire inside. The resistance between the two Nb wires comes from 
two such resistors connected in series. This resistance is given by the 
formula : 


R 5 ' 


Pca <4 - 2 K) 


it A 


An 


Cu 


Cu 

r Nb 


(4) 


We used Nb wire of diameter 2r Nb = 0.128 mm with copper coating which gives an 
overall diameter of 2rQj = 0.166 mm. The resistivity of copper at 4.2 K is 

_Q 

PCu “ 10 ® m * Therefore, when a wire of length = 10 mm is used, we 

obtain = 10 Q. This resistance has been verified in a decay time 
measurement of magnetic flux. The other two ends of the copper-clad Nb wires 
have their copper layer removed with concentrated HNO 3 so that superconducting 
joints can be made there. In the rest of the superconducting circuit, Nb 
wires without copper coating are used to avoid ac losses due to eddy currents 
[9]. 

In the original circuit used for tests, an additional SQUID was coupled 
to the superconducting circuit to measure the common acceleration simul- 
taneously. Two identical superconducting transformers were inserted into the 
paths that connect the two sensing loops to (A and B in Fig. 3), with one 
on each accelerometer side. The primaries of these transformers completed the 
gradiometer circuit. Their secondaries were connected in series with the in- 
put coil of the second SQUID. The current induced in this loop was propor- 
tional to the sum of the signal currents from the two accelerometers, per- 
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mitting a simultaneous readout of the common acceleration in the "gradiometer 
mode". Although such a readout is useful for error compensation and active 
vibration isolation [12], we have omitted, for simplicity, this additional 
circuit in the later tests in favor of the basic circuit shown in Fig. 3. 
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III. CONSTRUCTION AND PERFORMANCE OF TEST APPARATUS 


A. Mechanical and Thermal Design of Experimental Space 

The superconducting gravity gradiometer is operated in a very quiet cryo- 
genic vacuum space with a stable temperature. A low boil-off liquid helium 
(He) dewar was designed and constructed with great care. Radiation shields 
with superinsulation are used in preference to a liquid nitrogen (N 2 ) shield, 
because boiling liquid N 2 is a much noisier source of vibration than boiling 
liquid He. Pressure pads are sandwiched near the bottom of the dewar between 
the fiberglass inner tube and the sturdy outer A1 wall. The relatively weak 
inner tube is then damped against swinging motions. The inside dimensions of 
the dewar are 0.378 m in diameter and 1.778 m in height. 

A low heat-leak insert has been constructed to complement the dewar. An 
A1 vacuum can is supported by three equally spaced fiberglass tubes and one 
center fiberglass tube which permits access to the vacuum space from the dewar 
top. The vacuum space is 0.33 m in diameter by 0.61 m in height. The can 
remains totally submerged for one week with 60 liters of liquid He. A long 
hold-time of He is needed for an uninterrupted operation of sensitive 
measurements. The low boil-off rate is also important for lower boiling noise 
from the He and for lower thermal drift of the experiment. 

For most of the experiment, the gradiometer assembly was suspended inside 
vacuum by means of a fiberglass rod from the room temperature end of the dewar 
insert. All the lead wires are heat-sunk in liquid He before they enter the 
vacuum space through specially made feedthroughs. In this setup, long term 
temperature variations of the gradiometer assembly over many weeks average 
about 10 mK. In a time scale of several hours, the temperature variations are 
within the 2 mK resolution of a germanium (Ge) resistance thermometer mounted 
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on the gradiometer . The heat capacity of the gradioraeter and the weak thermal 
coupling to the He bath form a low-pass filter for the temperature fluctua- 
tions between the gradiometer and the He bath. The response time is estimated 
to be on the order of minutes. 

B. Electromagnetic Shielding 

External magnetic field fluctuations and electrical pickup can produce a 
high noise level in comparison with the very weak gravitational force. Rf 
interference will also deteriorate the noise performance of the rf SQUID and 
can unlock the SQUID feedback loop. The first stage of magnetic field isola- 
tion is provided by a double wall mu-metal shield which surrounds the cryo- 
stat. Two additional superconducting stages of magnetic shielding are pro- 
vided by a 1 pm lead (Pb) film evaporated on the inner wall of the vacuum can, 
and by the Nb body of the gradiometer itself. Also, the A1 walls of the dewar 
and the vacuum can aid in shielding high frequency electromagnetic radiation. 

The success of the electromagnetic shielding is evidenced by the achieve- 
ment of the optimum SQUID amplifier noise when the SQUID input is coupled to 
the gradioraeter at full sensitivity. The magnetic shielding has also been 
tested when the gradioraeter is balanced. A magnet producting 600 Gauss being 
flipped just outside the dewar produced a response less than an equivalent 
flux change of 2*10”^ or 4x10”^ Wb at the SQUID input. 

Each current lead for storing a supercurrent has an inductive choke 
inside the vacuum can. The choke is a high-impedance path to attenuate the 
electrical noise that the wire can pick up outside the Pb shield and 
especially outside the dewar. The other leads pass through two stages of rf 
filters, one inside the Pb shield and the other outside it. The outside one 
is an RC low-pass filter with the manganin wire down the dewar providing the 
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resistance and the inside one is a capacitor in parallel with the load. 

C. Vibration Isolation 

Vibration proved to be the most troublesome source of noise. The test of 
the gradiometer has therefore been conducted in a deep underground laboratory, 
which is secluded from the traffic of people. The most critical noise data 
have been obtained in the evenings to minimize the mechanical and electrical 
disturbance of the experiment. In spite of these precautions, the gradiometer 
required careful mechanical isolation from the environment. 

The vibration isolation system consists of a three stage low-pass filter 
schematically shown in Fig. 4. The gradiometer assembly is suspended by a 
fiberglass rod and a length of latex rubber tubing in vacuum which form the 
first stage of vibration isolation. The vacuum is extended above the top of 
the dewar insert by means of a thick-wall plexiglass tube. The rubber tubing 
inside this plexiglass tube hangs from a connecting rod which is attached to a 
massive brass plate that seals with the plexiglass at the top. The height of 
the connecting rod can be adjusted from outside the vacuum so that the gradio- 
meter can be lowered to the bottom of the vacuum can to short out the filter. 
When the gradiometer is suspended by the rubber tubing, the vertical resonance 
frequency is 0.9 Hz. 

The fiberglass and rubber also give wideband isolation from tilt noise . 
The gradiometer suspended by a long fiberglass rod constitutes a pendulum. A 
natural property of a pendulum is its tendency to align itself along the 
direction of gravity. The gradiometer is therefore decoupled from a crucial 
error source : the ground tilt . which would cause an error signal due to 

modulation of Earth's gravity for a hard-mounted gradiometer. 

Another important property of the pendulum is that the mass at the end of 
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Fig. 4. 



the rod is under a two-dimensional "free fall" in the horizontal plane. An 
oCCciciOmcLci luOuuteu Oii uue pendulum bub therefore cannot experience any 
horizontal acceleration relative to the platform to the first order. A 


gravity gradiometer can, however, couple to the angular velocity and the 
angular acceleration which result from the pendulum response to the horizontal 
acceleration (see Section V and Appendix A for detailed discussions of the 
pendulum action). The net effect of the pendulum action for the horizontal 
vibration is that of a low-frequency filter, in series with a wideband 
attenuator. 

Additional resonant modes of the suspension structure can deteriorate 
vibration isolation at the resonance peaks. In spite of these limitations, we 
have chosen the pendulum suspension for the gradiometer because of its isola- 
tion characteristic for low frequency horizontal acceleration and its simpli- 
city in construction. The pendulum frequency (swinging mode) is approximately 
0.3 Hz and the rubber tubing gives a torsional mode at 0.08 Hz. In addition, 
the massive gradiometer suspended by the long and therefore bendable fiber- 
glass rod gave a troublesome high-Q peak at 2.5 Hz. This peak was the result 
of a rocking mode of the gradiometer assembly around a horizontal axis. 

The massive brass top of the plexiglass column is suspended with a set of 
rubber tubings. The bottom of the plexiglass column is connected to the top 
of the dewar insert by a flexible vacuum bellows. The rubber tubing and this 
bellows form the second stage of the filter with a vertical resonance fre- 
quency of 2.8 Hz. This filter can be shorted out by connecting clamps across 
the bellows . 

The third stage of the filter is formed by three sets of rubber tubings 
that lift the dewar off the floor. The dewar is suspended from a framework by 
these rubber tubings and has a vertical resonance frequency of 0.5 Hz and a 
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swinging frequency of 0.3 Hz. The upper joints of the rubber tubings in this 
suspension are adjustable so that the dewar can be levelled. Again, this 
filter can be eliminated by lowering the dewar to the cart. 

The advantages of using latex tubing as a passive filter are its high 
elasticity and low Q. A latex tubing of initial length can be extended to 
a final length A f which is a few times A t . The cross-sectional area will 
decrease from A^ to Af = A^A^/Af) in order to keep the volume constant. We 
assume that the modulus of elasticity is approximately constant even with such 
large extension. One can easily show then that the vertical resonance 
frequency is 


f 


J_ f g E a 1/2 Al/2 

2it '■A f - A i ) V A f j 


(5) 


where A f is the new length determined by gravity g E acting on the load. The 
factor corrects for a change in the cross-sectional area as the 

rubber elongates. By having large extension Af - Af, a low resonance fre- 
quency is obtained. The low Q provides self damping of the own resonance of 
the filter. Our latex tubing suspending the gradiometer in vacuum has Q 55 10. 

The disadvantages of rubber are in relaxation and drift which worsen as 
the latex extends. A compromise is to keep Af £ 2Af. We used Af = 1.6 Af for 
the rubber inside vacuum. Over a period of three months under load, it 
stretched by an additional 3 cm. 
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IV. GRADIOMETER PARAMETER MEASUREMENTS AND CIRCUIT OPTIMIZATION 


A. Determination of Geometric Asymmetry 

The main asymmetries between the two component accelerometers in a 
gradiometer come from the spacings of the sensing coils and mechanical spring 
constants. A high degree of matching for these parameters was not attempted 
because the mismatch can be compensated for by supercurrents. As can be shown 
from Section III of Paper I, the low frequency common mode balance condition 
in the small stored current limit with only two currents Ij and I 2 reduces to 


L 2b //L 2a t0 2M 
L lb //L la W IM 



( 6 ) 


The currents and I 2 can therefore be adjusted to compensate for the 
mechanical asymmetry, and the ratio thus obtained measures the original 

asymmetry before it was balanced* We shall determine the ratios on both sides 
of Eq. (6) to check against the theory. This check will also indicate how 
good a "symmetric component-accelerometer approximation” is. Such an approxi- 
mation shall often be used later for simplicity. 

The parameters measured in a particular cool-down are summarized in Table 
1. The uncoupled mechanical frequencies of the proof masses were measured at 
4.2 K with the gradiometer assembly firmly attached to an effectively infinite 
mass platform. This platform is simply the dewar and the ground. We have 

2 2 

measured aJ 2M / ^ a) lM = an< * the mas s ratio is estimated to be close to 

unity. Therefore, the ratio of the measured mechanical frequencies reflects 
the mismatch in the spring constants for the two proof masses. 

The experimental values of the various inductances shown in the circuit 
diagram of Fig. 3 can be measured only when the system becomes superconduc- 
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Parameter 


Symbol 


Value 


Mass of each proof mass 
Mass of gradiometer assembly 
Mass ratio 

Mechanical proof mass frequency 
Mechanical spring constant ratio 
Gradiometer Baseline 
Pancake coil area 
Pancake coil wire density 
Inductance per spacing 
Coil-to-proof-mass spacing 
Inductance of one coil 

Measured coil spacing 
asymmetries 

Sensing loop asymmetry 
Geometric asymmetry factor 


“k 

0.40 kg 

M 

7.7 kg 

ro^/M 

0.052 

f 

19 Hz 

kM 

W 2M /W 1M 

0.92 

a 

0.16 m 

a l 

2019 mm 2 

n L 

11 mm - ^ 

2 

\ = Pq n L A L 

0.306 H nf 

A 

0.14 mm 

L 

L 

42 pH 

Hb - ha 

0.15 

L lb + L la 

L 2b " L 2a 

0.34 

L 2b + za 

L 2b 11 ^2a 

0.91 

L lb " L la 

L 2b 11 L 2a W 2M 

0.83 

L, v // L, 2 

lb la » m 


Table 1. Experimental parameters of the gravity gradiometer 
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ting. The method of stored flux measurement described in Section IIB was em- 
ployed to measure various loop inductances s The measured series i nHurtanr^s 
L lb + L la and + I- 2 a agree with the inductances calculated from Eq. (2). 

In order to measure the individual inductances or the parallel inductance 
L| ca //L lcb , we used an acceleration-response method to measure the ratios 
Lib/l«ia and This method does not require knowledge of the behavior 

of Lq. A very small stored current is used so that negligible shifts are pro- 
duced in the equilibrium position of the proof mass and in the proof mass 
resonance frequency. When a current is stored on only one side, say on side 
one, the accelerometer sensitivity is proportional to the effective current: 


r . I + hb ~ hi i 

1 X 1 L, . + L, X 1 


J lb 


la 


L lb T la L la I lb 

L,. + L, 
lb la 


(7) 


as can be seen from Eqs. (34) and (46a) of Paper I. The accelerometer sensi- 
tivity is first measured with Ij a => I t , Ij ^ = 0; and is then remeasured with 
lla “ *lb “ *t* rat ^° these two sensitivities is Iq^/Lia Eq» 

(7). The asymmetry factors (L lb - L la ) / ( L lb + L la^ and ( L 2b “ ^a^^b + 
L 2 a ) are then found to be 0.15 and 0.34, respectively, and the ratio of 
L2b//l‘2 a to Iq b //Iq a is 0.91. The geometric asymmetry factor can now be 
calculated as 


L 2b^ L 2a M 2M 
L lb //L la W 1M 


0.83 . 


The supercurrent ratio, required to compensate for this asymmetry is 

found to be 0.84, so that the two measurements are in agreement within the 
errors. 
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B 


Accelerometer Node and Gradiometer Node 


In the gradlometer mode, the supercur rents are stored such that only 
differential accelerations at the two component accelerometers will produce 
signals. A common acceleration is "current-dif ferenced" out before detection. 
If the relative polarity between the supercurrents at the two component 
accelerometers are reversed, the differential acceleration signals will then 
be current-differenced out. The gradiometer is, thus, activated in a “common 
accelerometer mode” . Being an identical electromechanical system, the 
accelerometer mode has the same characteristics toward common accelerations as 
the gradlometer mode has toward differential accelerations. 

In the experiment, the application of a pure acceleration is much easier 
than that of an acceleration gradient. Interchangeability of the roles 
between common and differential accelerations therefore provides us with an 
alternate method to evaluate and calibrate the gradiometer. This method is 
highly sensitive because a known acceleration can be applied with a much 
higher amplitude than a gravity gradient. 

In each mode of circuit operation , the two coupled component accelero- 
meters have two degrees of freedom and will have two normal modes of motion . 
The two proof masses move with the same phase in one mode at a frequency a> + 
and with opposite phases in the other mode at a frequency These normal 

mode frequencies and Q f s can be measured experimentally and provide informa- 
tion on the gradiometer parameters. We will first show how they are related 
to the gradiometer parameters. 

The equations of motion and their constraints have been given in Eqs. 

(41) and (42) of Paper I. Strictly speaking, the normal mode frequencies w + 
and 03_ are obtained as eigenf requencies of these equations. This can be done 
by first eliminating the current coordinates i ( j(w) and i c (to). For free 
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oscillations, i.e., when g d (io) = g c (u)) “ 0, the dynamical equations become 


2 x r (L 0 +L lp )m l 


2 , L 0 m 2 


(-0) + W. )[ + l] X. + (-0J + 03, )-= 

1M k 2 ,2 1 2M a. 2 Ill’ 2 


2 , l 0 "l 


2 


(“<*> + U 1 M )~2 X 1 + + W 2 M ^ 2 — 2 + ^ x 2 

1M A. 2 !! II 1 2M A. 2 !'? 2 


where the frequency dependence of x^ and X 2 have been omitted for simplicity. 
The resulting eigenf requencies are given in Ref. 10. We shall obtain approxi- 
mate solutions in simpler forms to illustrate the physics. 

The two component accelerometers in our gradiometer are partially matched 
geometrically. We use the symmetric component-accelerometer approximation: 


m l = m 2 > (9a) 

U 1M = * (9b) 

L ls = Ms * <10a) 

Mp - Mp • < 10b > 


In this case, a wideband balance could be achieved without introducing i^ and 
i 2 » If we now take Ij = -I 2 in the gradiometer mode, the eigenf requencies are 
just given by 

u+ - wi c - w 2c , (11a) 

w_ - w ld * u 2d , (lib) 
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to the leading terra in the above approximation. In the accelerometer mode in 
which Ij ■ I 2 , the frequencies are reversed: 

u + = ^ld * u 2d » (12a) 

o)_ * w lc - w 2c • (12b) 

The frequencies o)^ c and u)^ are effective frequencies for the 
"complementary'* current (i c ) and "additive" current (ij) modes, as is 
evidenced from their definitions in Eq. (43) of Paper I, and are theoretical 
parameters for describing the operation of the gradiometer. The normal mode 
frequencies of motion, u) + and a)_, on the other hand, are physically measurable 
quantities that provide us with approximate values of o>^ c and 

When the gradiometer assembly is suspended by a weak spring system of a 
low resonance frequency Q, the gradiometer platform has a finite mass M and 
will couple with the pair of masses. The three frequencies, Q, eo + , and u)_, 
will all be shifted, and the equations of motion are modified. If we denote 
the platform position by X and measure the proof mass coordinates relative to 
it, the kinetic energy of the system is 

T = j [MX 2 + m^X + x^ 2 + m 2 (X + x 2 > 2 ] , (13) 

1 2 2 

and the potential energy V has an additional terra y MQ X . The linearized 
equations of motion, Eq. (37) of Paper I, are then modified into 

-co 2 [(M + mj + m 2 )X + + m 2 x 2 ] + MQ^X = 0 , (14a) 
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(14b) 


2 2 
A if 


_ 9 1 

-arm^X + xj) + mi ( oj 1m + — •] Xj + AIJ ijCto) » 0 , 

1 Is 


. 2.2 

2 . K 2 


-arn^X + x 2 ) + m 2( w 2M + — j- — •) x 2 - AI£ i 2 («>) = 0 . (14c) 

2 2s 


These equations are subject to the constraints given by Eq. (38) of Paper I. 

We again make the approximation of Eqs. (9) and (10), along with i^ = i 2 
= 0. In the gradiometer mode (1^ = -I 2 ), the eigenf requency of the "out-of- 
phase mode" (xj = -x 2 ) of motion is unchanged: 



(15a) 


whereas, the "in-phase mode” (xj = x 2 ) eigenf requency is shifted to 

2 2m v 2 

■ <* + -r> \c ■ (15b) 

2 2 

Here we have made an approximation based on » Q . In the accelerometer 
mode (Ij = 12 )* we find 


b) 


A- 


0), 


kc * 


(0 


A+ 


2m k 2 

' (1 + TT* “kd 


2 2 2 
The parameters and u)^ are functions of 


2 

V 


2 

\d 


2 *** 


\M + « k <L. * 2L„) > 


(16a) 


(16b) 


(17a) 


(17b) 
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as obtained from Eqs. (44) of Paper I. 

The shifted frequencies, fg_, f q + , f^_ and f^ + , measured as functions of 

2 

I k are presented in Fig. 5. The intercepts to 1^ ■ 0 are given by 

f G- (0) ■ 4- <0) ■ f W > <I8a) 

4 <°> - 4 + (0) - (1 *^r> 4 • (18b) 


Taking the ratio of the two intercepts enables a measurement of the mass 
ratio. We obtain m^/M ~ 0,052 ± 0.006, which is in agreement with the value 
given in Table 1. 

2 2 2 

The slopes of the graphs for co A _, w A+ , Wq 

A 2 2m u A 2 A 2 2m v a 2 

-A (1+ ^ A n+ _ fey. A. _ 

W M W 2 V ’ \ ( V 2 V ’ M "VV ' 

respectively. From these, we obtain an inductance ratio: 


U> 


G+ 


versus 


* 


are 


— = 1.58 ± 0.14 . (20) 

V 

This result indicates that the output impedance has been chosen sufficiently 
large. Using the value of A = 0.306 H m“* calculated from Eq. (15) of Paper 
I, we find 1^ = 23 pH, which is consistent with inductance values in Table 
1. Also, we obtain Lq = 18 pH. Consistency in the measurements indicate that 
the symmetric component accelerometer approximation is sufficiently good for 
diagnostic purposes. 

The output inductance Iq is the equivalent inductance at the transformer 
primary when the secondary is loaded by the SQUID input of approximately 2 pH. 
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SQUARE OF NORMAL MODE 
FREQUENCIES, f 2 (IOO Hz 2 ) 


20 



SQUARE OF STORED CURRENT, lj[(A 2 ) 


Fig. 5. Normal mode frequencies of the gradiometer as functions of the 
stored current. 


-Ill- 



The transformer has a turns ratio of 69:19 with calculated inductances of 42 

|aH and 3.2 pH, respectively, to a 20 % accuracy, can then be calculated 

from Eq. (75) of Paper I, and would have the above measured value of 18 pH if 

2 

we take the transformer coupling factor to be = 0.93. However, the power 
transfer function of the transformer would then only be H*’ ** = 0.51 ± 0.10, 
as one can see from Eq. (76) of Paper I. The transformer has not been 
optimized in the past assembly, but this is not a problem until one reaches 
the amplifier-limited sensitivity. The resulting value of the effective 
electrical energy coupling to the SQUID is t) = 0.31. 

C. Circuit Optimization and Transfer Functions 

Equation (53) of Paper I shows that the gradioraeter transfer function, 

d 2 

Hgi(co), of an equivalent accelerometer is proportional to I ' /(jJ q in the low 

signal- frequency limit. In the symmetric component-accelerometer mode, co 0 is 

identified with 

2 

In Fig. 6, we plot a measured curve of lljJ/f* versus |l^| to indicate 

the functional dependence of the dif ferential-acceleration-to-response-current 

transfer function on the stored current 1 1^ | . This plot shows that the 

gradiometer sensitivity first increases as | 1^ | is increased from 0, and then 

slows down due to an increase in f Q caused by an increasing |l^|. The gradio- 

2 2 

meter is expected to reach its maximum sensitivity when = 2f^ at I I ~ 2 

A, However, this maximum is very smooth so that a smaller current can be used 
without a substantial reduction in sensitivity. 

If the gradiometer sensitivity is SQUID noise limited, optimizing the 
above transfer function yields a maximum signal-to-noise ratio for the gradio- 
meter. However, if the sensitivity is limited by the Brownian motion noise, 

1/2 2 

one would have to maximize the function Q_ 1 1-^ I / f 0 instead, where Q_ is the 
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Fig. 6. Sensitivity scaling of the gradiometer as a function of the 
stored current. 
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quality factor of the out-of-phase mode: fg * f Q . 

In Section III of Paper I, we used the response current in Lq to define 
transfer functions in harmony with the model of a reduced equivalent accelero- 
meter. Here we prefer using the response current at the SQUID input a 
quantity which can be directly calibrated using the SQUID transfer function. 

The in-phase mode of motion has the following parameters: 


m + = 2mj < . , 


x 2 + X 1 


x + = 2 


( 21 ) 


( 22 ) 


In the calibration experiment, a common acceleration is provided by applying a 
sinusoidal displacement X(oo) to the gradiometer assembly. In this case, Eqs. 
(14) can be solved to obtain 


(-to 2 + u^) x + (w) = U) 2 XO.o) = g (u). 


(23) 


This acceleration produces a response current i + (w) which is proportional to 
x 2 + X 1 = 2x + . Th e displacement-to-current transfer function is 




i + (w) 2i + 

x + (w) (x 2 + Xj) + 


(24) 


At a low frequency w « to , the acceleration-to-current transfer function is 


A , , i + (w) 2 

H gi (w) _ / X 2 


( 25 ) 


g c (w) “o (X 2 + X l ) + 


The parameters for the out-of-phase motion are 
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( 26 ) 


m. 



J 


X_ - X 2 - Xj . 


( 27 ) 


Notice that m_ is the effective mass m of the gradiometer defined in Section 
HIE of Paper I. In the gradiometer mode, a differential acceleration g d (to) 
produces a differential proof mass displacement x_ with 


(-or + w ) x_(u) = g d (w). 


( 28 ) 


The response current i_(w) is proportional to just x 2 - Xj = x -* The dis “ 
placeraent-to-current transfer function is therefore 




i_(w) i_ 

X_(U) “ (Xj - Xj)_ * 


( 29 ) 


Th e gravity-gradient-to-current transfer function for w « t») 0 is 


H^ i (o>) 


i_(<o) ^ i_ 

r(w) (x 0 - x. ) 

O l 1 - 


( 30 ) 


where A is the baseline and r(o>) = g d (o j)/A. 

If the same magnitudes of stored currents are used in the two modes of 
circuit operation, the following equality is obeyed by the two equivalent 
circuits : 


i 



(x 2 + x l } + 


( 31 ) 


The transfer functions in the two circuit modes are therefore related by 
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- 1 *&<“> 


(32) 


This simple relationship follows because Eqs. (25) and (30) contain the same 
resonance frequency u) 0 , as is expected from Eqs. (lib) and (12a). 

The function nhw is measured by applying a sinusoidal acceleration 
g c (u) = to^X(o)) to the gradiometer assembly at a) below 1 Hz. The response 
current is measured with the SQUID amplifier which has a current-to-voltage 
transfer function of 0.20 V pA - *. The acceleration is applied to the gradio- 
meter in the vibration isolation system described in Section IIIC. An 
electromagnetic shaker is coupled with a weak latex tubing from one side to 
the other of a vibration isolation stage (see Fig. 4) to provide a vertical 
acceleration to the gradiometer. The weak tubing acts as an attenuator, 
permitting a good signal-to-noise ratio without degrading the vibration 
isolation. The platform acceleration produced by the shaker has been 
calibrated by measuring its displacement and later confirmed with a commercial 
accelerometer. 

Figure 7 shows the calibration of the accelerometer mode measured in 
terms of the current response i + at the SQUID input as a function of the 
applied acceleration g c . Very small stored currents are used because of the 
extremely high sensitivity of the accelerometer. From the graph, we obtain a 
transfer function of (1*55 ± 0,05) mA per m s for the accelerometer with 
only Ij^ = 1 mA stored on each side. The transfer function should be linear 
with the stored currents for low current values that produce negligible shifts 
of proof mass frequencies. This linearity is confirmed by Fig. 8, which is a 
plot of the measured acceleration-to-current transfer function as a function 
of I k . 

At high currents, the calibration curve in Fig. 7, which takes into 
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Fig. 8. Sensitivity scaling in the accelerometer mode as a function of the 
stored current. 
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account frequency shifts, is used to extrapolate our transfer function 


measurements c We then obtain 




i + (u>) 

g (y) “ (l* 6 ± 0.1) A (m s -2 ) 1 


(25*) 


for a near optimum current value of I k = 1.5 A. Converting this value to the 
transfer function in the gradiometer mode, we have 


H^Cw) - f — -2 = (0 * 13 * °- 01 ) nA E" 1 . (32') 

lms 

The SQUID amplifier noise, when loaded with the fully charged gradio- 
meter, was measured to be S I 1 ^ 2 (f) = 8.9 pA Hz -1 / 2 , which corresponds to E A (f ) 
= 7.9x10 2 ^ J Hz ^ 2 for = 2 pH, at a high frequency (> 60 Hz) where the 

proof mpsfj rH flnl acpmpnf noi ao becomes notr1_£cH kl o TV»^e o<a l Q , T al T . T/ ^,,i A 

correspond to an amplifier-limited gravity gradient noise of 0.069 E Hz -1 / 2 . 

This value is in close agreement with the theoretical sensitivity of 0.068 E 
- 1/2 

Hz which is obtained by directly substituting the measured SQUID noise and 
the measured values of q = 0.31 and 8 = 0.53 into Eq. (79) of Paper I. In our 
tests, the gradiometer noise at low frequencies was limited by seismic noise 
at 0.3 ~ 0.7 E Hz ^ 2 level. Therefore, we would not lose sensitivity by 
using a somewhat smaller gain. Stored currents of 0.55 A were used in the 
tests reported in Section VI. The transfer function deduced from the 
accelerometer mode calibration gives 62 pA E" 1 at this current level. In 
Section VIA, we will check this value with the direct calibration of the 
gradiometer which employs gravity gradient signals. 
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V. INVESTIGATION OF INSTRUMENT NOISE AND ERRORS 


A. Environmental Noise and Common Mode Balance 

Vibrations were found to be the most serious environmental noise. Our 

early tests with the gradiometer rigidly mounted inside a liquid N 2 shielded 

—8 

He dewar have encountered a vibration noise energy of 10 J at the proof mass 
resonance frequency of 19 Hz. Such a noise level corresponds to an equivalent 
Brownian motion at a noise temperature of 7 x 10^ K, which is 10 times the 4 
K Brownian motion in vibration amplitude! 

When the gradiometer is balanced against common accelerations, its 
response to the vibration noise is greatly reduced. In the experiment, 
however, the gradiometer can be balanced only if we can use its acceleration 
response as a guide. Sufficient vibration isolation to achieve a noise level 
of the unbalanced gradiometer well within the amplifier dynamic range is 
therefore essential. A minimum of 60 dB isolation is required at the proof 
mass resonance frequencies in order to be able to operate at the full 
sensitivity of the accelerometers. 

If vibrations have a white acceleration spectrum, the low frequency noise 
will be Q times less than that at the proof mass resonances; thus, this noise 
will not overload the amplifier. Therefore, we chose to start with vibration 
isolation at higher frequencies where the SQUID amplifier was being over- 
loaded. This isolation would enable us to balance the accelerations and 
obtain a working gradiometer. 

If the proof mass acceleration noise amplitude can be reduced by four 
orders of magnitude, one will be able to observe the actual Brownian motion at 
the resonances in a gradiometer with a 60 dB balance. Although such an obser- 
vation may not improve the acceleration noise at the low frequency end, it 
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could prove to be very informative. This observation would determine the 
■Rrnwnfan motion level of the instrument, and show whether the observed 
acceleration noise comes from the environment or from the intrinsic noise of 
the instrument. 

Although an electronic “cold damping" technique [8,10] could have been 
applied to selectively filter the vibration noise at the proof mass 
resonances, there were other modes in the environment which require wideband 
isolation. Therefore, we chose to develop a passive vibration isolation 
system with which the gradiometer can be tested in its simplest form. 
Satisfactory vibration isolation of the gradiometer was an essential step in 
leading to a successful demonstration of its performance and the operation of 
the device as a detector in the inverse square law experiment [7,5]. Various 
combinations of low-pass vibration filters were tried to isolate the noise 
coming from the ground and the dewar. Filters were improved in steps with the 
guidance of vibration data obtained in each step. 

With the three-stage vibration isolation system described in Section 
IIIC, the noise at the proof mass resonance decreased to 2x10^ K, corres- 
ponding to a displacement of 5 x 10 A m. The achieved improvement at 20 Hz 
was therefore better than 100 dB. At frequencies above 60 Hz, the spectrum 
indicated the SQUID-limited noise level. 

We now turn to the wideband common mode balance method described in 
Section HID of Paper I. In this method, calibrated sinusoidal acceleration 
signals are provided in the vertical direction by means of an electromagnetic 
shaker* Two out of the three stored currents 1^, I 2 and ±2 are iteratively 
adjusted for balance at 0.9 Hz (« u> Q / 2it ) and 70 Hz (» goq/2h;) in turn. This 
is a painstaking procedure since a balance point is searched for in a two- 
dimensional parameter space and sufficient time should be allowed after each 
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adjustment of currents for the mechanical disturbance of the system to damp 
down. 

Figure 9a shows the 0 - 50 Hz noise spectrum of the gradiometer balanced 

in the wide band to 1 part in 10 3 . The transfer functions are obtained from 

the accelerometer mode calibrations using the relationships discussed in 

Section IV. This procedure gives [H^. (w)] * = /i = 1.29 m s -2 A - * . The 

gi o - - 

gradiometer is less sensitive to external common accelerations now. Shown in 

the spectrum is the differential mode of the masses at f Q = ( 0 o / 2 it = 24.8 Hz of 
2 -? 

an rras amplitude u) q x_ = 2.75 nm s . The noise temperature for this mode is 
given by 


1 “kv 22 . . _ 

(x— ) to x - 4.5 K 
k B 2 ° " 


(33) 


Thus, this spectrum represents the observation of the Brownian motion of the 
proof masses in the differential mode at the liquid He temperature. All other 
peaks are attenuated by about 10 3 , which is the degree of balance obtained. 
This data confirms that the observed noise is almost entirely environmental 
and is not intrinsic to the gradiometer. 

Figure 9b shows the 0 - 10 Hz spectrum. The vertical scale of the plot 
has been calibrated by using the accelerometer mode, and confirmed by a direct 
calibration with a known gravity gradient signal which will be explained in 
Section VI. The noise in the 5 - 10 Hz band reaches below 1 E Hz~^ 2 but the 
background seems to be close to 0.3 E Ilz -1 ^ 2 . Excessive vibration noise has 
degraded the spectrum at several peaks. The noise in the 1 - 4 Hz range is 
heavily contaminated by the rocking mode of the gradiometer (f r ) at 2.5 Hz and 
other modes associated with the vibration isolation system. The 0 - 1 Hz 
spectrum shown in Fig. 9c exhibits the swinging modes (f sl , f g2 ) around 0.3 
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Hz, the vertical dewar suspension mode (f vl ) at 0.5 Hz, the vertical gradio- 
raeter rubber suspension mode (f v 2) at 0.9 Hz, and a torsional mode (f t ) at 
0.08 Hz. 

Our low-pass filters introduced numerous peaks at their own resonances. 
These filters were needed initially to prevent the SQUID from overloading at 
the higher frequency peaks. After the common mode balance, the gradiometer 
became intrinsically insensitive to the environmental noise. Therefore, in 
the gradiometer mode, the filters could be reduced to a minimum. We found a 
single-stage filter provided by the rubber tubing inside the vacuum was 
sufficient for a balanced gradiometer and eliminated many undesirable peaks 
from the signal bandwidth. 

B. Gradiometer Noise 

Figure 10 shows the noise spectra of the gradiometer obtained with the 
single-stage filter. The upper trace is the noise spectrum measured with the 
gradiometer in the accelerometer mode. The lower trace is the gradiometer 
output after the common mode balance to 3 parts in 10^ has been achieved 
toward an applied acceleration at 0.18 Hz. The SQUID noise limit of 0.069 E 
Hz with the associated 1/f noise is plotted by a dotted line. The 
vertical scale shows both the acceleration and gradient noise calibrations, 
which are related by Eq. (32). The ratio of the two spectra at any frequency, 
therefore, represents the actual reduction of noise outputs at that particular 
frequency as a result of the balance. 

First, comments are due for the measured acceleration noise spectrum (top 
trace). The general structure of this spectrum, which contains a broad peak 
around 0.25 Hz and a dip around 0.1 Hz, and the noise level of approximately 
2.5 x 10 ms Hz ' at 0.5 Hz are obtained consistently from run to run 


- 124 - 


GRADIENT ACCELERATION 

-l/o -2 -1/2 

NOISE (E Hz ) NOISE (/xm s Hz 



FREQUENCY (Hz) 


Fig. 10, Low frequency noise spectra of the gradiometer on a single-stage 
vibration filter. 
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after waiting sufficiently long for resonances in the gradioraeter platform to 
damp down. This noise data has been reproduced by using conventional 
accelerometers mounted at the room temperature end of the pendulum suspension 
and on the laboratory floor, with the exception of the dip at 0.1 Hz for which 
the room temperature accelerometers did not have enough resolution. Since the 
data represents the residual acceleration noise of a deep underground 
laboratory, the upper curve in Fig. 10 must reflect geophysically significant 
seismic activity, except for the broad peak near 1 Hz which is due to the low- 
Q vertical mode of the rubber suspension. It is noteworthy that, while 
acceleration noise spectra measured at the top end of the pendulum in the 
vertical and horizontal directions looked similar, the noise output of the 
superconducting accelerometer correlated closely (> 90%) with the vertical 
acceleration. This correlation indicates that the horizontal acceleration is 
indeed attenuated by the pendulum action, as has been predicted by the 
analysis shown in Appendix A. Since the gradiometer is pointing at an 
umbrella angle © u - tan -1 /T « 54.7°, the noise spectrum of Fig. 10 should 
then be raised by a factor of sec 0 U = /3 to obtain the proper calibration of 
the vertical acceleration noise at College Park, Maryland. 

The lower trace of Fig. 10 represents a typical residual noise spectrum 
of a balanced gradiometer, with the exception of the 0.1 Hz peak of the 
torsional mode which disappears into the background noise in the quietest 
situations. Unlike the torsional mode, however, the swinging mode at 0.32 Hz 
is always present with its peak at around 100 E Hz - */ 2 level. Also 
persistently present is the low frequency excess noise below 0.1 Hz. The 
noise level is approximately 0.7 E Hz -1 / 2 between 0.5 and 0.9 Hz and reaches a 
minimum value of 0.3 E Hz */ 2 in a narrow window around 0.15 Hz. Comparison 
of the two traces reveals that the noise was reduced only by a factor of 3 x 


- 126 - 


q 

10 , which represents a discrepancy of a factor of 10 from the value expected 
from the common mode balance. Another noticeable feature is that, although 
the general structure of the background spectrum of the balanced gradioraeter, 
with the resonance peaks removed, appears similar to the shape of the 
acceleration noise, the two spectra are not well correlated. Therefore, the 
residual noise of the gradiometer must be coming from a source other than 
vertical acceleration. 

It is desirable to determine the Browian motion noise level of the 
instrument. The power spectral density of the thermal noise of the gradio- 
raeter is given in Eq. (87) of Paper I: 


S r>B (f) 


8 

m^A 


2 T(f) • 


(34) 


it is difficult to measure the effective damping time x(f) at f « f 
directly. If one assumes the "white acceleration noise model", t(f) in Eq. 
(34) can be replaced by the measured quantity, t(f 0 ). The measured Q values 
are of the order of 5 x 10^ and are probably pressure limited. This implies 
x(f) = 3.2 x 10 2 s * and Sp^(f) - 0.01 E Hz”*/ 2 for our gradioraeter. If, on 
the other hand, the effective Q is assumed to be independent of frequency, the 
Brownian motion noise would become frequency dependent: 


sj/g(f) = 0.01 E Hz" 1/2 (|-l l/2 . (35) 

* o 

In either case, the thermal noise would be negligible compared to the SQUID 
noise level of 0.069 E Hz”*/ 2 . 

We therefore conclude that the observed noise spectrum is extrinsic to 
the instrument and is due to error mechanisms which couple environmental noise 
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into the gradiometer output. Error sources responsible for the observed 
residual noise of the gradiometer are identified in the following section with 
the aid of Appendices. 

C. Instrument Errors 

The dynamics of a platform suspended as a pendulum is analyzed in 
Appendix A. The pendulum decouples the gradiometer from angular motions of 
the external world. The horizontal acceleration transmitted to the gradio- 
meter is almost exactly cancelled by the modulated Earth's gravity which is 
caused by the induced tilt of the instrument. It is found, however, that the 
resulting angular velocity and angular acceleration can couple to the gradio- 
meter through centrifugal acceleration and through the misalignment of the 
effective gradiometer sensitive axis with respect to the baseline. Since the 
vertical component of linear acceleration is balanced out to a sufficient 
degree, the horizontal acceleration can produce a limiting instrument error by 
driving a secondary angular motion. 

Experimentally, the noise of the unbalanced gradiometer correlates with 
the external acceleration in the vertical direction. This correlation and the 
absence of a pronounced resonance peak at the pendulum frequency fp = 0.32 Hz 
in the upper trace of Fig. 10 prove that the horizontal vibration isolation by 
the pendulum action really works for an accelerometer . A close examination of 
the noise spectrum reveals a barely visible resonance peak at fp above the 
background noise. This result renders a remarkable proof for the theory 
developed in Appendix A. According to Eqs. (A13), the horizontal acceleration 
sensed by an accelerometer under the pendulum suspension should be peaked at 
Its resonance with a net gain of unity. Since the measured frequency spectra 
of the horizontal and vertical accelerations agree within 10 dB, the supercon- 
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ducting accelerometer pointing at an umbrella angle should measure a hori- 
zontal acceleration with a maximum amplitude of approximately 3 times that of 
the vertical acceleration at f p . Therefore, the horizontal acceleration peak 
sticks out barely above the vertical acceleration background. 

The analysis in Appendix A also leads to a prediction that, for a 
gradiometer. the overall effect of the pendulum action is that of a wideband 
attenuator for a horizontal acceleration, with its attenuation improving as 

A 

the frequency approaches dc. In particular, when co > Q E /|6n A |, the dominant 

+5L 

coupling mechanism for the horizontal acceleration is the induced angular 
acceleration, as given by Eq. (24). Below this frequency, the centrifugal 
acceleration error dominates. The cross-over frequency is 10 - ^ Hz when 

A A 

|5n A | = 10 • The frequency response of the angular acceleration, |h (co) I , 

+£ a 

has been plotted in Fig. A2. One can see a strong resemblance between this 

curve and the gradiometer output, the bottom trace of in Fig. 10, above 0.2 
Hz. In fact, the agreement improves when the spectrum of the driving back- 
ground acceleration is taken into account. In Fig. 11, we plot typical ob- 
served spectra for the vertical (top dotted line) and horizontal (top dashed 
line) accelerations ^(w) and ^(w) as well as the gradiometer noise output 
(solid line). The residual vertical acceleration expected from the vertical 
common mode balance of 90 dB is plotted by the bottom dotted line. The bottom 
dashed line is a fit obtained by multiplying the spectrum for X^(w) with its 
transfer function through the angular acceleration. The misalignment error 

A . 

required for this fit is |5n A | = 7 x 10 , which is in close agreement with 

the value determined from the observed centrifugal acceleration: |5n A | * 6 x 

-4 

10 (see Appendix B). Notice that the gradiometer noise above 0.2 Hz is 
completely accounted for by this single error mechanism: the horizontal 

acceleration coupled through the induced angular acceleration . 
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Fig. 11. Various error contributions to the gradiometer noise spectrum. 
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The value of 6n A (the misalignment of the average sensitive axis of the 
+i 

gradiometer with respect to the baseline vector 1) found above is rather 
large in view of the fact that the parallelism of the mating surfaces of the 
mechanical components of the gradiometer was specified to 5 x 10 -5 . This 

A 

mechanical precision implies that the experimental value of 6n_ (the relative 
misalignment between the sensitive axes of the component accelerometers) is of 

W.A A 

the order of 10 • A similar value could have been obtained also for 6n A if 

the center lines of the two accelerometers had been located within 1CT 2 mra 

from each other, which is certainly feasible. However, without realizing the 

seriousness of the induced angular acceleration problem, we had left the 

centering of the accelerometers to crude positioning by the mounting screws. 

This must have produced a mismatch of the center lines by approximately 0.1 

^ * 

ram, causing the observed misalignment |6n A | = 6 x 10 . 

We now turn our attention to the excess noise below 0.1 Hz, seen in Fig. 
10. During the experiment, the SQUID electronics was shielded against rf 
■f ni-Arf<»r»rw'<» fn avoid electronics drift. Although instrument drift can be 
caused by relaxation of trapped magnetic flux [13], the observed low frequency 
noise did not decrease significantly when the ambient magnetic field was 
reduced by a factor of 10^ by adding mu-metal shields. With these electrical 
sources of low frequency noise eliminated, we suspect the temperature drift 
and various nonlinear effects as the most likely sources for the observed 
excess noise. As has been discussed in the Appendix of Paper I, the supercon- 
ducting gradiometer can be rather sensitive to temperature noise and 
nonlinearities in the system can down-convert the wideband vibration noise to 
a gradient noise near dc. These low frequency error mechanisms are analyzed 
in detail in Appendix B of this paper. Owing to the limited resolution of the 
Ge thermometer used to monitor the temperature of the gradiometer assembly, 
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the low frequency noise spectrum due to temperature fluctuation is computed 
from an independent measurement of the temperature noise of a liquid He bath 
of a different experimental setup [8] with the aid of a model for heat 
conductivity between the He and the gradiometer. The resulting spectrum is 
found to be 


[S r>T (f)] 1/2 = 1.0 x 10 -2 E Hz“ 1/2 ( -• ■ 1 - f - Hz l 2 . (36) 

The scale-factor nonlinearity is computed from the second order circuit 
equations for component accelerometers, generalized from Eqs. (29) and (33) of 
Paper I. The low frequency noise expected from this is of comparable 
magnitude as the temperature effect at 0.1 Hz, but has a slightly different 
frequency dependence: 

tS I\cc (f) ] 1/2 1 x 10 ~ 2 E Hz_1/2 ) 3/2 • (37) 

The effect of the centrifugal acceleration is found to be about two orders of 
magnitude smaller than this value. 

The temperature noise of Eq. (36) is plotted in Fig. 11 by a dash-dot-dot 
line. The actual low frequency noise, denoted by X(u>), is about 30 times 
larger than this calculation. Since the temperature drift of the He bath de- 
pends on the exact geometry and thermodynamics of the particular cryostat in- 
volved, and on account of additional uncertainties in the thermal model of the 
instrument, the discrepancy by a factor of 30 does not exclude a possibility 
that the observed low frequency noise is due to temperature fluctuations. 
However, because a comparable level of noise is also predicted from the non- 
linearity of the instrument, we cannot be conclusive, without further experi- 
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mental Investigations, as to the origin of the observed excess noise below 0,1 


Hz. In future experiments, the tempera Lure of the He bach will be regulated, 
and the gradiometer will be linearized by means of a force rebalance feedback 
with a hope to eliminate the low frequency excess noise. 


Another interesting error mechanism discussed in Appendix B is the 
second-harmonic generation by the centrifugal acceleration. This effect has 
been observed for large resonance peaks when the system is disturbed 
mechanically (Fig. B2). During the quiet operation of the instrument, how- 
ever, no second harmonics were observed above the background noise spectrum of 
the gradiometer. 
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VI. GRAVITY GRADIENT MEASUREMENT AND SENSITIVITY CALIBRATION 

We shall now describe the gravity gradient calibration with the observa- 
tion of signals in two types of experiments. A precise calibration of the 
Instrument can be obtained by using a periodic gravity gradient generator and 
averaging the signal over a number of periods. For such an ac experiment, 
however, any electrical pickup or mechanical cross coupling at the signal fre- 
quency must be eliminated. A ”dc" experiment, In which a compact object Is 
brought near the gradiometer and removed, fares better in this regard. 

Although the use of signal averaging to improve signal the signal-to-noise 
ratio is less applicable here, any cross coupling problem will appear as a 
transient which can be distinguished from the dc shift in the gradiometer 
output caused by the gravity signal. We have performed both experiments 
successfully, as well as an inverse square law test in which the gradiometer 
outputs in three orthogonal orientations of the sensitive axis are summed. 

A. AC Gravity Gradient Experiment 

The ac calibration of the gradiometer was obtained with a large (M = 1.6 
x 10 3 kg) Pb pendulum. The Pb sphere was located approximately in the same 
horizontal plane with the gradiometer at an average distance r of 2.3 m away. 
The gradiometer sensitive axis was at an angle of cot ^ ^2 with respect to the 
horizontal plane and was in the plane of swing of the pendulum. The gradio- 
meter thus measured a gravity gradient component given by 

Tj! = -^|[i + 3 cos (2 cot -1 /2)] = - • (38) 

r r 

When the pendulum swung by Ar(w), it generated a change in the gravity 
gradient by 


- 134 - 


( 39 ) 


Ar u (co) - Ar(w). 

r 

The SQUID amplifier output was measured with a spectrum analyzer. 

In order to assure that the observed signal was truly gravitational, 
various precautions were taken. When all three vibration isolation stages 
were used, excess signals at various harmonics of the pendulum frequency were 
seen. Correlation studies indicated that a strong second harmonic was 
generated by an air current driven by the swinging pendulum. This effect was 
eliminated by shielding the entire apparatus with a wooden enclosure and by 
using only one rubber filter which was inside the vacuum (see Fig. 4). The 
possibility of direct mechanical coupling between the pendulum and the gradio- 
meter has also been investigated. With the gradiometer charged in the 
accelerometer mode and the pendulum swinging at its full amplitude ±30°, we 

Hi H nnl* Gi»#» anv ar>r«a1 arof-l r\rt 1 J.. 1 -C _ i , 

- w jys.nuui.um awuvc LUC gUJUUU 

noise. The problem of such mechanical coupling has been alleviated by several 
factors. The pendulum was suspended from the 2 ra thick concrete ceiling; 
whereas the dewar sat on a firm floor, 12 ra below the ground level, in the 
basement of the building. Furthermore, the gradiometer was balanced against 
common accelerations, and the pendulum suspension used for the gradiometer 
gave isolation of tilt at all frequencies. In order to avoid electrical 
pickup at the signal frequency, a pneumatic driving mechanism was developed 
for the Pb pendulum and a natural resonant oscillation was used for most 
tests. 

The rras voltage measured at the pendulum frequency of 0.228 Hz is plotted 
in Fig. 12 against the calculated rras gravity gradient generated by the 
pendulum. This plot gives a calibration of 58 pA E - ^ and agrees with the 
accelerometer mode calibration, reported in Section IVC, to within 7 %. This 


- 135 - 



APPLIED GRAVITY GRADIENT, T(E) 


Fig. 12. AC calibration of the gradiometer. 
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agreement Is within our errors. The distance and orientation measurements in 
either calibration were not intended for high precision. 

B. DC Gravity Gradient Experiment 

The noise spectrum of Fig. 10 indicates that the noise is lower below the 
0.3 Hz swinging mode of the platform. In order to perform a dc experiment, we 
put in active low-pass filters which filtered off the peaks above 0.16 Hz from 
the amplifier output. The response function of the filters to a step dc 
voltage change is shown in Fig. 13a. The characteristic rise time and ripples 
of the filters will be seen when a step function, gravity gradient signal is 
applied. 

We then moved a Pb brick in and out manually near the dewar at 10 ~ 15 s 
intervals. The brick weighed 12 kg and was brought to 0.17 m from the gradio- 
meter. The gravity gradient generated was about GM/r J = 40 E, with a geome- 
trical factor of the order of unity which depends on the exact source-detector 
orientation. 

In Fig. 13b, the Pb brick was moved in, out and then in again, roughly in 
the horizontal plane of the gradiometer. The times and directions of the 
movement of the Pb brick are indicated by arrows. The line of sight of the 
source was approximately at an angle of cot - * /2 from the gradiometer axis. 

The expected signal is then 

3 (1 " 71 = 37 E * (40) 

r r 

o o 

The AVr term is the finite baseline term of the gradiometer. Here the Pb 
brick has been approximated as a point mass. Its actual geometry and size 
should be taken into account for a precise calibration. The measured peak-to- 
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Fig* 13. DC gravity measurements (dc level is arbitrary). 
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peak dc current change is 2.3 nA. This current corresponds to (38 ± 2) E, in 
good agreement with Eq. (40). if the accelerometer mode calibration of Section 
IVC is used. 

For Fig. 13c, the brick was moved in the same way but along a direction 
roughly perpendicular to the gradiometer axis. The expected signal, if the Pb 
brick were a point mass and were located accurately, is given by 


GM 2 

3 U ~ 8 2 


r 


r 


) * -34 E. 


(41) 


The value measured from Fig. 13c is (-25 ± 2) E. An error of 10 % in r and in 
angular orientation could produce this gradient error of 30 %. The brick has 
a long rectangular shape so that a substantial amount of its mass was at other 
orientations and at larger separations relative to the gradiometer. 

The experiments have been qualitatively repeated at a few other orienta- 
tions relative to the gradiometer axis. The results agree with what we expect 
from the tensor property of the single-axis gradiometer. The validity of a 
gravity signal has therefore been substantiated. 


C. Inverse Square law Experiment 

The most impressive demonstration that the gravity gradiometer output is 
a genuine gravity gradient signal comes from the inverse square law data [5, 
7] in which the angular pattern of the gravity gradient tensor is checked. 

For this experiment, the sensitive axis of the single-axis gradiometer was 
turned into three orthogonal directions by rotating the entire experimental 
cryostat on a turn table incrementally by 120° around the vertical axis while 
the gradiometer was suspended in the umbrella angle. The resulting three 
outputs and their sum are plotted in Fig. 14 after signal-averaging over 500 
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pendulum periods. The peak-to-peak amplitudes of the three outputs agree with 
the theoretical formula; 


Ar il = cos (42) 

where r 35 2.3 m is the average distance between the Pb sphere and the gradio- 
raeter, Ar = 0.84 m is the peak-to-peak amplitude of swing, and = 120 °(i - 
1) - 3.7° is the angle that the projection of the i-th sensitive axis onto the 
horizontal plane makes with the line of sight between the detector and the 
source. 

The sura of the righthand side of Eq. (42) over i * 1 to 3 vanishes as ex- 
pected from the Poisson equation for Newtonian gravity, Eq. (12) of Paper I: 

V a n ... r\ / / \ 

/, 41L ^ - V. 

i 

The sum of the three orthogonal measurements, the bottom trace of Fig. 14, 
verifies this condition within experimental uncertainty: 

l Ar u = (+0.15 ± 0.23) x 10 -9 s“ l 2 . (44) 

1 

The signal and noise in this result have been analyzed carefully and published 
as a null test of the gravitational inverse square law [5]. Within the 
inverse square law, the mutual cancellation of the signals in three arbitrary 
orthogonal directions is strong proof that the instrument is really measuring 
genuine gravity gradient signals since it will be highly unlikely that other 
error signals also satisfy the condition (43). 
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VII. SUMMARY AND CONCLUSION 


A superconducting gravity gradiometer has been constructed and evaluated 
very carefully. Both the signal transfer and noise characteristics of the de- 
vice agree closely with the theoretical model. Practical error sources have 
been identified experimentally as well as analytically. Although the instru- 
ment has only been tested on a passive platform without any feedback or error 
compensation, it has shown a performance level of 0.3 ~ 0.7 E which 

surpasses the value achieved in a sophisticated room temperature gravity 
gradiometer [14] by three orders of magnitude in power. Further reduction of 
noise should be possible by improving the alignment of component accelero- 
meters or by actively stabilizing the gradiometer platform against horizontal 
accelerations. The temperature of the He bath could also be controlled to 
reduce the low frequency noise of the instrument. 

The instrument has been used to carry out a preliminary null test of 
Newton’s inverse square law of gravitation. The successful operation of the 
superconducting gradiometer in a rather simple setting demonstrates that the 
macroscopic quantum mechanical phenomenon exhibited in superconductivity can 
be taken advantage of to meet many practical challenges posed by sensitive 
gravity measurements. The careful modelling of the instrument and its experi- 
mental confirmation establishes the feasibility of constructing a more ad- 
vanced superconducting gravity gradiometer and carrying out precision gravity 
experiments. 

Guided by the theoretical and experimental work reported in Papers I and 
II, we are in the process of constructing a three-axis superconducting gravity 
gradiometer of an advanced design [6, 15]. The new gradiometer incorporates 
additional technical innovations such as a "superconducting negative spring" 
and a "three-dimensional residual common mode balance" as well as the cold 
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damping and force rebalance feedbacks to achieve an Instrument noise level of 

a _ __ -i _ . . . . . .... 

1U ' K Hz iioj . Such a sensitive and stable gradlometer will find 

Important applications in Inertial navigation and geophysical survey as well 
as In fundamental physics experiments. A series of precision inverse square 
law tests as well as tests of general relativistic effects are planned with 
this instrument both In the terrestrial laboratory and in space [17, 18]. The 
work presented here lays the ground work for further development of supercon- 
ducting technology for gravity and acceleration measurements. 
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APPENDIX A: VIBRATION ISOLATION BY PENDULUM ACTION 


Isolation of horizontal vibrations by pendulum action has been analyzed 
in Ref. 7 by treating the angular motion of the pendulum in time domain. Here 
we treat the problem in frequency domain. For a platform suspended in a 
stationary laboratory, the angular amplitude of the pendulum is of the order 
of 10 6 rad at its resonance frequency. Effects that arise from second order 
contributions in angle will therefore be negligible in general. 

1. Vibration Isolation of an Accelerometer 

We first consider a single accelerometer mounted on a platform suspended 
by a symmetric pendulum of effective length Ip, the distance between the pivot 
point and the "center of percussion" of the platform (see Fig. Al(a)). The 
effective length of the pendulum is defined from ujp, the degenerate (angular) 
resonance frequency of the swinging mode, by 

*p H Sg/Up * (A1) 

Let the center of mass of the accelerometer proof mass be located at a 
vertical distance 1 Q below the pivot point with its sensitive axis defined by 

A 

n. We further assume a "perfectly rigid” pendulum in which the pendulum mass 
and the suspension rod form a perfectly rigid body and the pivot is perfectly 
rigid for linear motions while it is perfectly soft for angular motions. This 
condition could be approached in practice by designing the pendulum with the 
resonance frequencies of unwanted modes high compared to the swinging mode 
frequency u) p . 

Under these assumptions, the accelerometer platform is completely iso- 


- 145 - 





lated from the angular motions of the celling from which the pendulum is sus- 
pended. This isolation is independent of frequency although, in practice, 
finite rigidity for the torsional mode tends to make the pendulum a low-pass 
filter for the torsional jitter of the external world. No attenuation is ex- 
pected for vertical vibrations. Since linear accelerations can be balanced 
out in one direction in a gradiometer, our interest is in the extent to which 
horizontal accelerations a^(t) can be rejected by means of the pendulum 
action. 

Let us denote the horizontal displacements of the pivot point and the 
center of percussion of the platform by ^(t) and r h (t), respectively. Then, 
the equation of motion of the pendulum can be written as 





x *(r, - $. ) 
p h h 


2,+ 
M n (r 
P h 


‘ *h> 


0. 


(A2) 


Taking a Fourier transform of this equation, one finds 


r h (w) 2 2 


u* - + jut 1 

J l -K m - 


(i) - (0 + jut 

P P 


(A3) 


The horizontal acceleration of the platform is therefore 


2 -1 

u + jut 

a h (u>) = -w 2 r h (w) - ~2 — E__ 


U - 0) + jut 

p p 


r r\ _y 

where A^(u) = - u^ R^(u) represents the external acceleration. 

The resulting angular displacement of the platform in the Earth 
coordinate system is given by 


(A4) 


$ (w) “ — [r h (u) - Vu)] x z , 
P P 


(A5) 
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to the first order in "§ . Substitution of Eq. (A3) into Eq. (A5) leads to 


*,<•> ' 2 2 


-1 


0) - U) + jWT P 

P P 


ZjifL 


r V"> x 2 • 


(A6a) 


jp ^(w) x z * (A6b) 


B (w) ” jw t (») = -z , 

P P 4 ^ » j ^ 

v v 0) - w + iojt p 

p p 

2 

a (u>) = -w 2 $ (u) = — ^ -y ^(w) x z . (A6c) 

P P to - 03 + jtOT p 

p p 

Notice here that & (w) represents the angular velocity of the platform with 
respect to Earth. 

The total specific force acting on the accelerometer proof mass is 


g'(t) = n(t) • [g E - a ho ( t ) ] , 


(A7) 


where a^ Q (t) is the horizontal acceleration of the platform at the proof mass, 
which is given by a weighted average of £^(t) and a^(t) : 


%, (t > - <■ -r > V*> + r\< t) 

P P 


(AS) 


The Fourier transform of Eq. (A7) is 


V A y A 

g(o)) = 0 (to) x n • g E - n • a hQ (o)), 


(A9) 


to the first order in , where the vector identity (All) of Paper I has been 

A 

utilized to obtain n(w) and 


So <u) 


W 2 + jWT ^-(l-!l/i.)a3 2 

2 2 

2 2 -1 
w - w + iwx 
p p 


V u >' 


(A10) 
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■> * 

Substitution of * - gg z into Eq. (A9) gives the form: 

A A 

g(u>) - n • [g E e p (w) x z - a ho (co)]. (All) 

Now, we substitute Eqs. (A6a) and (A10) to obtain 

(1 - * /* ) <o 2 - jcot" 1 A 

g(w) - 2 2 31 n * M»). (A12) 

<o - to + i(OT 
p p 

where the main contributions to the two terras in Eq. (All) have cancelled each 
other by Eq. (Al). Namely, the horizontal acceleration of the proof mass is 
almost exactly balanced out by the modulation of Earth's gravity which results 
from the tilt of the platform. 

According to Eq. (A12), residual errors in the horizontal vibration iso- 

latlnn arl RPR from fwn CAM mac • ~ ~ c 

- - - — r w - U V k w v/i. Ill CL ij O 

(A 0 * ^p) anc * damping (Qp = tOpTp < °°). When = Ap is chosen, the attenua- 
tion obtained by the pendulum action goes as 

(A13a) 
(A13b) 
(A13c) 

At resonance, the platform is actually driven to an amplitude Qp times larger, 
as can be seen from Eq. (A10), and the pendulum action balances the horizontal 
acceleration to Qp^ so that the net gain will be unity. Away from resonance, 
the horizontal acceleration is attenuated by times a frequency factor 


(u)/o) p ) Q p X ,w « (Op, 


g((p) 

In • *„(“>! 


1 , CO = 0) p , 


(Wp/to) Q p ^ , <0 » (Op. 
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which helps the isolation further. Thus, we have proved that the pendulum 
suspension indeed provides wideband isolation of an accelerometer from 
horizontal accelerations as well as tilts except in the immediate vicinity of 
its resonance frequency. 

We have assumed symmetry of the pendulum frequency in two horizontal 
directions for simplicity. Equation (A12) shows that this assumption is not 
necessary for the operation of the pendulum action. Since the accelerometer 

A 

is not sensitive to a horizontal acceleration perpendicular to n, the fore- 
going analysis will be valid as long as A Q is chosen to be equal to the effec- 

A 

tive length of the pendulum A p in the vertical plane which includes n. 

2. Vibration Isolation of a Gravity Gradiometer 

We now turn to a gradiometer mounted on a platform under pendulum suspen- 
sion (Fig. Al(b)). In general, the two component accelerometers have slightly 

A A 

different sensitive axes, n^ and n£, and are located at different distances, 
and A 2 , from the pivot point. The distance between the pivot point and the 
center of percussion, which is near the midpoint between the two proof masses, 
is denoted by A Q . Since a very small differential acceleration signal is 
sought for in general, errors coupled through misalignments of the sensitive 
axes will have to be examined carefully. The result of a general analysis of 
the misalignment errors is given by Eq. (A6b) of Paper I, which can be written 
in the form: 


6I\ „(t) = j 5n 
n+A 


[g E - Ut)] 


(6n . 

+n 


+ 6 h • ff E - &(t) £(t)] 


- 5n A x n 

+i 


a(t), (A14) 
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where A Is the baseline of the gradioraeter, and 6n_, Sn A , 5n A , 5 A are the 

+A +n 

misalignment and misoripnfflHnn i/ertny-c aof h*r v<i a ( _ fQ/. \ t> 

— *•'''• «'t w » \ ' — / \ -^ -r / ui i. a[ic J. 

I* Th e gravity gradient (T^) term can be ignored in comparison with the 
gravity (g^) term, and the correction term to centrifugal acceleration can be 
dropped since the full centrifugal acceleration is regarded as a separate 
error source in Paper II, 

Upon Fourier transformation, Eq. (A14) becomes, to the first order in % % 

P 

= X ^ n - # ^E^p^) x z - a^Cto)] + 6n A x n • a^(o)),(A15) 


where appropriate quantities for the pendulum suspension have been substituted 

into the linear and angular accelerations. Here a , (oj) is the horizontal 

ho 

acceleration of the midpoint between the two accelerometers. In deriving Eq. 

(A15). use has been made of the rel^M onslvtn ho t“ uonn t"h o flr<a f- I & ** #- 

1 “ — - ^ 

two accelerometers, a^Cto), and a^ o (a)) : 



*ho<*> 


2k = • a *“ "X 

*o sin [0^ - (-l) k <J,] 


(A16) 


to obtain a perfect balance of the accelerations along the direction A: 


t 


ft h2 ((o) 


a hl (u)] = a h2 (o)) sinO^ -<{>)- a hl (t*)) sin(9 A + <(>) = 0. (A17) 


The angles 9^ and <|) are defined in Fig, Al(b). In addition, we have assumed 
A Q » A to approximate 

1 A -► + * + 

2 6n _ • k* h2 (w) + a^Cco)] “ Sn _ * * (A18) 
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Substitution of Eqs. (A6) and (A10) into Eq. (A15) results in 


6 T\ A (w) 
n+A 


a - * o /y » 2 - 

2 2 -1 

w - o> + jorr 
P P 
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(A19) 


The first term in Eq. (A19) is again minimized by choosing Z Q = £ p ; i.e., 
by locating the center of the gradiometer at the center of percussion of the 
platform. In this case, however, the horizontal acceleration couples in 
through the second term which represents the effect of the induced angular 
acceleration. In order to balance out in the two terms, Z a should 

rather be chosen to satisfy 




■) 6n_ + ( 
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z) 6n A 
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(6n A • z) n * 0. 
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(A20) 


Since there are three independent components in this equation, it is impos- 
sible in general to satisfy this equation by adjusting a single parameter Z Q 

A A 

even if 6n and 6n A are predetermined. Therefore, in order to obtain a wide- 
+A 

band rejection of the horizontal accelerations in a gradiometer by means of 

A 

the pendulum suspension, a sufficiently small value of 6n A is required. 

+SL 

An additional error term is generated through the centrifugal accelera- 
tion of the gradiometer proof masses induced by ^(w). This error is given by 
Eq. (A10) of Paper I: 


6r (w) = -2 [(n • <$_,) {n • ^ («)} -&_•&. (<o)], (A21) 

£ E p E p 

where is the constant angular velocity of Earth. Substitution of Eq. (A6b) 

hi 
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into this equation yields 


6 r_(w) = ~2 2^ If Y~ * ^g) n “ * ^ h (w) x z * ( A22 ) 

C to - <o + jure p 1 

P P 


Combination of Eqs. (A19) and (A22) gives the total residual error arising 
from l^(co) : 


6r h (u>) 
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P P 



-r— X (<o) x z .(A23) 
x n 
P 


Equation (A23) shows that the pendulum works as an attenuator with an attenua- 

A 

tion factor proportional to (A/Jt_) |5n „ | and (H/Si ) (Qg/o)) , respectively, for 

V +i V 

each error term, plus a high-pass filter with the cut-off frequency to p . 
Therefore, for a gradiometer, the pendulum action attenuates the horizontal 
accelerations at all frequencies, with the isolation becoming perfe ct as the 
signal frequency approaches dc . This is usually sufficient because the signal 
bandwidth satisfies to < <o p in most applications of a gradiometer. It is truly 
remarkable that the pendulum action accomplishes a nearly perfect , passive 
vibration isolation in two dimensions near dc because vibration isolation 
below 1 Hz is very difficult to achieve by means of conventional spring-mass 
isolators. 

/\ 

In the frequency range 0 ) > Qg/|Sn * | , the centrifugal acceleration terra 

+1 

can be ignored in comparison with the angular acceleration term so that 


J A A . A 

6r h (w) = h (to) -r— fin A X n • A. (to) x z , 

p +1 


(A24) 


a 


where 
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(A25) 


h (w) = 


w 

2 2 -1 
0) - (i) + j(j)T 

P P 


Is the normalized transfer function for coupling of IC^Oo) through induced 
angular acceleration. The absolute value of this function is plotted in Fig. 
A2 (solid line) for parameter values of <d p /2it = 0.32 Hz and Q p = (J p T p = 300. 
Also plotted in the figure (dotted line) for comparison is the absolute value 
of the transfer function for direct transmission of 1 ^(oj) to the pendulum 
mass : 


h t (w) = 


2 . , -1 
tO + j0)T 

E 1_2 

2 2 -i » 

a) - to + iorc 
P P 


(A26) 


which is read from Eq. (A4). Thus, the pendulum itself responds to the 
horizontal acceleration as a low-pass filter with a frequency characteristic 
given by Eq. (A26). Equation (A6a) shows that the induced angular displace- 
ment has an almost identical frequency response which cancels the direct 
transmission term to approximately Q p ^* What is left over after this cancel- 
lation is the frequency characteristic of a high-pass filter given by Eq. 
(A25) . 
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TRANSFER FUNCTION 



FREQUENCY (Hz) 


Fig. A2. Transfer functions of the pendulum suspension for horizontal 
accelerations . 
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APPENDIX B: LOW FREQUENCY NOISE OF THE GRAD IOMETER 


Gravity gradient signals of interest are usually at very low frequencies 
(1(T 4 ~ 1 Hz). Excess low frequency noise (red noise), which raises the 
overall instrument noise above the Brownian motion noise, is therefore of 
great concern in constructing a sensitive gravity gradiometer. Red noise in 
the gradiometer output can be generated by thermal and mechanical drift of the 
gradiometer, nonlinearities in the instrument, centrifugal acceleration of the 
platform, the 1/f noise and drift of the detecting electronics, etc. One way 
of overcoming certain types of such low frequency noise is an up-conversion of 
the gravity gradient signal by rotating the instrument followed by a 
heterodyne detection of the modulated signal [19]. The rotation, however, 
increases the centrifugal acceleration and introduces additional dynamically 
induced noise sources. In this Appendix, we analyze how various error 
mechanisms contribute to the excess low frequency noise for a non-rotating 
gravity gradiometer. 

1. Temperature Drift 

Temperature coefficients of the superconducting gravity gradiometer have 
been derived in Appendix of Paper I. It has been found that the dominant 
effect comes from the functional dependence of the penetration depth \(T) on 
temperature T. The good thermal conductivity through the gradiometer body 
makes the differential temperature fluctuation Tj(o)) negligible compared to 
the common temperature fluctuation T c (to). The temperature-induced error can 
therefore be written as 


^£>d,T^ ,A 0 - h^ c (co) T c (to), 


(Bl) 
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where 


h Tc («) = f«2^ W 2\ “ ^ " «X X " (i + O > (B2) 


T 2 

and a^, co^ are defined by Eqs. (A42) and (A45) of Paper I, respectively, and 
C < 1 is the modulation efficiency of the penetration depth for the pancake 
coils. 

In the experiment, i^ = 0 was chosen so that 


T „ Hia Hcb 

“ k = hcs 


(B3) 


^kX r a) kM* 


(B4) 


O..U « t- 4 *... *. -f — ~ -f — - T? /DON *-1 — ~ J 1 - r *-1 - 4 . 4 — 

L»UUOUXUUti.Ug JL lib V uv i * V XJ*~ / U. lit. lUCOOUL^U VOXUCO VJ. pci X ailll^ UV, L O XXO CV^U XU 


Table I and 


dX 

dT 


2X(0) 

T 

o 


- . .3 

VA/ V 

[i - (t/t q ) 4 ] 3/2 


9.7 x 1<T 10 m/K 


(B5) 


for Nb at T = 4.2 K, one obtains 


e^(oj) = j h Tc (co) = - 2.9 x 10 4 E K -1 , « « WfeM » (B6) 


for the present gradiometer. 

Thus far the temperature-induced error has been analyzed in the frequency 
domain. Although this description is complete, it is often convenient to 
measure and compensate for the temperature effect in real time. Notice that 
hij. c (w) becomes independent of u when w << Therefore, for a slow drift 
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with t » 


one expects a linear relationship : 


“ ^Tc ^c^^ * 


(B7) 


where 


h.pc - hi^g (u - 0) • (B 8 ) 

The observed noise T c (t) in time domain can be related to the power 
spectral density S^(f) by a standard formula for the digital FFT (fast Fourier 
transform) [ 20 ] : 


S T (f) = < | \ X J^ n T c (t) e * 2%ft dt | 2 >, (B9) 

where x is the sample length. For a linear drift, 

T c (t) * at, (BIO) 


this equation yields 


S T«) -72 <J J- • f ■ (BU) 

4x f 

Notice that the result depends on the choice of x. Combination of Eqs. (Bl) 
and (B9) leads to 


Sr >T (f) 



<a 2 > 


4it 


(B12) 
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Thus, a linear temperature drift causes a 1/f 2 power spectrum at the gradio- 
■ueter output. 

In the experiment, the gradioraeter was thermally isolated from the He 
bath through vacuum. The large heat capacity of the gradiometer body and the 
low conductivity through the electrical leads and the residual He gas formed a 
low-pass RC filter for heat flow with a transfer function: 


h F ( f ) 


i 

1 + J2xfT F ’ 


(B13) 


where x F is the filter time constant. In high vacuum, x F exeeded 10^ s. In 

fact , we had to maintain a pressure level of approximately 10 - ^ mm Hg to bring 

the gradioraeter assembly into thermal equilibrium with the He bath in a 

reasonable time after turning on heat-switches on the gradiometer. A typical 

response time uaeu was i F = 300 s. A Ge t’nerraoraecer was mounted on the 
gradiometer assembly to measure T c (t). However, the thermometer did not have 
enough resolution to track the stable temperature of the gradioraeter. All we 
could determine directly is a rough estimate of the experimental upper limit, 

2 x 10"* K hr -1 , for dT c /dt so that 


a rras i 5 ‘ 6 x 10 ~ 8 K s_1 * 


(B14) 


Substituting x = 125 s, the sampling time used for the spectrum analyzer, and 
Eqs. (B6) and (B14) into Eq. (B12), we obtain 

tS r, T (f > J 1/2 £ 2 * 9 x 10 ~ 2 E Hz“ 1/2 (B15) 

In a separate experiment [8 ] , the temperature of the He bath was measured 
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by dipping a "superconducting penetration depth thermometer" of the type 

similar to the one described in Ref. 21 into a storage dewar. The Fourier 

analysis of this data is shown in Fig. Bl. This data indicates that, in the 

—2 —1 

frequency range of interest between 10 * Hz and 10 1 Hz, the spectrum can be 
described by 

[S T (f > ] B^h = 1A X 10 ~ 4 K HZ_1/2 ( - - f 1 • (B16) 


Since the exact behavior of the temperature noise should depend on the 
characteristics of the dewar and the insert as well as the barometric pressure 
fluctuation in the laboratory, this data should be interpreted as an order-of - 
magnitude estimate for the temperature noise in our experiment. Combining Eq. 
(B16) with Eqs. (B6) and (B13), one finds 


l8 r,i«> ]1/2! 2E: E ! ls t«)lMh - ‘•° * 


10" 2 E Hz 


- 1/2 


0.1 Hz \2 


V. (B17) 


This estimate falls within the approximate limit given by Eq. (B15). 


2 m General Consideration of Nonlinear Effects 

Since nonlinearity is a higher order effect, the most important error 
coupling mechanisms will involve the second order nonlinearity . Thus we 
restrict our discussion to a nonlinear behavior of the form: 


6g d , N (t) = f kj^(t “ t 1 ) g^Ct 1 ) dt 1 • (B18) 

The driving acceleration g(t) could be either g c (t) or g d (t). In this 
section, we wish to obtain a relationship between the Fourier components g(oo) 
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and 6g<i fl(w), and investigate specifically how the wideband noise g(w) 
contributes to a low frequency divergence in Sg^ jj(co) through the squaring 
process. 

The Fourier transform of Eq. (B18) is 

= h N (co) G(uj), (B19) 

where G(u>) is the Fourier transform of 

G(t) = g 2 (t). (B20) 

In experiments, one usually measures one-sided power spectral densities as 
functions of f - 0)/2it > 0 rather than oo. Equation (B19) implies a 
relationship for such power spectral densities: 

S T\N< f > " 4 M f >| 2 S G ( f ), (B21) 

X 

where hjg(f) = hjj(to) and SQ(f ) is related to Sg(f) by a convolution [22] : 

S (f) = S (f ' ) S (f - f’) df'. (B22) 

o o 

The tilde (~) represents two-sided spectral densities with a property: 

S(-f) = S(f) = y S(f) , f > 0. (B23) 

For a sinusoidal function with random phase. 
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I 


g(t) = g cos (27if Q t + 9), 


(B24) 


it is straightforward to show [23] 

S g (f) - J g 2 t6(f " f Q ) + 6(f + f Q )], (B25) 

S G ( f) = y|- g 4 [6(f - 2f 0 ) + 46(f) + 6(f + 2f 0 )]. (B26) 

In practice, the 5-functions in these equations are replaced by continuous 
functions with finite peaks and widths due to the finite sample length t used 
in the computation of the spectral densities* Denoting this narrowband 
response function with a subscript r, the one-sided spectral densities can be 
written as 

S g (f) = i g 2 6 x (f - f Q ), (B27) 

S G (f) = i g 4 [26 x/2 (f) + 8 t/2 (f - 2f 0 )], (B28) 

where the line broadening by a factor 2 that takes place in the squaring 

operation [22] has been indicated by x/2. Notice that the factor 1/2 in front 
o 

of g comes from the assumption that the phase 9 of the signal is random. For 
a signal in phase with the reference signal, g 2 /2 must be replaced by g 2 in 
Eqs. (B27) and (B28). 

The exact line shape 6 T (f) depends on the window function employed in the 
Fourier analysis. A simple box-car shape window yields 

S T (f) = 8®<f) = e- j2 ’ £t dt - 8ln J’ lfT) . (629) 
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In the experiment, a Hanning window function [23] with the characteristic 


5 "< E > 5 i - 7> + 1 * i 5 '< £ + ^ < B30 > 

was used. Substitution of Eq. (B29) into Eq. (B30) leads to 

6,(f) = 6*1(0 = Sl o - ■ >- fT) T* (B31) 

T T 2nf 1 - (ft) 2 

_3 

Notice that the Hanning filter produces a tail that goes as f . 

In the frequency range 2 /t < f < f c , Eq. (B28) can now be approximated as 

4 

S G (f)3-S— -4. (B32) 

2 xt r 

r\ 

In the general case when Sg(f) is continuous, g^/2 must be replaced by an 
integral of Sg(f) over all frequencies so that 

S G (f) 2 -j [ ^o Sg(f’) df’] 2 . (B33) 

XT f J C b 

Substituting this into Eq. (B21), one obtains the final result: 

ISr,N«» 1/2 5 (|> 1/2 H p7T l h N (f >l C s g (f) df. (B34) 

Notice that Sp jj(f) can re duced by increasing the sample length t. 

The nonlinearity coefficient |hjj(f)| can be measured by driving the 
system with a sinusoidal acceleration at the frequency f/2: 

gp(t) = g F cos (xft + 9p) , (B35) 
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and observing the resulting gradient at Its second harmonic f: 


r s (t) = J gg cos (2mft + 0g). (B36) 

From the foregoing analysis, it is easily seen that 

g S = i IV f) l 4 • (B37) 

The value of | hj^( f ) | obtained from this relationship could then be substituted 
into Eq. (B34). 

3. Nonlinearities of the Instrument 

A general expression of the nonlinearity errors was given in Paper I: 

Sgd.N^) " h cc (w) G cc («) + h dd (w) G dd (w) + h cd (co) G cd (u)), (B38) 

1 where o cc M, Gdd(w) ^cd^) are t1ie Four i er transforms of [g c (t)]^, 

o 

an d gcCOgd^t)* derivation of the nonlinearity coefficients in 

| this equation requires knowledge of experimental details. 

We divide nonlinear mechanisms into three classes: 1) the instrument 

scale factor nonlinearity internal to the gradioraeter, 2) dynamical nonlinear 

i 

effects in which second order terms are produced by the finite compliance of 
( the gradioraeter structure in undesirable directions, and 3) nonlinear effects 
of the platform. The first class of nonlinearity can be overcome by 
linearizing the gradiometer by means of a "force rebalance" feedback; i.e. by 
operating the gradioraeter as a null detector in a feedback circuit. This 
method, however, will not suppress the latter two classes of nonlinearity 
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errors inasmuch as these nonlinear effects are produced external to the 


instrument. In the pendulum-suspended platform, the third class of non- 
linearity is produced by the centrifugal acceleration of the platform. In 
this section, we discuss the first two classes of nonlinearity, deferring the 
treatment of the third class to the following section. 

The scale factor nonlinearity can come from a failure of Hooke's law in 
the mechanical front end Ox/9g) , from a nonlinear response of the transducer 
(ai/&x), and from nonlinearity in the gain of the amplifier OV/5i) (see Fig. 

1 of Paper I). Since the SQUID operates in a negative feedback mode, it is a 
highly linear device and 9V/91 can be regarded as a constant. Further, the 
"cantilever” suspension spring for the proof masses has been specially de- 
signed to achieve a high degree of linearity along the sensitive axis. How- 
ever, the transfer function 9i/3x of the superconducting transducer is 
intrinsically nonlinear, as can be seen in Eqs. (19), (25), and (32) of Paper 
I. The second order terras in the transducer transfer functions vanish if L a - 
1^ (symmetric coils) and i = 0 (absence of symmetric current). The first 
condition, however, is difficult to satisfy in practice because the spacings 
between the coils and the surfaces of the proof masses are not easily matched. 

In order to obtain relationships between the nonlinearity coefficients of 
the component acceleration transducers and of the gradiometer, we expand the 
current responses i^( t ) of individual transducers as functions of driving 
accelerations g^(t): 

M'> ' h g t %<'> + h 01 g k (t) + •" 1 h gi Ig k< c > + °k g k (t) + , " 1 - (B39) 

Here an instantaneous response of the circuit has been assumed because of the 
low frequency nature of the signals under investigation. The coefficients of 
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the linear terms are matched by the common mode balance of the gradiometer so 
that the gradiometer output can be written as 


g* ( t) 


gd(t) + (C 2 - Cj) fg c (t) + ^ + (cf + C 2 N ) g c (t)g d (t) + •••. 


(B40) 


The transition from Eq. (B39) to Eq. (B40) is approximate because the 
electrical coupling between the two component transducers affects C^, as they 
are combined into a gradiometer. Taking the Fourier transform of this 
equation and comparing it with Eq. (B38), one can identity 


h cc (w) = Ah dd ( U ) 5 C 2 - , 


(B4la) 


M XT 

hcjCw) s Cj' + • 


(B41b) 


Notice that nonlinearity arises in the gradiometer even when the nonlinearitv 
coefficients of the individual acceleration transducers are matched: i.e., 

C 1 = The onl y wa y t0 eliminate the nonlinearity from the gradiometer 

scale factor is by making = 0. 

In the low frequency limit, one can derive an approximate relationship: 



1 A ^*ka Hcb 

Hip + 2L 0 Ha + Hb 


s ( 


1_ 

2 ^ 
“m d I 


L - 
a 

L + 




v k ’ 


(B42) 


where d L is the average coil-to-proof-mass spacing. Substituting values 
listed in Table 1, we find c” s - 7.5 x 10" 2 (m s -2 ) -1 and C 2 “ 

~1.7 x 10 1 (m s 2 ) - ^ so that 
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(B43a) 


h cc (u>) 5 = -1.0 x 10" 1 (m s“ 2 ) -1 , 

hcd( u ) - “2.5 x 10-1 g-2^-1^ (B43b) 

These values must be considered only as order-of-magnitude estimates in view 
of the approximations used in the derivation. 

The dynamic nonlinearity of the instrument arises from the finite com- 
pliance of the gradioraeter structure which allows the misalignment angles Sn_ 
and 5n A to be driven by external accelerations a(t) or a(t). From Eq. (A6b) 
of Paper I, it is clear that a nonlinearity error arises if 6n or 5n A is 
proportional to a(t) or a(t). A bending mode of the structure can be excited 
by the perpendicular component of a(t) or oc(t) applied to the midpoint between 
the two accelerometers. The symmetry of the induced motion is such that the 

A A 

sensitive axes of the component accelerometers, n^ and , go out of alignment 

A ^ 

with respect to the baseline vector 1, with opposite signs for a(t), and with 
an equal sign for a(t). Therefore, the induced misalignments have functional 
forms given by 


5n_(t - t') = h^(t - t') a^f), (B44a) 

aa 

5n At - t’) = h^(t - t') a ± (t ' ) , (B44b) 

+i aa 

where the history dependence of the coupling coefficients has resulted from 
the resonant response. Equation (A6b) of Paper I is then generalized into 

6g d N (t) = - [6n_(t-t ' ) • a(t') + 1 fin *(t-t’) x n • a(t')] dt'. (B45) 
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The general form (B18) is obtained if Eqs. (B44) are substituted into Eq. 
(B45). The corresponding equation for frequency domain is 


5g d>N (o>) = h^(w) G^(w) + h^(a>) G +> (w), (B46) 

’ aa aa oca aa 

2 2 

where G ) > (oj) and G +> (w) are the Fourier transforms of a^(t) and a^(t), 
aa aa 

respectively. 

The coefficient h ( ) (a)) has been computed [24] for a circular tube that 
aa 

connects the two accelerometers. We only quote the result here: 


K>>l " ~2 [l + 

aa iUg 



(_“)2jl/2 


f(f ), 

U B 


(B47) 


where 


f (r) = [(1 - r 2 ) + Q~ 2 r 2 ] 1 ^ 2 . 


(B48) 


Here u>g and Qg represent the (angular) resonance frequency and the quality 
factor of the bending mode involved. A similar result is expected for 

K>)|. 

aa 

Experimentally, the mounting cube for the gradioraeter was made rather 
rigid so that u)g/2 it 300 Hz. This gives upper limits for the nonlinearity 

error coefficients: 


|h^(u>)| = <,4.7 x 10" 6 (m s 2 ) -1 , « < Ug. (B49) 

aa aa 


This is completely negligible compared to the coefficients of the scale factor 
nonlinearity, evaluated in Eqs. (B43). Unlike the latter, however, the former 
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cannot be reduced by means of feedback operation of the instrument. It can 

therefore become a dominant source of nonlinearity when a very sensitive 

gradiometer is operated in the force rebalance mode. Further, larger values 

of h ? ) (to) and h > ^(a>) may arise under certain conditions from undesirable modes 
aa aa 

of the suspension springs of the proof masses which are much more compliant 

than the mounting structure of the gradiometer. 

Experimental values for the scale factor nonlinearity coefficients can be 

obtained in principle from the curvature of the calibration curve of the 

instrument. In the calibration data displayed in Figs. 7 and 12, however, the 

stochastic measurement error was too large to determine a systematic departure 

from a linear response. Only an upper limit for + dj can be obtained from 

—0 _0 

Fig. 7. For an applied acceleration g c 2 4 x 10 ‘ i s , the current response 
lies within 5 % from the linear fit of the data. Equation (B39) then implies 

T (C 1 + C 5 ) g c ~ °* 05 g c (B50) 


or 


cj + dj <; 0.1 g" 1 = 2.5 (m s" 2 )" 1 . (B51) 

The predicted value for dj* + dj, given by Eq. (B43b), is an order of magnitude 
smaller than this experimental limit. 

Using the nonlinearity coefficients predicted by Eqs. (43), one can com- 
pute the contribution of the scale factor nonlinearity to the low frequency 
noise spectrum of the gradiometer. The total common acceleration measured in 
the experiment is approximately 
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(B52) 


r S_(f') df' < 1CT 10 (m s“ 2 ) 2 . 

Substituting this and |h cc (f)| * 0.1 (m s“ 2 )"*, i - 0.16 m, x - 125 s into Eq. 
(B34), we find 

[S r )] 1/2 < 1 x 10" 2 E Hz" 1/2 ( 0,1 v H - Z -) 3/2 . (B53) 

Effects coming from terms involving g^Ct) can be shovm to be many orders of 
magnitude smaller* 

4. Centrifugal Acceleration 

Here we consider nonlinear effects of the platform which convert external 
accelerations into low frequency acceleration noise of the platform, which 
then couples into the by tb#> linear transfer function of the 

instrument. A well-known nonlinear mechanism of the platform is its angular 
motion driven by an external acceleration. The resulting centrifugal 
acceleration is a second order effect. 

An accelerometer mounted on the pendulum-suspended platform will ex- 
perience a centrifugal acceleration a 0 (t), in the radial direction away from 
the pivot point, with a magnitude: 

a c (t) - i p Qp(t). (B54) 


For sinusoidal motions : 


Qp(t) = Q p cos (u)t/2 + 0 q), 


(B55) 
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a C (t) - a c cos (ut + 0 a ) + a Co , 


(B56) 


the amplitudes are related by 


a C 


— Jl 

2 X P p 


(B57) 


By using Eqs. (A6b) and (A12), Qp can be related to the apparent acceleration 
gp measured by the accelerometer at the fundamental frequency f/2: 




|n • (Q x z) | 


X w 
P P 


§F> 


(B58a) 


where 9 = ^p/*9p is the unit vector along the direction of ^ and X Q = X^ has 
been assumed. Likewise, a^ can be related to the acceleration gg measured by 
the accelerometer at the second harmonic f : 


a C 


1 


A 



«s* 


(B58b) 


Substitution of Eqs. (B58) into Eq. (B57) leads to Eq. (B37) with the identi- 
fication : 


lV f )l 


4 


A A A O 

n • (Q x z) | g E 


(B59) 


where Eq. (Al) has been used. The quantity |h N (f ) | is the coefficient of 
apparent nonlinearity of an accelerometer suspended by a pendulum. 

When X Q A X for the accelerometer, Eq. (B58a) must be replaced by 


9 

P 


1 

A A A 

|n • (9 x z) | 


1 

(1 - X/X) n if g F» 
op p 


(B58a') 
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so that 


|h N ( f >l 


i • Z 1 1 

(Q x z) | 2 (1 - A /A ) 2 
o p 



(B59 ' ) 


Therefore, the magnitude of apparent nonlinearity coefficient depends on how 
well the horizontal acceleration is rejected by the pendulum action. 

For a gradiometer mounted on the same platform, Ap should be replaced by 
A in Eq. (B54) so that 


a r = -i- A Q 2 
^ L p 


(B60) 


With the aid of Eqs. (A6b) and (A19) , Q p and a c can be related to the apparent 
differential acceleration signals g F and g s at f/2 and f, respectively: 


Q. 


1 


|5n J 
+A 


itAf g F’ 


a. 


1 


1 - (n • Q) 


2 g S- 


(B61a) 


(B61b) 


where A Q = Ap is assumed so that the second term dominates in Eq. (A19). Sub- 
stituting Eqs. (B61) into Eq. (B60), we find 


i h (f) . . i - jl hi A A) 2 . 
In()I | 6 ;.i g E ^ 

+ A 


(B62) 


This is the coefficient of apparent nonlinearity of a gradiometer suspended by 
a pendulum. 

We now evaluate |hjj(f)| numerically for the gravity gradiometer in the 
umbrella orientation for three modes of operation: 1) common accelerometer 
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mode, 2) single accelerometer mode, and 3) gradiometer mode. For the particu- 
lar orientation chosen, one finds n • z - 1//3", n • Q - 0, n • (Q x z) » /2/3. 
The condition l Q = i p is satisfied for the common mode and the gradiometer 
mode. For the single accelerometer mode, one can show that 


|1 





(B63) 


Substituting 1 = 0.16 m, 
(B59), (B59 ' ) , and (B62), 


ip = 2.4 m, Q = 300, and |6n A | * 7 x 10“* into Eqs. 

+i 

we obtain the following theoretical predictions : 


O _ o — I 

8.0 x 10 J (m s~^)~ , common accelerometer, 

2f 

|h[j(f)|th = { 2.4 x 10 2 (ms 2 ) * (— single accelerometer, (B64) 

2f 

3.0 x 10^ (m s -2 ) -1 (— j^-) 2 , gradiometer. 


In the experiment, we were able to observe |h N (f)| at f - 2f p by de- 
tecting the second harmonic amplitudes for the swinging mode of resonance fre- 
quency fp = 0.32 Hz. Figure B2 shows the data for the common accelerometer 
mode (crosses), the single accelerometer mode (dots), and the gradiometer mode 
(triangles). The experimental values of |h N (2f p )| are determined from these 
graphs with the aid of Eq. (B37): 

1 o —0 i 

5.6 x 10° (m s z )~ , common accelerometer, 

4.6 x 10 2 (ra s -2 ) -1 , single accelerometer, (B65) 

4.4 x 10^ (m s -2 ) -1 , gradiometer. 
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5 - 


Fig. B2. Second harmonic generation by the centrifugal acceleration in 






Notice that all three of these values agree with the corresponding theoretical 
values within a factor of 2. 

This excellent agreement gives us confidence that the observed second 
harmonics are indeed the centrifugal acceleration of the platform. Since only 

a 

one parameter, |6n A |, is relatively unknown in Eq. (B62), one can use this 
equation and the measured value of |h N (2fp)| for the gradiometer to obtain a 
more reliable value of the misalignment : 


|6n A | = 6 x 1(T 4 . (B66) 

+A 

^ / 

This represents a better measurement of |6n A | than the estimate 7 x 10 - ^ 

+A 

obtained from the observed wideband noise spectra in Section VC. 

We now compute the low frequency noise of the gradiometer which is 
generated by the centrifugal acceleration. The formalism developed in Section 
B2 of this Appendix can be directly applied to the nonlinear equation: 


6g _(t) * A[1 - (n • Q) 2 ] Q 2 (t). (B67) 

d , C P 


A A O 

Substituting A[1 - (n • Q)^] 6(t - t’) into hjj(t - t') and Qp(t) into g(t), 
one finds the power spectral density of the centrifugal acceleration error: 


S (t) = [1 - (n • Q) 2 ] ~ [f S (f») df'] 2 , (B68) 

T,C kt f Q 


where S (f') is the power spectral density of ^ (t). 

Q p 

For the gradiometer under pendulum suspension, Eq. (A6b) leads to 


T S ff') df' 
o ± 


2nX 


h r 0 du 


-1 2 
T W 

P 


2 2 X 2 . -2 2 S * ( 2h* • (B69) 

a) A 


(<d - + T 

P P 
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where S^Cf) is the power spectral density of the horizontal acceleration 1^(0 

of the pivot point. Thus, the pendulum acts as a narrowband filter centered 

at fp so that S + (u>/2it) can be evaluated at fp and taken out of the integral. 

A 

The integral then has the well-known value it/2 so that 


r °y f ' )df '-^Y f ^ 


(B70) 


Substituting this and n • Q = 0 into Eq. (B68), we obtain 


s r;c (f> ^ ( V4[4wi 2 - 


f 3 V T‘ p> 


(B71) 


When the experimental values t p = 150 s, x » 125 s, A p = 2.4 m, and S^(f p ) 
7 x 10 -13 m 2 s -4 Hz~* are substituted, this yields 


[S (f)] 1 / 2 - 2 x 10“ 4 E Hz -1 / 2 ( P . tA— — )^/2 . (B72) 

r,c 

The derivation of this spectrum does not involve a detailed model for the 

apparent nonlinearity coefficient | h^j( f ) | . Neither does it depend on the 

estimate of integrated acceleration noise such as Eq. (B52), since the high 

frequency part of S (f) is cut off sufficiently fast by the pendulum. There- 

A 

fore, the resulting formula (B72) must give a reliable estimate of the low 
frequency noise generated by the centrifugal acceleration. 
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PART 3 
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SUPERCONDUCTING GRAVITY GRAD IOMETER 


- 181 - 


A SUPERCONDUCTING GRAVITY GRADICMETER FOR SPACE AND TERRESTRIAL APPLICATIONS 


M.V. Moody, H.A. Chan and H.J. Paik 
Department of Physics and Astronomy, 
University of Maryland, College Park, MD 20742 


A three-axis superconducting gravity gradiometer with a potential sensi- 
tivity better than 10“^ Eotvos Hz“ ly ^ 2 is currently under development for 
applications in space. Though such a high sensitivity may be needed for only 
a limited number of terrestrial applications, superconductivity offers many 
extraordinary effects which can be used to obtain a gravity gradiometer with 
other characteristics necessary for operation in a hostile moving-base envi- 
ronment. Utilizing a number of recently devised techniques which rely on 
certain properties of superconductors, we have produced a design for a sensi- 
tive yet rugged gravity gradiometer with a high degree of stability and a 
common mode rejection ratio greater than 10^. With a baseline of 0.11 m, a 
sensitivity of 0.1 Eotvos Hz -1 / 2 is expected in an environment monitored to a 
level of 10“ 2 m s“ 2 Hz“^ 2 for linear vibration and 7 x 10**^ rad s * Hz ^ 2 
for angular vibration. A conventional stabilitized platform can be used at 
this level. The intrinsic noise level, which is two orders of magnitude 
lower, could be achieved by monitoring the attitude with a superconducting 
angular accelerometer which is under development. In addition, the new 
gradiometer design has the versatility of adapting the instrument to different 
gravity biases by adjusting stored dc currents. 
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I. INTRODUCTION 


The Equivalence Principle of Einstein snakes it impossible, even in 
principle, to separate gravity and acceleration by a local measurement. How- 
ever, by making a differential measurement over a baseline, one can cancel out 
acceleration and detect gravity without being confused by platform motion. 
Although torsion balances have been used to detect gravitational force 
gradients for over two centuries, only in the most recent two decades have we 
seen serious efforts to develop moving— base gravity gradioraeters. ^ ^ Re- 
search on superconducting gravity gradiometers started more recently as an 
outgrowth of the superconducting transducer work for low temperature gravita- 
tional wave detectors.^ 

In a superconducting instrument, the inconvenience of cryogenic operation 
is offset by the opportunity of utilizing many exotic properties of supercon- 
ductors to improve Hie sensitivity and stability of gravity sensors. Iu addi- 
tion to the obvious reduction of the thermal noise of the instrument, the 
quantization of magnetic flux can be used to obtain "perfectly" stable means 
of signal transduction, scale factor matching, and proof mass levitation. The 
availability of SQUIDs (Superconducting Quantum Interference Devices) at 
liquid helium temperatures is another important factor that makes the super- 
conducting device attractive. SQUIDs are highly sensitive flux measuring 
devices which are based on the concepts of Josephson tunneling and fluxoid 
quantization in superconducting loops.® The commercial SQUIDs we employ are 
coupled to input coils for measuring small currents. The sensitivity of these 
instruments is 1.5 x 10 ^ a atu j the dynamic range is 10®. 

A three-axis superconducting gravity gradioraeter with a potential sensi- 
tivity better than 10 ® E Hz~^^ (1 E = 1 Eotvos = 10 - ^ s - ^) is currently 
under development at the University of Maryland.® This instrument has been 
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designed primarily for applications in space. However, a sensitive gravity 
gradiometer would also have a number of terrestrial applications if a satis- 
factory method of rejecting the high levels of environmental noise can be 
found. This paper describes a design for a superconducting gravity gradio- 
meter which incorporates several new features to help deal with the problems 
of a dynamically noisy environment. This design maintains a high sensitivity 
along with the convenience of a short base line. A versatile magnetic levita- 
tion is applied to the proof masses so that the same hardware can be operated 
in any gravity environment from 0 gg to 1 g£ (gg is the earth's gravitational 
acceleration) by adjusting persistent currents in a number of superconducting 
coils. 

Although the analysis in this paper will be confined to an in-line (i.e. 
a gradiometer which is sensitive to the diagonal components of the gravity 
gradient tensor, component gradiometer, it can be extended to a cross- 

component gradiometer (i.e. one which is sensitive to an off-diagonal 
component of the gravity gradient tensor, r^j, j t i). 


II. PRINCIPLE OF OPERATION 

An in-line component superconducting gravity gradiometer consists of a 
pair of spring-mass accelerometers coupled together by a superconducting 
circuit to measure differential acceleration.^ Each accelerometer consist of 
a superconducting proof mass confined to move along a single axis and a spiral 
superconducting sensing coil located near the surface of the proof mass (see 
Fig. 1). An acceleration will cause a displacement of the proof mass which, 
because of the Meissner effect, will modulate the inductance of the coil at 
frequencies down to dc. The sensing coil is connected to the input coil of a 
SQUID amplifier forming a closed superconducting loop. Since the flux in this 
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Fig. 1. Schematic diagram of a superconducting accelerometer. 



loop must remain constant, the change in the inductance of the sensing coil 
results in a current change through the SQUID input coil* In this manner very 
small accelerations can be detected* 

In the present design, each proof mass is confined to move along a single 
axis by a pair of low-loss cantilever spring systems* In practice, the 
dynamic axes of the proof masses cannot be perfectly aligned and this mis- 
alignment can cause various error signals to couple to the gradiometer output. 
These effects will be discussed in Section VI. As in any two-mass system with 
only one degree of freedom, the motions of the proof masses can be decomposed 
into a common mode (i.e. the displacements of the proof masses are in the same 
direction) and a differential mode (i.e. the displacements are in opposite 
directions). By coupling the two proof masses together by persistent currents 
I^jj and 1 ^ 2 * flowing in the closed superconducting loops shown by the solid 
lines in Fig. 2, and adjusting the ratio of 1^ to 1^2 > the sensitivity of the 
system to common-mode accelerations can be balanced out. By using a similar 

c o 

design, a balance of 2 parts in 10 has been demonstrated. Though this 
degree of balance should be sufficient in a low noise space environment, a 
higher degree of rejection to common-mode noise is necessary for terrestrial 
applications and ground tests of the instrument. 

By incorporating additional superconducting coils, shown by the dotted 
lines in Fig. 2, into the circuitry, the frequency of the common-mode 
resonance can be increased and the resonance peak passively damped without 
affecting the differential-mode resonance. This effect may be understood by 
noting that the flux in each of these loops, must remain constant. The 
electromagnetic energy in these two loops is given by, 


E - 


r. 

cl 

2(L 1 + L 2 ) 


* % 

c2 

2(L 3 + L 4 ) 


(i) 
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where and $ c 2 are the trapped fluxes. In a zero-g environment, = $ c 2 
is chosen. In the earth f s gravity environment, one of these fluxes can be 
greater than the other. When the gradioraeter experiences a common-mode 
acceleration, the two inductances in each loop change in a like manner re- 
sulting in a change in E. However, during a differential acceleration, the 
changes in the two inductances cancel and E remains constant. Increasing the 
common-mode frequency decreases the sensitivity of the gradioraeter to common- 
mode accelerations while making isolation of the common-mode resonance peak 
from environmental noise a simpler task. This isolation, along with the 
passive damping, limits the amplitude of the signal produced by the common- 
mode peak at the input of the SQUID amplifier and allows greater dynamic 
range . 

The upward shifting of the common-mode spring constant not only increases 
the rejection to common-mode accelerations, but it also increases the 
linearity of the gradiometer by confining the motions of the proof masses. 
Though the springs are designed for a high degree of linearity, the spring 
constant still contains higher order terms which may become significant for 
large displacements. 

A conventional approach to the linearity and dynamic range problem of an 
inertial instrument is to use an active feedback network which senses and 
cancels the response of the proof mass. One disadvantage of this approach is 
the possible introduction of a significant noise source from the feedback 
signal. Since the circuit discussed above is passive and superconducting, 
there are no additional noise sources to increase the fundamental noise level 
of the gradiometer. If necessary, active "force rebalance" feedback can of 
course be applied to both common and differential modes in addition to the 
simple passive circuitry. 
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III. DYNAMICS OF THE INSTRUMENT 


The inductance of a spiral coil located at a distance d from a supercon- 
ducting plane is given by 


L = P Q n 2 A d, 


( 2 ) 


where n is the turns density and A is the area of the coil. With damping 
ignored, the equations of motion for two proof masses, m^ and ra 2 , coupled by 
the superconducting circuits of Fig. 2 are 


X^t) + WjqIXjU) + x 1Q ] 


_L_ r— L2. x 2 (t) — — i 2 (t) + — — i 2 (t)] 
2m 1 1 cl U; d 3 c2 U ' d $ X dl u;j 


gl(t) 

(3a) 


x 2 (t) 


+ a) 20 [x 2 (t) + x 20 l 


1 

2m 2 



4 (t > 



i 2 c2 <t> 



i 2 2 (t» 


g2^>» 

(3b) 


where u> 10 and u) 20 are the uncoupled (angular) resonance frequencies of the 
proof masses and x^q and x 2 q are their equilibrium positions. The driving 
specific forces for the two proof masses are denoted by g^(t) and g 2 (t). 

Also, d^ is the equilibrium spacing between the i-th coil and the proof mass, 
and Lxq = (i () n 2 A i d i . 

The requirement that the flux in a closed superconducting loop must re- 
main constant imposes four constraints: 


{Ljq[1 + xj(t)/di ] + I^Q f 1 + x 2 (t)/d 2 ]} I c x(t) = $ cl , (4a) 

{LjqU - Xx(t)/d3] + L^gtl “ x 2 (t)/d4]} I c 2(t) ■ $ c 2> (4b) 
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{LjqI 1 + x 1 (t)/d 5 ] + Lg} I dl “ Lg I d2 - $ dl , 
{L 60 [i + x 2 (t)/d 6 ] + Lg} I d2 - Lg I dl “ $d 2 . 


The constraint equations give, to the first order in x/d. 


„ „ L 10 x l /d l + ho x 2 /d 2, 

*cl ~ ^lO^ 1 t — >y 


L 10 + L 20 


X c2 “ Z c20 (1 “ 


L 30 X l /d l * L 40 x 2^ d 4 
L 30 + L 40 


), 


T dl “ ^10 f 1 " 


T d2 ~ X d20 f 1 " 


L 50^ L 60 + x l^ d 5 + ^20^10^ 

(L 50 + L S ) (L 60 + L S ) " L S 
L 60^ L 50 + ^ x 2^ d 6 + ^10^20^ L 50 L S X l < ^ d 5 

<L 50 + V <L 60 + V - 


Upon substitution of these results into the equations of motion, 


finds 


+ 


*1 + sr i*io + *1 * *3 + *s» + r it )] 

1 S 

+ sj- [(K!K 2 ) 1/2 + (K 3 K 4 ) 1/2 + (K 5 K 6 ) 1/2 ] x 2 

_ _L_ 'hi t-2 _ ^30 2 L 50 2 

nij 10 2m j 1 dj clO d 3 c20 d $ dlO J 

•• 1 

*2 * iT *20 + *2 + *4 + *6< l + -r 11 x 2 

^ s 

+ ^ f(K l K 2> 1/2 + (K3 K 4> 1/2 + (K 5 K 6 ) 1/2 ] Xj 


«1» 


(4c) 

(4d) 

(5a) 
(5b) 
, (5c) 

. (5d) 

one 


(6a) 
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(6b) 


K 20 

+ -^ X 20 


1 f L 20 2 L 40 ,2 

20^ d 3 clO 


d, ‘c20 
4 


"60 t 2 

>6 


1 + d I d20^ ~ g 2» 


where 


K c = 


Ki ,-2 


2 2 
T cl0 L 10 


d l (L 10 + L 20 ) 


K 3 5 T2 


2 2 
1 c 20 1j 30 


d 3 (L 30 + L 40> 


2 2 
a50 L 50 


K 2 5 72 


2 2 
T T 

clO 20 


d 2< l l0 + So> 


K 4 = ~2 


2 2 
c20 b 40 


d ; a 30 + So> 


VSoSo'S + So + L 60> 


K 6 5 ~2 


4o L 60 


d 6 (L 50 L 60 /L S + L 50 + W 


Equations (6a) and (6b) can be rewritten in a simpler form: 


" V 1 v 3 2 

X 1 + ^-*l + 5p2 + “10 *10 - C 1 - 81 . 


v 2 v 3 2 

x 2 + ^ x 2 + ^ X 1 + u 20 X 20 " c 2 " 82* 


by defining 


Vj = K 10 + Kj + K 3 + K 5 (1 + L^/Lg), 

v 2 h k 20 + k 2 + k 4 + k 6 (i + l 50 /l s ), 
v 3 = (K^)^ 2 + (K 3 K a ) 1/2 + (K 5 K 6 ) 1/2 , 


and 


- 1 t 2 ^30 t 2 , L 50 t 2 >> 

C 1 = 2mj 1 dj clO " d 3 c20 d $ A dl0 ’ ’ 


(7a, b) 
(7c, d) 
(7e,f) 

(8a) 

(8b) 

(9a) 

(9b) 

(9c) 

(10a) 
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(10b) 


o = j_ (hi t 2 Ho .2 ^60 2 % 

2 " 2^ d 2 clO " d 4 c20 d 6 <120 J * 

Thus the equilibrium positions are given by 

x 10 " < C 1 + SlO^ /a) ?o » (lla) 

2 

x 20 = ( c 2 + 820 )/to 20 » (llb) 


where g^Q and g 2 g are constant bias forces. The solutions to Eqs. (8a) and 
(8b) are of the form: 


Xj^ = Ae iwt , x 2 = Be iu)t . 


( 12 ) 


Making these substitutions gives 


1,2 2 lm^ '-^m^ m^nu,-* 


(13) 


where the general solutions are 


iu^t -ia>it iw 9 t -iw 9 t 


xj = Aje + A_ 1 e + A 9 e + A_ 9 e 


(14a) 


iooit -itOit iwot -iuot 

x 2 = B^e + B_^e + B 2 e + B_ 2 e , 


(14b) 


As we will show later, if co^q and io 2 q are small, the amplifier noise is 
negligible and the coupling to the amplifier can be reduced. In this case 
and 1^2 are small and the terras involving K 5 and K 6 can be neglected. By 
combining Eq. (13) with Eq. (6a) or (6b), one can show that the eigenvalues 
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approach Xj + x 2 an <l x 2 ~ X 1 on ly if 


Kj K 2 K 3 K 4 

^ • 5 ^' ^ an,i “10 “ “ 20 - 


Thus in order for the coupled resonance frequencies to correspond to the true 
differential and common modes, the coil geometries of and l^, and L 3 and 
as well as the masses and the uncoupled resonance frequencies must be matched. 
Under these conditions, the differential-mode resonance frequency is given by 





), 


(16) 


and the common-mode resonance frequency is 


0 ) 


oc 


0 ) 


od 


m. 


(Ki + Ko). 


(17) 


f The eigenvalues are then 

f 

ito-jt “itlJ-jt 

X d = x 2 - X 1 = 2 ^ A d e + A -d e ° d (18a) 

1 iaK„t -iu rto t 

X c = y ( Xl + x 2 ) - 2(A c e oc + A_ c e oc ). (18b) 

The equations of motion can now be rewritten in the form of two simple 
harmonic oscillators : 


X d + “od (x d - x do> ■ «d. »«•) 

X c + “oc «c - x co> • *C. < 15 ») 
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where 


g d (t) 5 §2^) " 8i<t), (20a) j 

I 

g c (o s j ig^o + g 2 (t>]. ( 20 b) ; 

Since the displacement of a driven harmonic oscillator at frequencies below 

the resonance frequency is inversely proportional to the square of the 

resonance frequency, the sensitivity to common-mode acceleration is reduced by 
2 2 

a factor of a) /a) , . 

oc od 

The signal through the input coil of the SQUID is, from Eqs. (5c) and 

(5d), 


( 21 ) 

Thus the sensitivity to common-mode excitations may be further reduced by 
matching 1^20^6 anc * ^ act » this simple model, perfect common- 

mode rejection may be obtained in principle by adjusting 1^20 an< * 1^10* h° w “ 
ever, experimentally it is often easier to match several sets of parameters to 
moderate accuracy than to match one set to very high accuracy. 


t t - t - t + rr I<120 - y 

M2 Ml “ M20 M10 + d. ' X d 


-2 (i|10 + id20 ) ] [i ( J,^)]-^ 

d 5 d 6 L 50 L 60 


IV. SUPERCGNDUCING NEGATIVE SPRING 

A description of the superconducting negative spring has been presented 

Q 

in a previous paper . That paper, however, gave only a numerical solution. 

In this section, after a discussion of its application to the gradiometer, we 
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present an analytical solution which allows the data to be more easily related 
to other geometries. 

The noise power spectral density of the gradiometer can be expressed 10 as 

8 u 2 

S r (f) = — 2 l2nf k B T R(f) + E A (f)], (22) 

where ra, A, E A (f), and 0 ti are, respectively, the mass of each proof mass, the 
baseline, the amplifier noise energy (called the "input energy resolution"), 
and the energy coupling factor for the amplifier. The function R(f) is a 
f requency— dependent damping factor, which becomes equal to the inverse of the 
quality factor at the resonance frequency f * u> od /2-ru The first terra on the 
right hand side of Eq. (22) is due to the Brownian motion noise and the second 
Lenu oil the rignc is due to the noise ot the amplifier. This version of the 
sensitivity equation is different from the version which appeared in Ref. 

10. In its present form, the equation has been modified to include the fact 
that the magnitude of the force fluctuations at the signal frequency f is, in 
general, different from that at the resonance frequency. Namely, the Brownian 
motion noise has a frequency dependence which is governed by the nature of the 
loss mechanism in the spring. 

The superconducting gravity gradiometer at present has its fundamental 
noise limited by the amplifier noise rather than the Brownian motion noise. 8 
Equation (22) indicates that one of the most obvious ways to increase the 
resolution of the gradiometer is to lower its resonance frequency w Q( j. 

Lowering the mechanical spring constant while maintaining rigidity along the 
nonsensitive axes is a difficult task. One method of overcoming this dilemma 
for a superconducting gradiometer has been previously demonstrated.^ This 
method uses a superconducting negative spring to counteract the positive 
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mechanical spring. Each negative spring consists of a disk with a semi- 
circular edge located in a solenoid with a length less than the thickness of 
the disk. The proof mass is shaped to contain several of these "disks" (see 

! 

Fig. 2). The negative spring constant can be adjusted by changing the per- 

i 

sistent current in the solenoid. The lower limit for the resonance fre- 
quency will be determined by a number of factors. These factors include the 
required bandwidth of the gradiometer, the 1/f noise of the SQUID amplifier, 
and the magnitude of the higher order terras in the spring constant. 

When two large spring constants are balanced to obtain a low frequency 
spring, the stability and linearity requirements for each of these springs can 
become much more important. In the present design, the stability of the 
spring constants is maintained by a number of beneficial features which are 
available at liquid helium temperatures. These features include the stability 
of materials, the stability of persistent currents in superconducting loops, 
and the stable temperature environment. The linearity requirement is substan- 
tially decreased by the common-mode rejection coils which confine the motions 
of the proof masses. Also, if necessary, the stability and linearity of the 
system can be further enhanced by standard feedback techniques. 

To estimate the negative spring constant (see Fig. 3), we approximate 
that the field between the solenoid and the superconducting surface is cons- 
tant for a given displacement. Then, the change in the magnetic field energy 
with displacement is due to a change in the effective volume of the solenoid. 
The volume as a function of displacement x is 

V(x) = y itD [L(d 0 + R) - R 2 9] - V Q - j xDR 2 9, (23) 


where 
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(24) 
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r L/2 
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-1 r L/2 ~ 


"d + R 




Here, D, L, R and d Q are defined In Fig, 3. The magnetic potential energy Is 


E(x) 



o 


2L 


* 


2 

o 


2p Q n V(x) 


(25) 


Expanding E(x) In powers of x gives 


E(x) = j p Q n 2 I 2 (V Q - cjx 2 + c 2 x 4 + •••), 


(26) 


where 


xDR 2 L 


2(R + d ) 3 [1 + (L/2) 2 /(R + d ) 2 ] 2 ‘ 
o o 


(27) 


Dropping the higher order terms and differentiating twice with respect to x 
leads to the spring contant : 


5 2 E 

ax 2 


2 t 2 
- \i n I c 
o 


'-i • 

n 1 


(28) 


For the geometry in the aforementioned demonstration, this expression gives 
k ■ - 1100 N m - *, which is in good agreement with the experimental value of 
- 1180 N m -1 . 


▼. A PRACTICAL DESIGN 

If the effect of the resonance frequency is excluded, the determining 
quantities for the Brownian motion and amplifier noise terms of Eq. (22) are 
R(f) and E^(f), respectively. Presently, the most sensitive commercially 


- 198 - 


available SQUID 11 has an energy resolution of E A (f) = 3 x 10“^° J Hz -1 ^ 2 down 
to 0.1 Hz at which point 1/f noise becomes important. The damping factor of 
the proof mass motion contains contributions from both the mechanical and 
electromagnetic spring constants with the latter generally dominating for 
large coupling between the electrical and mechanical systems. Though dif- 
ficult to obtain, R(f) < 10 ^ has been observed in the superconducting coils 

l o 

of similar design near 1 kHz, For the purpose of our sensitivity calcula- 

tion, we assume that R(f) - 10 can be achieved at low frequencies (f 1 

Hz). 

With the above values in mind, and keeping with the objectives stated in 
the introduction, we propose the following parameters for a practical design: 
a differential-mode resonance frequency of 1.6 Hz, a baseline of 0.11 m, and a 
hollow niobium (Nb) proof mass 0.038 m in diameter by 0.029 m long with a mass 
of 0.1 kg. This design would give a sensitivity of 2 x 10~^ E Hz -1 ^. The 
dimensions can be reduced further if the sensitivity goal is set at a more 
moderate level of 0.1 E Hz -1 ^ 2 . 

For the common-mode rejection (CMR) coils and negative spring coils, the 
primary limitation is the critical field of the proof mass material. For 
Stanford grade niobium at 4.2 K, the critical field is 0.12 Wb m -2 . 1 ^ A 
second limitation is the minimum spacing between the superconducting coils and 
the surface of the proof mass. These coils have been made in the past with 
niobium wire. With a wire coil, the minimum spacing is about 1 x 10"^ m. One 
should be able to reduce this value substantially using thin-filra coils. 

For the CMR coils, the field value at the surface of the proof mass 
should be kept approximately 10% below the critical field value. This pre- 
caution would allow the gradiometer to withstand a common-mode acceleration of 
up to twice the earth's gravity without the field exceeding the critical value 
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and trapping flux in the superconductors. In order to maintain the correct 

spacing between the proof mass and the coils, two CMR circuits (containing 

persistent currents I c1 q and I c 20> respectively) with one coil on each side of 

2 

each proof mass are necessary (see Fig. 2). For vertical orientation, I^q " 
0.8 1^20 to compensate for the earth's gravitational field. Then, with 

p o nI c i0 “ °* 10 m ~ 2 > H o nI c20 " °* n Wb m ~ 2 » d = 1 x 10 ~ 4 m and A = 

1.1 x 10" 3 m 2 , and with Eqs. (17) and (7), the common-mode resonance frequency 
is 230 Hz. With this increased stiffness, Eqs. (19) and (21) imply that a 
total CMRR (common-mode rejection ratio) of 1 x 10^ is achievable, if the 
gradioraeter is balanced to two parts in 10^. 

For the negative spring coils, the field strength can be near the 
critical value giving nl n =* 0.9 x 10^ A m The expression in Eq. (27) has a 

maximum value of 0.27 at L/ (R + d Q ) - 2//3. If 12 disks with a diameter of 
3.8 x 10“ 2 m are used, one obtains k n = - 4.3 x 10^ N m which can 
compensate a mechanical resonance frequency of 33 Hz at 4.2 K. By lowering 
the temperature to 1.1 K, the critical field will increase by 25% allowing the 
mechanical resonance and the common-mode frequencies to be increased to 40 Hz 
and 300 Hz, respectively. 

VI. ENVIRONMENTAL NOISE AND INSTRUMENT ERRORS 

—3 —1 / 2 

In order to realize an operational sensitivity of 2 x 10 E Hz with 

a CMRR of 1 x 10^ and a baseline of 0.11 m, a linear acceleration noise level 

better than 2 x 10""^ gg Hz“^ 2 is required. The seismic noise level in a 

relatively quiet place is less than 10*^ g^ Hz~*^ 2 ; consequently, the passive 

CMR will be sufficient for a stationary platform. For a moving base applica- 

—3 

tion, however, the platform vibration level can be as high as 10 gg 
Hz”^ 2 . With this vibration level, a sensitivity of 0.1 E Hz ^ 2 would still 
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To extend a moving base system to 10 -3 E Hz - */ 2 , the platform vibrations would 
have to be monitored to a level of 10 -5 g E Hz -1 / 2 to compensate the common- 
mode errors of the gradiometer. A vector measurement of the platform 
acceleration with this resolution could be made by using a triad of conven- 
tional accelerometers. Alternatively, in a three-axis gradiometer, an addi- 
tional SQUID could be coupled to the (MR circuit of each component gradiometer 
to obtain a simultaneous reading of the three linear acceleration components 
of the gradiometer. 

With the common-mode error removed, the second most important error 
source is the angular motion of the gradiometer with respect to an inertial 
frame. Angular motion about an axis other than its own sensitive axis pro- 
duces o narfDof I \t al-fcrnorl cry-arM omp»f-£>r t*h TrM icrh rhp 

o * ' — r— j o u - - k _. 


centrifugal acceleration, which is indistinguishable from an in-line component 

gravity gradient. In order to suppress this error to a level of 2 x 10 -3 E 

Hz ' , for a vertical or horizontal orientation, the attitude rate of the 

gradiometer must be known or controlled to 1.4 x 10 - ** rad s - * Hz - */ 2 . For 0.1 

E Hz - */ 2 , this value becomes 7 x 10 - ^ rad s - * Hz - */ 2 . The requirement for the 
—1 12 

0.1 E Hz 1 instrument could be met with conventional gyroscopes mounted to 

the platform. The measurement of the attitude rate at the level of 10"® rad 
-1 - 1/2 

s Hz ' may be difficult for a conventional gyroscope. A superconducting 
"six-axis" accelerometer, which measures three linear and three angular 
acceleration components simultaneously with high sensitivity is under develop- 
ment^ and could be used for this purpose. 

Up to this point, this paper has dealt with a gravity gradiometer in 
which the sensitive axes of the component accelerometers are perfectly 
aligned. In a gradiometer whose sensitive axes are misaligned, linear and 
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angular motion will generate additional errors in the gradiometer output. 
Linear motion orthogonal to the direction along which a single-axis gradio- 
meter is balanced would couple directly to the gradiometer output at a level 
proportional to the degree of misalignment. Angular motion can couple in 
through axis misalignment in one of two ways. First, in the earth’s field, an 
angular displacement will result in a change in the dc bias level for the two 
accelerometers. When the two sensitive axes are misaligned with respect to 
each other, the change in the bias level will be different for the two 
accelerometers. Second, the misalignment of the average sensitive axis with 
respect to the baseline will result in a direct coupling of angular accelera- 
tion to the gradiometer output. We summarize below the error model associated 
with these misalignments. 

The gradiometer axis alignment errors can be described in terms of a 
misalignment between the sensitive axes of the component accelerometers : 


6n_ = n 2 - n 1 , (29) 


and a misalignment between the average direction of the sensitive axis and the 
direction of the base line: 


6n A = y (n ? + n. ) - l . (30) 

A. A 

In these equations, n^ and are the unit vectors in the direction of the 

A 

sensitive axes of the two component accelerometers and Z is the unit vector in 
the direction of the baseline. In addition to causing a gradiometer orienta- 
tion error, these alignment errors cause coupling to the gravity gradient out- 

A 

put from the common-mode acceleration component along the 5n_ direction and 
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from the angular acceleration component along the 5n A x n direction. 


+i 

In a terrestrial environment, a common-mode acceleration along the 6n 
direction is generated not only by linear motions, but also by angular motions 
which modulate the earth's gravity field g E « The error term along the Sn_ 
direction is then given by 


(f) = y 6n_ • [g x t (f) + a (f ) ] , 


(31) 


where $_(f) is the angular displacement noise and a (f) is the linear 
n n 

A A 

acceleration noise. The error term along the n A x n direction is given by 

+A 


61* A (f) = 5n A x n • a (f) , 
n + +i n 


(32) 


where a^(f) is the angular acceleration noise. A detailed derivation of Eqs. 
(31) and (32) will be published elsewhere. 

Using ordinary machining techniques and taking care to relieve stress in 

A A 

the mechanical components, the alignment errors Sn_ and 6n A can be reduced to 
the level of 10 . One possible method of improving the mechanical alignment 

is through the use of piezoelectric crystals. In such a system, a set of 
three or four piezoelectric crystal stacks would be used to adjust the rela- 
tive angle of the sensitive axes of the two accelerometers in a single-axis 

O A A 

gradiometer. An alignment of one part in 10° for both 6n_ and 6n A appears 

+A 

feasible by using this method. 

A 

A second method for reducing the alignment error 6n requires a three- 
axis gradiometer. In this method, additional superconducting circuits which 
are sensitive to the common-mode components of the acceleration along two axes 
are coupled to the proof masses of the third "orthogonal" axis. By adjusting 
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the persistent currents in these circuits, in a manner similar to the one- 
dimensional balance discussed in Sections II and III, the residual coupling 
between common-mode accelerations and the gravity gradient outputs due to axis 

i 

misalignment is balanced out. Thus a rigorous three-dimensional balance j 

I 

against sensitivity to linear accelerations is obtained. However, the angular 
acceleration error caused by the misalignment 6n A must be compensated for 
separately. 

Equations (31) and (32) determine the requirements for the attitude and 
attitude acceleration control/knowledge for a gravity gradiometer with a given ! 

i 

sensitivity. In order to achieve 2 x 10 -3 E Hz -1 / 2 , 9 n (f) and a n (f) must be 
controlled or known to 2 x 10 -6 rad Hz -1 ^ 2 and 2 x 10"^ rad s -2 Hz -1 ^ 2 , re- 

A A 

spectively, if 6n_ and 6n A can be reduced to the level of 10“ 8 . For 0.1 E 
Hz l ' , these requirements become 10~ 4 rad Hz** 1 / 2 and 10“ 2 rad s” 2 Hz -1 / 2 , re- 
spectively. Conventional gyroscopes could be used to satisfy these require- 
ments. If the alignment errors cannot be reduced sufficiently below 10”^, the 
gradiometer may be integrated with the superconducting six-axis accelerometer, 
which will have orders of magnitude improvement in attitude resolution over 
the conventional gyroscopes. 

We are also investigating a "pendulum suspension" 10 for the gravity 
gradiometer. If properly designed, a pendulum suspension can provide isola- 
tion in the three angular and two of the three linear degrees of freedom. 

Since the gradiometer can be balanced in the remaining linear degree of 
freedom, rejection of acceleration noise in all six degrees of freedom is 
achieved. Details of the pendulum isolation will be presented in a forth- 
coming paper. 

The extreme sensitivity of the gravity gradiometer requires careful iso- 
lation of the device from the thermal and electromagnetic fluctuations of the 
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environment as well. Below the larada point (T c = 2.17 K) , the liquid helium 
provides a stable and a gradient-free thermal environment. Also, the sensi- 
tivity of the instrument to temperature drift can be tuned out by employing a 
method similar to the common— mode acceleration rejection technique discussed 
in Section II. ^ The superconductor itself is a nearly perfect shield against 
fluctuating magnetic and electric fields. Thus, the superconducting gravity 
gradiometer can be isolated very effectively from the thermal and electro- 
magnetic noise, leaving the mechanical noise mechanisms discussed above as the 
most important error sources. 

VII. CGNCUJSIGNS 

The necessity of operating a very sensitive gravity gradiometer in an 
environment with a large common-mode acceleration background requires extreme 
stability in the acceleration-to-current transfer functions of component 
accelerometers and a reliable means of balancing out the common-mode sensiti- 
vity. The perfect stability of quantized magnetic flux in superconductors can 
be used to obtain a very sensitive gravity gradiometer with a high CMRR. Com- 
bining experiences obtained with a prototype superconducting gravity gradio- 
meter and new technological innovations, we have produced a design which gives 
a sensitivity of 2 x 10 - ^ E Hz - *^ an( j a oir i n excess of 1 x 10^. A three- 
axis in-line component gravity gradiometer, which incorporates many of the 
features discussed in this paper, is under construction for space applica- 
tions. This instrument has been designed for a relatively quiet environment 
and has been scaled up slightly to deliver a sensitivity of 10 - ^ E Hz - ^^ 

The new design utilizes magnetic levitation of the proof masses to null 
out the gravity bias, permitting operation of the instrument in an arbitrary 
orientation on the earth and in space. The low temperature environment gives 
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an opportunity to isolate the instrument from thermal and electromagnetic 
fluctuations in the survey vehicle. The inherent sensitivity of all gravity 
gradiometers to angular motion induced errors makes the attitude control of 
the gradiometer platform a challenge. However, superconducting techniques can 
again be employed to monitor the linear and angular motions of the platform 
with sufficiently high sensitivity and stability. The feedback and error com- 
pensation techniques which have been developed for conventional inertial 
navigation systems and gravity gradiometers could be adapted to the cryogenic 
instruments discussed here. 
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